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Analysis  of  the  Cost  of  Variable  Workloads 
on  Shipbuilding 

by 

A.  H.  Magnuson 
and 

R.  W.  Terry 

ABSTRACT 

The  effect  of  shipyard  workload  variation  on  the  cost  of  building 
ships  has  been  analyzed.  The  results  of  four  efforts  are  presented. 

The  first  major  effort  consists  of  an  analysis  of  the  effect  of  work 
density  (i.e.  worker  crowding)  on  shipbuilding  productivity  and  cost. 

The  results  show  that  an  optimum  least  cost  construction  time  and  work¬ 
force  level  exist  as  a  result  of  a  tradeoff  between  work  density  effects 
and  fixed  costs.  The  second  effort  was  an  attempt  to  identify  causes  of 
shipyard  productivity  variation  based  on  interviews  with  shipyard  super¬ 
visory  personnel.  The  third  effort  involved  development  of  a  framework 
for  estimating  transfer  functions  to  describe  how  workload  variation 
affects  cost.  This  work  is  to  be  based  on  historical  production  and  cost 
data.  A  description  of  adapting  the  Box-Jenkins  forecasting  methodology 
to  the  problem  is  presented.  The  fourth  effort  concerns  development  of 
a  shipyard  planning  system  to  minimize  cost  of  adjusting  to  workload 
variations.  A  review  of  current  approaches  to  multi -resource/multi -project 
planning  models  is  given  along  with  a  proposed  decomposition  of  the  planning 


problem  into  strategic  and  tactical  components.  The  strategic  or  long- 
range  planning  deals  with  aggregate  issues  such  as  organizational  goals, 
long-range  manpower  planning  and  facilities  expansion.  The  tactical  com¬ 
ponent  is  more  detailed  and  involves  workforce  allocation  on  a  trade  level 
to  the  various  activities  composing  the  construction  of  each  ship.  The 
tactical  planning  level  is  short-term  detailed  planning  that  takes  into 
account  inter-and  intra-trade  interferences,  precedence  relationship 
(proper  sequencing)  for  each  task  and  manpower  allocation. 
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Executive  Summary 


The  work  done  under  this  contract  is  described  in  detail  in  the  four 
Appendices,  each  of  which  is  a  self-contained  report  including  introductory 
discussion  material,  lists  of  references,  etc.  This  section  summarizes  the 
four  major  efforts  described  in  detail  in  each  appendix. 

The  first  effort,  Appendix  A,  deals  with  the  effect  of  work  density 
(worker  crowding)  on  shipyard  production  rate,  shipbuilding  cost  and  time. 

A  normalized  production  curve  is  developed  using  a  hyperbolic  tangent  function 
to  represent  the  dropoff  in  production  rate  as  work  density  effects  dominate. 
The  normal ization  is  accomplished  using  parameters  obtainable  from  historical 
data.  Total  cost  is  taken  as  the  sum  of  labor  costs  and  indirect  costs. 
Material  costs  are  not  taken  into  account  because  they  are  not  significantly 
affected  by  work  density  (workload).  A  minimum  cost  exists  because  of  the 
tradeoff  between  labor  and  fixed  costs.  The  minimum  cost  point  corresponds 
to  an  optimum  workforce  level  and  construction  time.  The  results  are  inter¬ 
preted  quantitatively  using  shipyard  data. 

The  second  major  effort  (Appendix  B)  is  entitled,  "A  Diagnosis  of  the 
Workload  Variation  Problem  in  Shipbuilding."  As  part  of  this  phase  of  the 
research,  a  telephone  interview  was  conducted  in  order  to  identify  the  primary 
preceived  causes  of  productivity  variation  in  shipyard  operations.  Several 
shipyards  were  also  visited  to  identify  specific  causes  of  variations  in 
productivity.  The  most  frequently  cited  reasons  for  productivity  vari¬ 
ations  were: 

1)  low  labor  progress, 

2)  low  manning  level , 
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3)  craft  incompatibility,  and 

4)  delayed  materials. 

The  first  three  of  these  causes  were  somewhat  expected  results  and  these 
topics  were  specifically  identified  on  the  interview  questionnaire. 

Material  shortages  is  an  item  that  was  not  specifically  listed  on  the 
questionnaire.  However  a  number  of  managers  listed  it  as  one  of  the  top 
three  causes  of  variations  in  productivity.  The  complete  results  and 
specific  questionnaire  is  available  in  the  report  along  with  an  analysis 
of  the  components  for  a  shipyard  planning  system. 

The  third  effort  (Appendix  C)  is  entitled,  "A  Framework  for  Analyzing 
How  Variations  in  Shipyard  Workloads  Impact  Shipbuilding  Cost."  This  report 
is  divided  into  three  principal  sections.  In  Chapter  II  of  the  report 
quantitative  definitions  of  workload  and  shipbuilding  costs  are  presented. 
Chapter  III  presents  the  causal  relationships  of  the  impact  of  workload 
variations  on  shipbuilding  cost.  The  costs  identified  include: 

1)  layoff  cost, 

2)  cost  of  hiring  new  workers, 

3)  rehiring  cost, 

4)  overtime  cost, 

5)  workloads  of  less  than  40  hours  per  week 

6)  subcontract  cost,  and 

7)  use  of  surplus  workers  to  perform  maintenance  work. 

In  Chapter  III  an  overview  of  the  Box-Oenkins  forecasting  methodology  is 
presented.  In  this  chapter,  a  transfer  function  model  is  introduced  and 
is  used  to  find  which  of  a  large  family  of  mathematical  models  best  de¬ 
scribe  how  a  change  in  input  affects  output  over  time.  This  methodology 
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is  developed  in  order  to  establish  a  set  of  difference  equations  that  can 
be  used  to  describe  how  workload  variations  dynamically  affect  shipbuild¬ 
ing  cost. 

The  fourth  report  (Appendix  D)  entitled,  "A  Shipyard  Planning  System," 
is  principally  concerned  with  the  development  of  a  framework  for  a  shipyard 
planning  system  which  shipyards  could  use  to  minimize  the  cost  of  adjusting 
to  workload  variations.  Chapter  I  provides  a  brief  introduction  to  the  plan¬ 
ning  problem  in  a  shipyard.  The  complexity  of  the  task  of  building  a  ship 
is  brought  out,  resulting  in  the  identification  of  the  system  as  a  multi¬ 
resource/multi-project  planning  system. 

In  Chapter  II,  a  brief  review  of  current  approaches  to  multi -resource/ 
multi-project  planning  models  is  given.  The  factors  identified  in  the  use 
of  these  models  are  (1)  the  data  requirements  which  requires  that  the  user 
specify  the  precedence  relationships  among  the  various  activities  in  building 
a  ship,  and  (2)  the  evolutionary  nature  of  the  shipbuilding  process.  The 
first  factor  results  in  a  combinatorial  problem  too  large  and  too  expensive 
for  any  single  heuristic  planning  model  to  handle.  The  second  factor  posses 
a  more  difficult  barrier  in  that  the  specification  of  the  sequence  in  the 
performance  of  some  important  construction  activity  cannot  often  be  identi¬ 
fied  a  priori  at  planning  time.  Such  precedence  relationships  tend  to  evolve 
over  time  as  the  work  on  the  various  ships  in  the  yard  progresses.  Thus, 
it  is  necessary  to  decompose  the  planning  prob'»m. 

Chapter  III  of  the  report  discusses  the  theoretical  framework  under  which 
the  Multi -Resource/Multi -Planning  Problem  can  be  decomposed  into  solvable 
parts.  The  planning  problem  is  broken  down  into  its  strategic  components 
and  its  tactical  components.  The  strategic  planning  problem  spans  a 


long-range  planning  horizon  and  deals  with  aggregate  issues  of  organizat¬ 
ional  goals,  long-range  manpower  planning,  facilities  expansions,  etc.  A 
clear  understanding  of  how  workforce  levels  and  project  completion  dates 
impact  cost  is  necessary  for  the  planning  performed  at  the  strategic  level. 
Such  costs  include  the  cost  of  late  deliveries,  materials,  and  labor  in¬ 
cluding  the  cost  of  hiring  and  firing.  The  strategic  planning  program, 
modeled  as  a  mixed  integer  quadratic  programming  problem,  is  discussed  in 
detail  in  Chapter  IV.  A  precedure  in  solving  this  problem  is  also  described 
in  this  chapter. 

The  tactical  planning  problem,  on  the  other  hand,  is  more  detailed  and 
would  involve  the  development  of  plans  on  how  each  trade's  workforce  should 
be  allocated  to  the  various  activities  necessary  for  the  construction  of 
each  ship.  The  planning  at  this  level  spans  a  shorter  planning  horizon 
and  requires  that  the  tactical  planner  recognize  that  (1)  both  intertrade 
and  intratrade  interference  can  cause  substantial  productivity  losses,  (2) 
failure  to  consider  strict  precedence  relationships  can  cause  unnecessary  de¬ 
lays  and  costs  to  be  incurred,  and  (3)  overtime  can  be  used  when  it  is  cost 
effective  to  do  so.  A  complication  that  is  inherent  to  the  tactical  planning 
process  is  the  random  variation  in  activity  duration  times  and  the  uncer¬ 
tainty  regarding  the  sets  of  activities  which  will  actually  be  necessary. 
Thus,  it  is  imperative  that  the  tactical  plan  be  revised  periodically.  The 
tactical  planning  problem  is  formulated  as  a  resource-constrained  cost-CPM 
type  of  problem  and  solve  via  a  procedure  based  on  the  model  developed  by 
Dar-El ,  et.  al .  Both  the  formulation  and  the  solution  procedure  are  dis¬ 
cussed  in  greater  detail  in  Chapter  V. 
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ABSTRACT 

A  quantitative  analysis  of  ship  construction  costs 
as  a  function  of  construction  time  and  workforce 
level  is  presented.  A  normalized  production  curve 
is  developed  using  a  hyperbolic  tangent  function 
to  represent  the  dropoff  in  production  rate  as 
work  density  (crowding,  worker  interference)  ef¬ 
fects  dominate.  The  normalization  is  accomplished 
using  parameters  obtainable  from  historical  data. 
Total  cost  is  shown  to  be  the  sum  of  labor  costs 
which  increase  with  work  density  and  indirect 
costs  proportional  to  construction  time.  A  mini¬ 
mum  cost  exists  as  a  result  of  the  tradeoff  be¬ 
tween  the  two  major  cost  groups.  The  minimum  cost 
is  associated  with  an  optimum  workforce  level  and 
construction  time.  The  results  are  interpreted 
quantitatively  using  shipyard  data. 

INTRODUCTION 

Any  meaningful  quantitative  analysis  of  the  ship 
construction  process  to  predict  costs  and  con¬ 
struction  times  becomes  very  difficult  because  of 
che  large  amount  of  factors  affecting  the  outcome. 
In  addition,  shipbuilding  is  an  activity  that  must 
be  differentiated  from  other  types  of  manufactur¬ 
ing  such  as  automobile  production.  The  large  body 
of  knowledge  on  mass-production  has  only  limited 
applicability  to  ship  construction  where  the  "pro¬ 
duction  runs"  are  restricted  to  tens  of  units  at 
most  (Reference  1).  Therefore,  the  analyst  must 
depend  primarily  on  his  own  resources  in  developing 
a  methodology  to  attack  a  specific  problem.  To 
make  the  present  problem  tractable,  the  limited 
scope  of  the  work  will  be  delineated,  terms  will 
be  defined  as  needed  and  simplifying  assumptions 
will  be  stated. 

This  paper  is  not  concerned  with  ship  design  or 
design  optimization.  The  starting  point  of  the 
analysis  is  the  following:  consider  a  ship  or 
series  of  identical  ships  to  be  built  in  a  given 
shipyard  where  the  construction  process  or  pro¬ 
cedure  and  facilities  are  specified  and  held  con¬ 
stant  throughout  the  construction  period.  The 
ship  design  is  fixed:  that  is,  the  configuration 
of  the  object  to  be  constructed  is  fully  defined 
by  suitable  drawings  to  the  necessary  level  of 
detail.  Once  the  design,  the  number  of  ships  to 
be  built,  the  construction  process  and  available 
facilities  are  known,  then  one  may  make  the  basic 
assumption  that  a  definite  cost  is  associated 
with  the  operation  (as  differentiated  from  the 
price).  The  construction  time  may  also  be  pre¬ 
dicted  assuming  the  shipyard  planners  schedule  the 
work  and  it  is  performed  subsequently  according 
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to  the  usual  work  pace  arrived  at  as  a  result  of 
experience.  This  construction  time  can  be  assumed 
to  be  the  nominal  optimum  construction  time.  (The 
construction  time  for  the  purposes  of  this  paper 
is  the  time  interval  between  contract  award  or 
other  authorization  to  commence  construction  to 
the  date  of  delivery  and  acceptance  by  Che  custom¬ 
er.) 

If  for  some  reason  the  normal  construction  process 
must  be  accelerated  or  slowed  down,  then  one  may 
expect  the  costs  to  Increase.  For  example,  the 
shipyard  planners  may  wish  to  stretch  out  Che  work, 
lengthening  the  construction  time,  to  avoid  laying 
off  workers  to  maintain  workforce  stability  or 
continuity  in  times  when  Che  overall  workload  is 
light.  Conversely,  the  construction  time  may  be 
shortened  when  the  yard  has  an  extensive  backlog 
of  work.  Accelerating  the  normal  work  pace  will 
require  an  increase  in  the  WQrkforce.  The  in¬ 
creased  workforce  will  not  be  able  to  work  as  ef¬ 
ficiently  because  of  limited  access  to  the  fixed 
production  facilities,  which  are  geared  to  che 
normal  workload.  In  addition,  simple  overcrowding 
may  occur  in  limited  workspaces,  particularly  in 
the  outfitting  phase.  The  overload  problem  then 
is  characterized  by  the  concept  of  work  density, 
a  term  first  used  by  Frisch  in  Reference  1.  Work 
density  (or  more  precisely,  worker  density)  can  be 
defined  as  Che  number  of  workers  per  unit  work 
area.  Most  if  not  all  phases  of  ship  construction 
are  performed  in  a  restricted  or  at  least  well 
defined  work  area.  Under  normal  workload  condi¬ 
tions  definite  values  can  be  assigned  to  the  work 
density  for  each  type  of  activity  such  as  pipe- 
fitting,  hull  fabrication  or  electrical  wiring. 

One  should  be  able  to  compute  a  global  (shipyard- 
wide)  mean  work  density  for  a  given  type  of 
ship  if  the  skill  or  trade  mix  is  fixed. 

If  the  work  period  is  extended  beyond  the  nominal 
optimum,  cost  will  increase  because  of  indirect 
or  overhead-related  costs.  These  fixed  costs, 
which  are  proportional  to  the  time  of  construction, 
Include  items  such  as  taxes,  depreciation,  utili¬ 
ties,  etc.  Indirect  costs  include  all  costs  that 
cannot  be  directly  charged  to  any  one  contract 
(Reference  2) .  It  is  clear  then  that  the  minimum 
cost  and  its  associated  optimum  construction  time 
are  both  the  consequence  of  the  tradeoff  between 
labor  costs  affected  by  work  density  effects  and 
fixed  carrying  costs.  The  purpose  of  this  paper 
is  to  analyze  this  relation  in  a  quantitative  way. 

Before  starting  with  the  analysis,  a  review  of 
the  usual  cost  estimating  procedure  used  by 
planners  or  Naval  Architects  is  instructive.  The 
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cost  estimation  procedure  used  in  bid  preparation 
and  design  studies  is  described  by  Benford  in 
Reference  2.  To  summarize  briefly,  labor  and 
material  costs  are  estimated  on  the  basis  of 
weights  (or  horsepower  in  the  case  of  propulsion 
machinery)  for  four  major  categories: 

1.  Steel  Hull 

2.  Outfitting 

3.  Hull  Engineering  (Deck  machinery  &  Equip¬ 
ment  ) 

4.  Machinery  (Main  Propulsion  and  Auxiliary) 

The  weights  are,  in  turn,  estimated  using  the 
cubic  number  and  other  parameters  describing  the 
snip  size,  configuration  and  proportions.  Benford 
(Reference  2)  includes  a  miscellaneous  direct 
labor  cost  to  cover  drafting,  purchasing,  schedul¬ 
ing,  material  handling,  etc.  This  is  computed  as 
a  percentage  of  the  direct  labor  costs.  Over¬ 
head  or  indirect  costs,  defined  as  Including  all 
costs  that  cannot  be  directly  charged  to  any  one 
contract  are  taken  as  a  percentage  of  the  direct 
plus  miscellaneous  labor  costs.  Benford  (Circa 
1960)  uses  an  overhead  rate  of  70Z  on  this  basis. 
Combining  the  above  four  functional  categories, 
three  major  cost  groups  can  be  identified: 

1.  Materials 

2.  Direct  labor 

3.  Overhead 


Mack-Forlist  (Reference  3)  gives  the  following 
percentage  breakdown  of  costs  for  three  types  of 
ships: 


Tanker 

Cargo  Ship 

Destroyer 

Material 

55 

55 

40 

Labor 

25 

23 

30 

Overhead 

20 

22 

30 

T'ne  overhead  rate  given  by  Mack-Forlist  is  higher 
than  Benford's  probably  because  Benford's  miscel¬ 
laneous  labor  cost  category  was  included  in  Mack- 
Forlist's  overhead.  The  exact  cost  breakdown  will 
vary  from  shipyard  to  shipyard  depending  upon  the 
accounting  system  used.  At  any  rate,  the  over¬ 
head  rate  at  a  given  yard  is  usually  given  as  a 
fixed  percentage  of  the  labor  costs.  This  rate  is 
based  on  an  average  workload  at  the  yard.  In  this 
study  a  more  general  view  of  overhead  or  indirect 
costs  will  be  used  that  explicitly  takes  into  ac¬ 
count  the  variable  workload.  One  would  expect 
that  in  an  underutilized  yard  fixed  expenses  would 
have  to  be  prorated  over  a  smaller  workload  thus 
increasing  the  overhead  rate.  Conversely,  an 
overloaded  yard  should  have  a  lower  overhead  rate. 

This  study  will  be  concerned  with  labor  costs  as 
affected  by  work  density  and  overhead  costs.  Ma¬ 
terial  costs  will  be  disregarded  as  they  are  not 
affected  by  workload  or  construction  time.  (Any 
carrying  charges  for  materials  can  always  be  in¬ 
cluded  in  the  overhead  category.) 

One  major  problem  in  an  analysis  of  this  kind  is 
how  to  define  the  governing  parameters  or  variables. 
The  variables  must  be  measurable  quantities  that 
have  some  clearly  understood  meaning.  In  addition, 
cause  and  effect  relationships  must  be  established. 


That  is,  one  must  distinguish  between  independent 
and  dependent  variables.  For  this  problem  the  in¬ 
dependent  variables  are  the  workforce  level  (W) 
and  the  construction  time  duration  (T).  These  are 
the  key  variables  the  planner  may  wish  to  manipu¬ 
late  to  compute  costs  when  stabilizing  the  work¬ 
force  or  to  expedite  work  in  an  overload  condition. 
The  workforce  level  is  taken  as  an  average  over 
the  construction  period  T  and  refers  to  production- 
related  workers  (those  not  charging  to  overhead). 
The  dependent  variables  (the  ones  the  planner  wants 
to  predict  as  a  function  of  W  and  T)  are: 

Cj  -  Total  Cost 

CL  -  Labor  Cost 

Cj  -  Indirect  Cost 

R  -  Production  Rate 

K  -  Unit  Labor  Cost 

Key  parameters  characterizing  the  ship  or  ships  to 
be  constructed,  the  production  facility  and  the 
construction  process  are: 

Ko  -  Low  work  density  unit  labor  cost 

R  -  Maximum  production  rate 

W°°  -  Critical  workforce  level  (“K  R  ) 
cr  o  * 

A  -  Work  area 

a  -  Critical  workforce  density  (-  W^/A) 

N  -  Number  of  units  (ships)  to  be  constructed 

W  -  Equivalent  workforce  level  for  indirect 
eq  costs 

B  -  Dimensionless  parameter  for  indirect 

costs  («  W  /W  ) 
eq  cr' 

T  -  Minimum  construction  time  (•  N/R  ) 

mm  ® 

C^p-  Zero  density  labor  cost  (»  NKq) 

Many  of  the  above  parameters  are  used  to  normalize 
the  dependent  and  independent  variables  by  forming 

dimensionless  ratios,  e.g.  W/W  ,  R/R  . 

cr  • 

NOND1MENSIONAL  PRODUCTION  RATE 

To  start,  consider  the  case  where  some  relatively 
simple  construction  task  is  being  performed.  For 
example,  fabrication  of  a  simple  subassembly  for 
a  ship  such  as  a  thrust  bearing  foundation  or  a 
support  pillar.  That  is,  some  activity  where  the 
labor  force  is  homogeneous  so  that  little  or  no 
specialization  is  required.  Later,  the  results  of 
the  analysis  will  be  applied  on  a  shipyard-wide 
basis  where  the  workforce  is  non-homogeneous .  On 
a  shipyard-wide  level,  the  workforce  skill  trade 
mix  for  each  yard  will  vary;  however,  each  yard 
will  have  a  characteristic  production  rate  reflect¬ 
ing  its  workforce  distribution.  In  addition,  the 
study  will  assume  that  the  workforce  is  already 
thoroughly  trained  so  chat  effects  of  learning  can 
be  ignored.  First  consider  the  case  where  the 
work  density  is  very  low,  i.e.,  where  there  is  no 
worker  interference  or  crowding.  The  input  is 
the  number  of  workers  (W) ,  say  in  worker-days  per 
day;  that  is,  the  workforce  level.  The  output  is 
the  number  of  units  produced  per  unit  time  (pro¬ 
duction  rate),  which  will  be  designated  R.  For 
example,  R  may  be  the  number  of  foundations  fabri¬ 
cated  per  day.  The  larger  in  mag  tude  or  more 


complex  the  construction  task,  the  lower  the  rate 
of  production.  For  a  given  construction  task,  the 
production  rate  is  proportional  to  the  workforce 
level  because  in  the  low  work  density  limit  there 
is  no  worker  interference.  For  example,  if  the 
workforce  is  doubled,  the  production  will  also 
double.  This  relationship  nay  be  expressed  as  a 
single  equation: 

R  -  W/K  ,  (1) 

o  o 

where  Kc  can  be  considered  an  index  of  the  cask 
magnitude  and  complexity.  This  equation  states 
that  the  production  rate  (R0)  is  proportional  to 
the  workforce  level  (U)  and  inversely  proportional 
to  the  complexity  (Kq).  Equation  (1)  can  be  re¬ 
arranged  to  solve  for  Kc: 

K  «  W/R  (2) 

o  o 

From  Equation  (2)  one  sees  that  the  complexity  K0 
has  the  units  of  W  divided  by  the  units  of  R.  For 
example  if  W  is  in  man-hours  per  day  and  R  is  in 
units  produced  per  day,  K0  will  have  units  of  man¬ 
hours  per  unit  produced.  That  is,  K0  is  the 
amount  of  labor  required  to  produce  one  item  in 
the  low-density  limit.  The  index  Kq  is  then  pro¬ 
portional  to  the  low  density  labor  cost. 

In  practice,  almost  all  tasks  are  performed  in  a 
more  or  less  confined  work  space,  where  worker 
interference  will  increase  as  the  number  of  work¬ 
ers  increases.  This  will  lower  the  increment  in 
productivity  as  W  increases.  This  effect  is  indi¬ 
cated  in  Figure  1  where  the  productivity  curves 
level  off  or  bend  downward  for  increasing  W.  This 
Figure,  taken  from  Ref.  1,  provided  the  basis  for 
this  study. 

h'ote  that  because  of  the  downward  curving  of  the 
lines  in  Figure  1  the  production  rate  no. longer 
doubles  if  the  number  of  workers  is  doubled. 

The  type  of  curve  shown  in  Figure  1  has  the  same 
general  shape  as  a  hyperbolic  tangent  (abbreviated 
tanh).  The  curves  in  Figure  1  resemble  hyperbolic 


tangents  because  they  start  out  linear  (straight 
lines)  for  very  small  values  of  the  argument  (W), 
and  eventually  level  off  to  a  final  constant  value 
for  large  values  of  Vi.  All  of  the  curves  in  Fig¬ 
ure  1  can  be  collapsed  into  a  single  hyperbolic 
cangent  curve  if  the  ordinate  (R)  and  the  abscissa 
(W)  are  suitably  normalised.  One  may  normalise 
the  ordinate  by  dividing  R  by  its  final  maximum 
(worker  saturated)  value,  denoted  as  R„.  Then  the 
final  value  of  the  normalized  production  rate 
(R/R.)  will  be  one.  One  may  normalize  the  work¬ 
force  level  W  by  dividing  by  Wcr,  the  critical 
number  of  workers.  The  parameter  Wcr  can  be  de¬ 
fined  (more  or  less  arbitrarily)  as  the  number  of 
workers  where  the  production  rate  reaches  76.2  per¬ 
cent  of  the  final  value  of  R,»,  as  Indicated  in 
Figure  2.  This  defines  Wcr  as  the  number  of  work¬ 
ers  where  the  argument  of  the  hyperbolic  tangent 
(W/Wcr)  is  equal  to  unity.  This  definition  of  WCr 
characterizes  the  falling  off  or  break  in  the 
hyperbolic  tangent  curve.  The  normalized  produc¬ 
tion  curve  is  shown  in  Figure  2.  This  curve,  in 
effect,  collapses  all  the  curves  shown  in  Figure  1. 
The  normalised  curve  in  Figure  2  is  written  in 
equation  form  as 

R/R_  -  tanh  (W/W^)  (3) 

The  maximum  value  of  the  production  rate  R„  can  be 
related  to  the  low  density  unit  labor  cost  Ko  and 
the  critical  workforce  level  Wcr  by  taking  Equation 
(3)  for  small  values  of  the  argument  (W/Wcr)  and 
setting  the  result  equal  to  Equation  (1).  Writing 
(3)  for  small  values  of  the  argument  gives  after 
rearranging: 

R  -  <R  J  «  <*) 

®  cr 

Comparing  (4)  with  (1)  indicates  that 


This  result  shows  that  Wcr  interrelates  the  low 
density  labor  cost  and  the  high  density  maximum 
production  rate. 


Figure  1.  Effect  of  Work  Density  on  Productivity 


n 
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One  naj  also  relate  K»  and  Ucr  to  measurable  param¬ 
eters.  These  parameters  are  the  low-density  com¬ 
plexity  or  cost  Ko,  the  work  area  and  the  critical 
work  density.  The  complexity  K0  has  been  defined 
earlier  in  Equation  (2).  The  work  density  will  be 
defined  more  precisely  below.  To  start,  consider 
R«, .  One  would  expect  the  maximum  or  final  value 
of  the  production  rate  to  decrease  with  the  com¬ 
plexity  Ko,  similar  to  the  way  it  did  in  Equation 
(1).  That  is,  one  would  expect  to  be  inversely 
proportional  to  K<>.  In  addition,  R«,  must  be 
proportional  to  the  available  work  area  A.  The 
work  area  is  an  easily  measured  quantity,  as  most 
construction  tasks  are  performed  in  a  restricted 
work  space.  Accordingly,  one  may  write  the  equa¬ 
tion  for  R,  as  follows: 

R^  -  oA/Kq,  (6) 

where  n  is  the  proportionality  constant.  The 
parameter  a  will  vary  depending  upon  the  construc¬ 
tion  task,  type  of  process,  equipment  used,  etc. 

One  may  relate  a  and  A  to  the  expression  for  Wcr 
in  (5)  using  (6)  to  obtain 

«cr  -  aA.  (7) 

Solving  Equation  (7)  for  the  parameter  a  gives 

a  -  Wcr/A.  (8) 

Equation  (8)  indicates  that  a  is  the  critical  work 
density  or  workforce  density,  l.e.,  the  critical 
number  of  workers  per  unit  of  work  area.  (The 
density  is  given  as  an  area  density  rather  than  a 
volume  density  as  it  seems  more  natural  in  this 
application.)  Equation  (8)  can  be  used  to  compute 
Wcr  in  (7),  as  the  critical  work  density  can  usual¬ 
ly  be  estimated  for  a  given  task.  For  example, 
in  the  case  of  non-mechanized  ditch  digging,  one 
would  expect  the  critical  work  density  to  be 
around  1.0  worker  per  square  meter,  based  on  the 
amount  of  space  needed  to  manipulate  a  pickaxe  and 
shovel  without  interfering  with  adjoining  workers. 

The  final  equation  characterizing  the  construction 
process  can  be  written  by  substituting  the  expres¬ 
sions  for  R„  (Equation  6)  and  Wcr  (Equation  7)  in¬ 
to  Equation  (3).  After  some  minor  rearranging. 


one  has 


R  -  (oA/K  )  tanh  (W/oA) .  (9) 

0 

Note  that  the  units  of  the  variables  in  Equation 
(9)  must  be  consistent.  For  example  if  the  units 
of  Ko  are  in  man-hours  per  unit  produced  and  the 
units  of  R  are  in  units  produced  per  day.  then  W 
must  be  in  man-hours  per  day. 

COST  ANALYSIS 

Referring  back  to  Equations  (1)  and  (2),  one  sees 
that  the  inverse  slope  of  the  productivity  curve 
gives  the  labor  required  to  produce  a  single  unit 
or  the  unit  labor  cost.  For  the  more  general 
case  where  the  work  density  is  taken  into  account, 
the  labor  coat  (denoted  K)  can  be  expressed  in  a 
form  similar  to  Equation  (2) 

K  -  W/R,  (10) 

where  R  is  given  by  Equation  (3) .  Substituting 
Equation  (3)  into  (10)  gives  after  some  manipula¬ 
tion 


K/K 

o 


(w/wcr) 

tanh  (W/W  )  • 
cr 


(11) 


One  sees  that  the  actual  labor  cost  relative  to 
the  low-denaity  cost  (K/K<>)  is  only  a  function  of 
the  ratio  W/Wcr.  Equation  (11)  is  plotted  in 
Figure  3.  One  sees  that  K/Ko  is  close  to  unity 
for  small  values  of  U/Wcr,  with  the  curve  increas¬ 
ing  to  1.313  for  W/WCT  -  1.  Below  0.5  the  effect 
of  work  density  is  fairly  weak.  At  the  critical 
workforce  level  the  labor  cost  is  31. 3!  higher 
than  the  low  density  level.  In  fact,  this  ob¬ 
servation  provides  an  alternative  method  of  defin¬ 
ing  the  critical  workforce  level. 

As  discussed  in  the  Introduction,  ship  construc¬ 
tion  costs  are  generally  broken  down  into  three 
major  categories:  direct  labor,  indirect  costs 
and  material  costs.  Since  material  costs  are 
not  affected  by  workforce  level  or  construction 
time  duration,  for  the  purposes  cf  this  study  they 
are  considered  to  be  constant,  and  will  be 
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DIMENSIONLESS  WORKFORCE  LEVEL  (W/W  ) 

cr 

Figure  3.  Relative  Labor  Cost  Versus  Workforce  Level 


disregarded.  The  simplified  cost  equation  may  then 
be  written  as  follows: 


C„ 


CL  +  CI’ 


(12) 


where  Cy  is  the  total  cost.  Cl  Is  the  direct  labor 
cost  and  Ci  are  indirect  costs.  One  may  write  the 
direct  labor  cost  for  constructing  N  identical 
units  (ships)  as 


CL  -  NK,  (13) 

where  K  is  the  unit  labor  cost.  One  may  also 
write  the  labor  cost  as 

CL  -  WT,  (14) 

since  the  labor  cost  is  proportional  to  the  con¬ 

struction  time  T  for  a  specified  workforce  level 
W.  Labor  cost  may  be  expressed  in  terms  of  work¬ 
force  level  using  (13)  as: 

C.  «  NK  (K/K  )  (15) 

L  o  o 


The  indirect  construction  costs  include  (for  the 
purposes  of  this  study)  all  costs  that  are  pro¬ 
portional  to  construction  time  duration.  As  men¬ 
tioned  in  the  Introduction,  indirect  costs  are 
usually  expressed  as  a  fixed  percentage  of  direct 
labor  costs.  Here  indirect  costs  are  expressed 
in  a  more  general  form  that  reflects  the  effect  of 
variable  workload.  Indirect  costs  also  reflect 
the  degree  of  capital  intensity  of  a  given  produc¬ 
tion  facility.  Indirect  costs  can  be  expressed  in 
terms  of  an  equivalent  workforce  level  in  a  form 
similar  to  Equation  (14): 


where  Weq  is  the  equivalent  workforce  level.  The 
total  cost  (excluding  material  costs)  can  be  writ¬ 
ten  from  (12),  (14)  and  (18)  as 

CT  -  (V  +  «)  T  (19) 

Fixed  costs  can  be  characterised  by  a  dimensionless 
parameter  6  as  follows: 


The  low  density  limit  of  (15)  is  from  (11): 


W  /U 

eq  cr 


(20) 


C.  «  NK  .  (16) 

Lo  o 

In  che  low  density  limit  one  sees  that  total  labor 
cost  is  independent  of  workforce  level.  The 
labor  cost  normalized  on  the  basis  of  the  low  den¬ 
sity  lir.it  is  then  (from  (15),  (16)  and  (ID): 


(W/W  ) 

C,  / C,  *  K/K  -  - — 7-  -Al — r  . 

L  Lo  o  tanh  (k/k  ) 

cr 


(17) 


Equation  (17)  indicates  that  the  normalized  total 
labor  cost  varies  the  s'.mc  as  the  unit  labor  cost 
as  would  be  expected  tor  the  same  workforce  level 


To  illustrate  the  variation  of  costs  with  construc¬ 
tion  time  T,  consider  first  what  happens  if  work 
density  effects  are  Ignored.  In  this  case  labor 
costs  are  constant,  and  equation  (12)  can  be  writ¬ 
ten  using  (18)  as 


W  2  +  C,  „ 

eq  o  LO 


(21) 


where  Clo  *  NK0  (from  16)  is  fixed  and  T0  is  the 
zero  density  construction  time.  If  total  cost  Cy 
is  plotted  against  construction  time,  one  sees 
that  the  costs  increases  linearly  with  time,  with 
the  minimum  cost  occuring  at  tine  T0  »  0.  This 
result  is  not  realistic  because  the  minimum  cost 
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can  only  be  achieved  wtch  an  Infinitely  large 
workforce.  This  can  be  seen  from  equation  (14) 
after  solving  for  W: 


K  *  CL0/To  (22) 

The  resulting  workforce  curve  is  a  nyperbola,  with 
V  going  to  infinity  as  T0  goes  to  zero.  It  is 
clear  that  the  minimum  cost  point  T0  *  0  cannot 
be  achieved  in  practice  because  work  density  ef¬ 
fects  will  decrease  the  productivity  of  the  work¬ 
force,  thus  increasing  labor  costs. 

To  obtain  a  more  realistic  characterization  of 
labor  costs  the  effect  of  diminished  productivity 
with  increasing  work  density  must  be  taken  into 
account.  For  a  production  run  of  N  units  (ships) 
the  production  rate  R  and  time  T  are  related  as 
follows: 


S  -  RT  (23) 

Recall  that  the  productivity  reaches  a  maximum 
denoted  R^  (see  Figure  2).  The  time  correspond¬ 
ing  to  maximum  production  rate  must  be  the  minimum 
denoted  Tnjn.  Since  N  is  specified,  one  has  from 

(23) 


T  .  -  N/R  .  (24) 

min  * 

That  is,  the  existence  of  a  maximum  production 
rate  requires  that  a  minimum  construction  time 
exists.  The  normalized  or  dimensionless  construc¬ 
tion  time  T/Tmin  can  be  obtained  from  (23)  and 
(24) : 

T/Tmin  *  <R/RJ_1  <25> 

From  equation  (3),  the  construction  time  is 


T/T  .  -  1/tanh  (V/W  )  (26) 

min  cr 

This  equation  relates  the  two  independent  vari¬ 
ables,  workforce  level  (W)  and  construction  time 
(T).  This  relation  (equation  26)  is  shown  in 
graphical  form  in  Figure  4.  One  sees  that  the 
workforce  level  is  infinite  at  T/Tmin  «  1,  i.e.  at 
a  finite  construction  time  (Tm^n)  instead  of  at 
zero  time  as  was  the  case  where  work  density  ef¬ 
fects  were  ignored. 


.Vow  returning  to  the  cost  equation  (19)  including 
work  density  effects,  one  has  after  normalization 

V<T»inUcr>  *  <»  +  W/Vcr>(T/W’  <27) 

where  (TminWcr)  is  a  normalizing  cost.  From  (24), 
(16)  and  (5)  one  has: 


T  .  W 
min  cr 


NK  “  C,_ 
o  LO 


(28) 


Equation  (27)  can  be  written  using  (28)  as 


CT/CLO-  (6+W/Ucr)(T/Tmin)  <29) 

One  may  eliminate  workforce  level  in  (29)  using 
(26)  which  gives 

VCtO  “{*  +  tanh'll(T/Tmin>‘1l}(T/T„in)'  <30> 

where  tanh  *  denotes  the  inverse  hyperbolic 


tangent.  Separating  the  total  cost  (29)  into  its 
components  gives 


VCL0  *  6<T/T«in> 

and 

CL/CL0  -  (tanh^KT/T^)-1))^^)  (32) 

The  two  components,  Cl  and  Cj  of  the  total  cost  Cp 
are  plotted  against  construction  time  in  Figure  (5). 
One  sees  that  labor  cost  is  no  longer  constant  as 
was  the  case  when  work  density  effects  were  ig¬ 
nored.  Instead  one  sees  that  labor  cost  Increases 
as  construction  time  decreases,  approaching  in¬ 
finity  when  the  time  equals  the  minimum  time  as¬ 
sociated  with  the  maximum  production  rate.  In¬ 
direct  costs,  varying  linearly  with  time,  are  shown 
for  three  values  of  8,  the  nondimens lonal  indirect 
cost  parameter.  Total  costs  obtained  by  adding 
the  components  in  Figure  (S)  are  plotted  in  Figure 
(6)  for  three  values  of  8.  One  sees  that  a  mini¬ 
mum  total  cost  exists  for  each  value  of  8. 

One  may  also  express  costs  as  a  function  of  work¬ 
force  level.  This  is  a  useful  representation  on  a 
shipyard-wide  level  where  one  may  want  to  compute 
the  optimum  workforce  level  needed  to  produce  a  given 
quantity  (N)  of  identical  ships.  The  corresponding 
optimum  construction  time  can  then  be  determined 
from  Figure  4  or  equation  (26).  Construction  time 
can  be  eliminated  in  (29)  by  using  (26),  giving 

8  +  (V/W  ) 

c  /c  ■ _ £I_ 

T/uLO  tanh  (U/W^)  ’ 

where  indirect  cost  is 


r  /r  m  - 2 _ 

V^LO  tanh(U/Wcr) 

and  labor  cost  is  given  by  (17). 

Total  costs  obtained  by  adding  fixed  and  labor 
costs  are  shown  in  Figure  7  for  three  values  of  8. 
It  is  apparent  that  there  exists  an  optimum  work¬ 
force  level  for  each  value  of  8  at  which  the  con¬ 
struction  cost  is  minimized.  The  minimum  cost 
workforce  level  can  be  computed  as  a  function  of  8 
by  differentiating  Equation  (33)  with  respect  to 
W/Ucr  and  setting  the  result  equal  to  zero.  The 
minimum  cost  construction  time  can  then  be  computed 
from  Equation  (26), 

RESULTS  AND  CONCLUSIONS 

In  Che  analysis  described  above  a  normalized  pro¬ 
duction  rate  curve  (Equation  3)  characterizing 
work  density  effects  was  derived.  Then  the  normal¬ 
ization  parameters  were  expressed  in  terms  of 
measurable  production-related  quantities.  Costs 
were  than  analyzed  using  normalized  construction 
time  and  normalized  workforce  level  as  independent 
variables.  Costs  were  broken  down  into  two  groups: 
labor  costs  and  indirect  costs.  The  results  are 
shown  graphically  In  Figures  6  and  7.  Total  cost 
has  a  definite  minimum  as  a  result  of  the  tradeoff 
between  labor  costs  and  indirect  costs.  The  op¬ 
timum  construction  time  and  workforce  level 


(33) 


(34) 
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correspond  to  the  minimus  cost.  The  optimum  point 
varies  with  the  dimensionless  parameter  6,  which 
characterizes  (but  is  not  the  same  as)  overhead 
rate. 

To  aid  in  interpreting  the  results  quantitatively, 
two  Tables  have  been  prepared.  Table  1  shows  the 
values  of  optimum  (minimum  cost)  time  and  work¬ 
force  level  for  the  three  values  of  6  indicated 
in  Figures  6  and  7.  Overhead  rate  was  computed  at 
the  optimum  point,  and  is  given  as  a  percent  of 
the  labor  costs.  This  permits  a  comparison  with 
Mack-For list’s  data  given  in  the  Introduction. 

Table  1.  Minimum  Cost  Values  for  Time  and  Work¬ 
force  Level  and  Corresponding  Over¬ 
head  Rate. 


£ 

T/T  , 
rnin 

w/w 

cr 

Overhead  Rate  (Percent) 

0.5 

1 . 45 

0.S7 

59 

1.0 

1.25 

1.06 

92 

1.5 

1.20 

1.19 

127 

Mack-Forlist 's  data  indicates  an  overhead  rate  of 
just  under  100  percent,  corresponding  to  the 
S  *  1.0  curve  if  one  assumes  the  yard  is  operating 
at  or  near  its  optimum  minimum  cost  point.  That 
is,  a  3  value  of  1.0  is  typical  for  U.S.  yards  ac¬ 
cording  to  Mack-Forlist' s  study.  It  should  be 
noted  that  the  optimum  construction  time  is  25 
percent  above  the  minimum  and  the  optimum  work¬ 
force  is  slightly  over  the  critical  workforce 
level.  This  means  that  work  density  effects  are 
quite  strong:  Figure  3  indicates  the  labor  cost 
is  33  percent  above  the  low  density  limiting  cost. 

The  curves  in  Figures  6  and  7  were  used  to  compute 
the  percentage  change  in  construction  time  and 
workforce  for  a  cost  increase  of  ten  percent  over 
the  minimum.  (Recall  that  cost  here  excludes 
material  cost.)  The  results  appear  in  Table  2. 

Consider  the  S  *  1.0  (92  percent  overhead)  case. 
The  increase  in  cost  can  be  achieved  by  either 
accelerating  or  slowing  construction.  The  first 
two  columns  slow  the  results  for  expediting  the 
work.  From  Table  2,  the  ten  percent  cost  increase 
will  reduce  construction  time  by  14  percent,  but 
will  require  a  31  percent  increase  in  production 
workers.  Alternatively,  the  construction  time  can 
be  increased  by  31  percent  with  a  ten  percent  in¬ 
crease  in  costs.  The  resulting  workforce  reduc¬ 
tion  in  this  case  is  40  percent.  It  is  apparent 
from  the  Table  that  it  is  more  expensive  to  ex¬ 
pedite  work  than  it  is  to  extend  it.  This  is 
because  labor  costs  blow  up  as  the  minimum  con¬ 
struction  time  is  approached.  (See  Figure  5.) 


In  addition,  note  that  as  B  (and  thus  overhead  rate) 
increases  the  optimum  minimum  cost  point  gets 
closer  to  Tnin-  This  means  that  yards  with  higher 
overhead  rates  reflecting  higheT  capital  intensity 
tend  to  be  less  flexible  in  adjusting  to  changes  in 
workload.  That  is,  the  incremental  costs  are 
higher  for  the  same  change  in  construction  time. 
This  effect  will  be  alleviated  to  some  extent  by 
the  lower  unit  costs  that  should  result  from  a 
higher  degree  of  automation.  Recall  that  costs 
have  been  normalized  using  low  density  labor  costs, 
which  should  decrease  as  automation  (capital  in¬ 
tensity)  Increases. 

Note  that  by  assuming  that  a  given  yard  is  operat¬ 
ing  at  a  particular  time  at  its  minimum  cost  point, 
a  particular  value  of  B  can  be  assigned  since 
the  overhead  rate  is  known  (Table  1).  This  then 
determines  the  yard's  cost  curve  as  a  function  of 
workforce  level  or  construction  time  (Figures  6  and 
7).  Then  the  workforce  can  be  determined  relative 
to  the  critical  workforce  from  Table  1  as  the  ratio 
(W/Wcr) .  The  following  identity  can  be  developed 
using  Equation  (7): 

W/Wcr  -  (W/A)/(Wcr/A) 

-  (VI/A )/a 


Since  (W/Wcr)  and  the  total  work  area  of  the  yard 
are  known,  the  yard's  critical  work  density  a  can 
be  found.  Similarly,  the  low  density  labor  cost 
Kq  can  be  determined  from  Figure  3  since  W/Wcr 
and  labor  cost  K  are  known.  That  is,  one  may  work 
backward  from  the  assumed  optimum  point  to  obtain 
quantitative  values  for  the  parameters  needed  in 
the  analysis. 

The  alternative  procedure  is  to  start  with  a  pro¬ 
duction  curve,  i.e.,  production  rate  versus  work¬ 
force  level.  This  curve  can  be  curve-fit  to 
Equation  (3),  the  normalized  production" curve  pro¬ 
vided  the  data  cover  a  sufficient  range  to  include 
the  dropoff  in  output  due  to  work  density  effects. 
Once  the  production  curve  has  been  established  the 
other  parameters  are  relatively  easy  to  obtain  from 
historical  data.  Obtaining  the  production  curve 
from  historical  data  is  especially  difficult  be¬ 
cause  the  data  is  "contaminated"  by  numerous  ef¬ 
fects  which  have  not  been  taken  into  account  in 
this  study.  These  effects  include: 

1.  Varying  types  of  ships  (Nonuniformity  of 
items  being  constructed) 

2.  Inflation 

3.  Learning 

4.  Changes  in  yard's  production  facilities/ 
processes 

5.  Variation  of  workforce  over  construction 
period. 


Table  2.  Increments  in  Construction  Time  and  Workforce 
Level  for  a  Ten  Percent  Increase  in  Costs. 


g/(Fercent 

Overhead) 


Percent 
Reduction 
in  Constr.  Time 


Percent  Increase 
in  Production  Workers 


Percent 
Increase 
in  Constr.  Time 


Percent  Reduction 
in  Production 
_ Workforce 


0. 5/(59) 
1.0/(921 
1.5/(127) 


22 

14 

13 


47 

31 

32 


47 

31 

32 


39 

40 
45 


Further  work  is  planned  to  verify  the  nodel  pre¬ 
sented  here  using  historical  data  after  compen¬ 
sating  for  the  factors  listed  above.  This  phase 
■••ill  require  extensive  computer  analysis. 
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PREFACE 


The  Industrial  Engineering  and  Operations  Research  Department's 
activities  were  divided  into  three  phases.  The  first  phase,  which  is 
reported  in  this  document  was  designed  to  identify  specific  problems  for 
further  analysis.  This  analysis  revealed  the  need  for  the  following: 

(1)  a  system  for  analyzing  how  variations  in  shipyard  workloads 
affect  shipbuilding  cost; 

(2)  a  shipyard  planning  system  which  could  be  utilized  by  shipyards 
to  minimize  the  cost  of  adjusting  to  workload  variations. 

Phase  2  focused  on  developing  the  framework  of  a  system  for  developing 
transfer  function  models  which  can  be  used  to  develop  the  difference 
equations  which  describe  how  variations  in  shipyard  workloads  affect 
shipbuilding  cost.  Phase  3  focused  on  developing  the  framework  of  a 
shipyard  planning  system  which  shipyards  could  use  to  minimize  the  cost 
of  adjusting  to  workload  variations. 


OBJECTIVES 


The  purpose  of  this  phase  was  to  determine  the  relative  importance 
of  the  problems  created  by  workload  variation  and  to  identify  specific 
problems  for  further  analysis  in  the  second  and  third  phases.  This 
research  has  focused  on  two  objectives.  The  first  objective  was  to 
develop  the  framework  of  a  system  for  analyzing  how  variations  in 
shipyard  workloads  affects  shipbuilding  cost.  Such  analyses  could  be 
used  by  those  responsible  for  establishing  Navy  shipbuilding  budgets  to 
analyze  how  alternative  funding  programs  influence  both  procurement 
cost  and  mobilization  potential.  The  second  objective  was  to  develop 
the  framework  for  a  shipyard  planning  system.  Such  a  system  could  be 
utilized  by  shipyards  to  minimize  the  cost  of  adjusting  to  workload 
variations. 

LITERATURE  REVIEW  AND  SHIPYARD  VISITS 

A  variety  of  techniques  were  utilized  to  acquire  this  information. 
Discussions  with  personnel  in  the  Naval  Sea  Systems  Command,  the  Office 
of  Naval  Research,  the  Maritime  Administration,  the  Census  Bureau,  and 
the  Bureau  of  Labor  Statistics  uncovered  a  large  number  of  potentially 
relevant  documents  (see  Appendix  A) .  These  documents  were  reviewed  for 
relevance  and  a  classification  system  (see  Appendix  B)  was  developed  to 
facilitate  locating  potentially  useful  information.  The  following  ship- 

m 

yards  were  also  visited: 

(1)  Norfolk  Naval  Shipyard;  Portsmouth,  Virginia;  December  21,  1978 

(2)  Norfolk  Shipbuilding  and  Drydock  Company;  Norfolk,  Virginia; 


March  20,  1979 
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(3)  Newport  News  Shipbuilding  and  Drydock  Company;  Newport  News, 
Virginia;  March  20,  1979 

(4)  Todd  Shipbuilding;  Seattle,  Washington;  February  2,  1979. 

The  documents  reviewed  and  the  confidential  discussions  with  shipyard 
executives  during  the  visits  revealed  that  the  following  were  potential 
causes  of  low  productivity: 

Unavailability  of  labor  at  the  time  needed  (low  labor  progress) ; 
Insufficient  skilled  personnel  (low  manning  level); 

Improper  order  of  sequence  of  assigned  labor; 

Interference  by  different  crafts  (incompatibility); 

Too  many  personnel  working  at  one  time  (work  density) ;  and 
Accelerated  schedules  (work  speedup) . 

STRUCTURED  TELEPHONE  INTERVIEWS 

In  order  to  determine  the  extent  to  which  these  potential  causes  actually 
occurred,  a  structured  interview  questionnaire  (see  Appendix  C)  was  design 
ed  to  be  administered  by  telephone  to  executives,  industrial  engineers, 
and  production  or  construction  scheduling  managers  at  various  shipyards. 
The  shipyards  were  chosen  so  that  a  representative  sample  of  large  and 
small  shipyards  in  several  Naval  districts  would  be  included  in  the  survey 
Twenty-one  shipyards  were  included  in  the  interviews  (see  Appendix  D) 
All  but  two  had  facilities  to  construct  ships  or  barges  larger  than  300 
feet,  and  all  but  four  had  employment  of  over  1,000  personnel.  Eleven  of 
the  shipyards  weie  categorized  as  "majors";  and  seven  Naval  districts  were 
included  overall.  Seven  of  the  shipyards  are  currently  constructing  Navy 
ships,  and  all  but  two  were  involved  in  commercial  ship  construction. 

20 

' . .  . . 


ANALYSIS  OF  STRUCTURED  TELEPHONE  INTERVIEWS 


The  results  of  these  structured  telephone  interviews  are  summarized 
in  Table  I.  The  entries  in  this  table  indicate  the  number  of  shipyards 
surveyed  that  ranked  the  causes  of  variation  in  productivity  in  the  top 
three.  Shipyards  were  classified  with  respect  to  size  and  whether  or  not 
they  built  Navy  ships. 

The  most  striking  feature  of  this  table  is  the  high  number  of 
responses  pertaining  to  low  labor  progress  and  low  manning  level  at  major 
shipyards,  particularly  those  performing  Navy  and  merchant  ship  construction. 
This  suggests  the  importance  of  labor  and  craft  availability  at  many 
shipyards  constructing  merchant  and  Navy  ships. 

Intermediate  size  shipyards  were  more  likely  to  indicate  problems 
related  to  improper  order  of  assigned  personnel,  craft  interference,  and 
work  density  as  the  major  causes  of  variation  in  productivity.  Smaller 
shipyards  selected  low  manning  levels,  craft  incompatibility,  work  density 
and  accelerated  schedules  as  major  causes  of  productivity  variation.  In 
the  category  cited  as  "other  causes",  a  number  of  shipyards  reported 
material  delays  and  shortages  to  be  of  major  concern.  Each  of  these  items 
is  discussed  in  more  detail  below. 

Taken  together,  it  is  evident  that  unavailability  of  labor  and 
insufficient  skilled  personnel  are  regarded  as  major  causes  of  variation 
in  productivity  by  shipyard  executives.  This  finding  corresponds  with  a 
1977,  Department  of  Defense  report  which  stated  that  the  majority  of  all 
completion  delays  were  due  to  low  labor  progress  resulting  from  inadequate 
manning  by  the  shipbuilders  (Department  of  Defense,  1977). 


Primary*  Causes  of  Variation  in  Productivity  Reported  by  Shipbuildini 
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The  survey  revealed  that  intercraft  interference  and  intracraft  work 
density  are  significant  problems  as  well.  Several  managers  pointed  out 
that  these  conditions  can  be  brought  about  by  accelerated  schedules  and 
improper  assignement  of  personnel  even  though  shipyard  manning  is  at  a 
relatively  low  level.  Some  find  it  necessary  to  stagger  personnel  in 
shifts  in  order  to  reduce  craft  interference.  (This  is  particularly  true 
on  smaller  vessels,  such  as  submarines.) 

Intercraft  interference  and  work  density  are  most  often  experienced 
in  the  engineering  spaces  aboard  ship  and  in  some  of  the  more  specialized 
compartments  during  outfitting.  Navy  ships,  with  their  complex 
electronics  and  weapons  systems,  often  experience  this  concentrated 
activity.  The  crafts  which  were  found  to  be  least  compatible  due  to  the 
nature  of  their  work  are  painters  and  finishers.  Also,  welders  can  have 
no  personnel  working  beneath  them,  and  forced  air  blowers  often  restrict 
the  amount  of  work  other  crafts  can  perform  at  the  same  time. 

The  problem  of  craft  interference  and  work  density  may  be  related  to 
a  problem  cited  by  shipyard  workers  in  a  MarAd*  study  of  nearly  1,300 
employees  representing  ten  shipyards.  The  most  common  spontaneous  com¬ 
plaint  among  production  workers  which  is  related  to  working  conditions 
concerned  inadequate  scheduling,  planning,  coordinating,  and  communicating 
between  crafts,  shifts  and  various  working  groups  (Muench,  1976). 

The  survey  also  revealed  a  problem  which  appears  to  be  growing  in 
dimension  at  many  of  the  shipyards.  A  number  of  managers  listed  material 
shortages  and  delays  as  being  among  the  top  three  causes  Of  variation  in 
productivity  even  though  late  material  was  not  listed  as  a  specific  item 


*  Maritime  Administration,  Washington,  D.  C. 
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on  the  questionnaire.  Late  material  deliveries  were  discussed  in  a  1977 
Department  of  Defense  report,  but  no  increase  in  lead  time  trends  was 
apparent  at  the  time  (Department  of  Defense,  1977).  Subsequent  reports 
may  produce  different  results.  A  number  of  shipyard  managers  believe 
the  problem  to  be  getting  worse. 

Another  outcome  of  the  shipyard  survey  summarized  in  Table  II  is  an 
illustration  of  the  degree  of  activity  in  ship  construction  and  the  stages 
at  which  the  greatest  concentration  of  personnel  is  likely  to  occur.  For 
most  ships,  this  occurs  during  machine  installation  and  outfitting.  A 
notable  exception  is  the  construction  of  tankers  in  which  the  most 
activity  occurs  during  the  steel  fabrication  phase.  This  is  particularly 
true  of  shipyards  using  modular  construction  or  block  techniques  which 
bring  needed  sections  of  the  hull  together  as  they  are  fabricated  and 
joined . 

NEED  FOR  SHIPYARD  PLANNING  SYSTEM 

The  results  of  the  above  survey  clearly  indicate  that  an  effective 
shipyard  planning  system  must  take  the  following  into  account: 

(1)  Labor  availability 

(2)  Material  availability 

(3)  Craft  incompatibility  and  interference. 

Failure  to  do  so  often  results  in  the  scenario  described  by  Chirillo  as 
"the  falling  domino  effect"  illustrated  in  Figure  1.  In  this  scenario, 
low  labor  progress  on  precedent  activities,  delayed  materials,  or  workflow 
disruptions  (caused  by  inclement  weather,  strikes,  etc.)  bring  about  a 
period  of  forced  idleness  and  low  labor  progress  which  results  in  an 
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Liquified  Natural  Gas  Carrier 
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and  Disruptions  in  Shipbuilding 


attempted  work  speedup  during  a  later  period.  When  the  schedule  is 


accelerated  to  make  up  for  lost  time  and  avoid  late  penalties,  inefficien¬ 
cies  occur  due  to  improper  sequence  of  assigned  labor,  interference  between 
crafts,  and  interference  within  crafts  as  too  many  personnel  require 
access  to  the  same  space.  The  review  of  documents  and  discussions  also 
revealed  a  number  of  statements  which  describe,  in  qualitative  terms,  how 
variations  over  time,  in  the  demand  for  ships,  adversely  impacts  cost  and 
productivity . 

The  following  statements  were  found  in  a  report  prepared  for  the 
Maritime  Administration. 

Introduction  of  new  or  improved  management  procedures 
[should  be  encouraged],  particularly  those  which  improve  ship 
construction  scheduling  and  thus  minimize  the  present  queueing 
problems  that  occur  when  certain  shipyard  occupations  are  used 
intensely  for  some  time,  only  to  stand  idle  later.  The 
lessening  of  insecurity  about  this  layof f /rehire  pattern  could 
attract  more  workers.  (Mark  Battle  Associates,  Inc.,  1974, 
p.  8) 


Generally,  the  rate  of  voluntary  quits  and  discharges 
increases  significantly  during  shipyard  employment  expansion. 
(Mark  Battle  Associates,  Inc.,  1974,  p.  45) 

Higher  turnover  rates  were  reported  by  shipyard  officials 
to  be  highly  correlated  with  manpower  expansion  in  the  ship¬ 
yards.  (Mark  Battle  Associates,  Inc.,  1974,  p.  55) 


NEED  FOR  MODEL  TO  PREDICT  HOW  VARIATIONS  IN  DEMAND  AFFECT  SHIPBUILDING  COST 


In  addition  to  the  above,  Mr.  E.  Karlson  of  the  Maritime  Administration, 
provided  the  following  statement: 

"Shipyard  productivity  is  difficult  to  measure  because 
of  the  uniqueness,  sophistication  and  diversity  of  vessels 
being  constructed.  There  are,  however,  several  conclusions 
of  a  qualitative  nature  that  can  be  stated  with  respect  to 
productivity  associated  with  a  declining  orderbook. 

1.  When  employees  know  that  they  are  working  themselves 
out  of  a  job,  there  is  certainly  no  incentive  to 
either  maintain  their  level  of  productivity  or  to 
improve  it.  This  is  a  known  fact  to  shipyard  manage¬ 
ment  but  difficult  to  document. 
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2.  Ongoing  shipyard  overhead  costs  will  be  distributed 

over  a  smaller  product  base.  This  adversely  affects 
the  shipyard's  capability  to  economically  compete  for 
new  work.  ' 

3.  Lower  paid  employees  with  lower  seniority  will  be 
separated  first  during  a  manpower  reduction.  This  will 
result  in  an  increase  in  the  average  wage  rate  for  those 
employees  retained.  Such  an  increase  in  average  wage 
will  have  an  adverse  impact  on  the  shipyard's  capability 
to  economically  compete  for  new  work. 

4.  A  declining  orderbook  will  result  in  under  utilization 
of  existing  capacity  adversely  impacting  on  efficient 
production  control. 

5.  There  will  be  no  incentive  for  capital  investment  to 
further  improve  productivity  when  a  shipyard's  order- 
book  is  declining.  In  shipyards  where  recent  large 
scale  private  capital  investments  were  triggered  by 

the  Merchant  Marine  Act  of  1970,  increased  productivity 
benefits  have  directly  accrued  to  the  government;  that 
is  to  subsequent  Navy  shipbuilding  programs  in  those 
yards . 

6.  Even  when  a  declining  orderbook  is  reversed  and  new 
contracts  are  received,  there  will  be  reduced  pro¬ 
ductivity  associated  with  new  hires  until  they  are 
fully  trained  and  capable  at  their  jobs.  The 
employee  recall  success  rate  in  shipyards  is  about 
50  percent.  Reduced  productivity  must  therefore  be 
recognized  for  the  other  50  percent  until  they  are 
fully  capable  employees. 

7.  There  will  be  direct  costs  of  training  new  employees 
during  a  recovery  period."  (Karlson,  1979) 

The  above  statements  suggests  the  impact  of  a  change  in  the  demand  for 
ships  on  productivity  and  cost  is  likely  to  be  in  a  number  of  time  periods 
as  opposed  to  being  concentrated  in  a  single  period.  Time  is  required 
to  recruit  and  train  the  new  workers  needed  to  meet  an  increased  demand. 

The  tendency  of  workers  to  slow  down  during  times  of  low  demand  is  likely 
to  persist  for  some  time  after  demand  increases.  These  statements  also 
suggests  that  both  delayed  and  anticipatory  reactions  could  occur.  Memories 
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of  past. layoffs  are  likely  to  make  potential  workers  reluctant  to 
accept  employment  with  a  shipyard.  This  can  cause  the  workforce 
expansion  plans  to  be  delayed.  Memories  of  past  layoffs  can  motivate 
shipyard  workers  to  voluntarily  quit  to  accept  more  secure  jobs  if  they 
suspect  that  the  level  of  demand  is  likely  to  fall  in  the  future.  The 
tendency  for  the  impact  to  be  felt  in  more  than  one  time  period  combined 
with  the  possibilities  of  both  delayed  and  anticipatory  reactions  strongly 
suggests  that  it  will  be  necessary  to  use  either  differential  equations  or 
difference  equations  to  analyze  how  changes  over  time  in  the  demand  for 
ships  impacts  productivity  and  cost. 


i 


29 


REFERENCES 


1.  Department  of  Defense.  Annual  Report  on  the  Status  of  Shipbuilding 

and  Ship  Repair  Industry  of  the  United  States  -  1977 .  Coor¬ 
dinator  of  shipbuilding,  Conversion  and  Repair,  Office  of 
Maritime  Affairs,  Naval  Sea  Systems  Command,  Washington,  D.  C., 
1977. 

2.  Karlson,  E.,  personal  communication,  1979. 

3.  Mark  Battle  Associates,  Inc.,  "Shipbuilding  Manpower  Study:  Executive 

Summary,  March,  1974. 

4.  Muench,  George  A.  Study  for  the  Improvement  of  Motivation  in  the 

Shipbuilding  Industry,  Phase  I,  Prepared  for  the  Maritime  Admin¬ 
istration  National  Shipbuilding  Research  Program  by  Dr.  George 
A.  Muench  and  Associates,  San  Jose  State  University,  San  Jose, 
California,  June,  1976. 


30 


4 


■ 


A 


APPENDIX  A 

SOURCES  AND  REFERENCES  ON  SHIPBUILDING 


31 


Annual  Report  on  the  Status  of 

the  Shipbuilding  and  Ship  Repair 

Industry  of 

the  United  States: 
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2-1  to  2-5 

Research  and  reports 
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Steel  fabrication  and  erection 
Blocking,  storing,  and  staging 
Machinery  installation 
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A 
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B 

45  to  71 
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C 

72  to  78 

Cost  determination  in  production  and 
construction 

D 

79  to  92 

Management  of  construction 

General 

Cargo  Ship  Economics 

• 

and  Design  by  H.  Benford,  University 

of  Michigan,  July  1962. 
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Pages 

Topics 

A 

61  to  70 

Weights,  estimating 

71  to  85 


Building  costs 
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21.  QNR  Project  -  Field  Visit  and  Interview,  by  Forrest  B.  Creen,  GRA , 

VPI  &  SU,  conducted  at  NOrfolk  Naval  Shipyard,  Portsmouth,  Virginia, 
21  December  1978. 

22.  Economics ,  3rd  ed.,  by  Lloyd  G.  Reynolds,  publisher  Richard  D.  Irwin, 
Inc.,  Homewood,  Illinois,  1969  —  Chapter  20,  "Factor  Markets  and 
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23.  Ideal  Labor  Progress  Curves:  New  Construction  and  Conversions, 

NAVSEA,  provided  by  Harold  Paul,  summer,  1978. 

24.  Shipyard  Detail  Information  (Chart),  NAVSEA,  provided  by  Harold 
Paul,  summer,  1978. 

’  • 

25.  Manhours  per  Ton  by  Ship/Type  (Table),  NAVSEA,  provided  by  Harold 
Paul,  summer,  1978. 

26.  Gantt  Charts  and  PERT  Networks,  Norfolk  Naval  Shipyard. 
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27.  Network  Analysis  and  its  Application  to  Shipbuilding,  by  Chambers 
et  al.,  Paper  presented  at  the  Royal  Institute  of  Naval  ARchitects 
meeting,  London,  England,  March  23,  1965. 

2 8 .  Principal  Shipbuilding  and  Repair  Facilities  of  the  United  States , 
DOD  and  DOC,  Office  of  the  Coordinator  for  Ship  Repair  and  Conver¬ 
sion,  June  1,  1975. 

29.  1977  Report  on  Survey  of  U.S.  Shipbuilding  and  Repair  Facilities, 
DOC,  MARAD ,  Office  of  Ship  Construction,  Division  of  Production, 
December  1977. 

30.  Cost  Report  Summary,  Chesapeake  Corporation  of  Virginia,  West 
Point,  Virginia,  October  20,  1977. 

31.  Labor  and  Material  Cost  Report,  prepared  by  the  Rust  Engineering 
Company,  A  five-part  contract  report  compilation  prepared  for  the 
Chesapeake  Corporation  of  Virginia,  October  18,  1977. 

32 .  Cost/Schedule  Control  Systems  Criteria  (C/SCSC)  Joint  Surveillance 
Guide,  Naval  Material  Command  Pamphlet,  NAVMAT  P5243,  July  1,  1974. 

33.  A  Statistical  Summary  of  Shipbuilding  Under  the  U.S.  Maritime  Com- 


mission 

During  W.W.II, 

by  Gerald  J.  Fisher,  Historical  Reports, 

No.  2, 

1949. 

Section 

Pages 
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Programs  planned  and  accomplished 

B 
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C 
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Shipbuilding  by  individual  yards 
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Speed  of  construction 
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90  to  106 
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F 

107  to  113 

Steel 

G 

114  to  147 
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H 
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34.  Military  Standard  Production  Management,  M1L-STD-1528  (USAF) , 
August  1,  1972. 

35.  Department  of  Defense  Instruction,  Number  7000.2,  June  10,  1977, 
Performance  Measurement  for  Selected  Acquisitions*. 

3 6 .  Military  Standard  Work  Breakdown  Structures  for  Defense  Material 
Items,  MIL-STD-881A,  April  28,  1975. 

37.  Naval  Analysis  Programs:  Description  of  Study  Areas,  FY  77-79, 
December  1975. 
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38.  Interview  with  Norman  Blnmberg,  Nuclear  Test  Engineer,  Norfolk 
Naval  Shipyard,  Portsmouth,  Virginia,  by  F.  B.  Green,  September 
8,  1978. 

39.  Confidential  Interview  with  Mr.  William  Nelson,  Fiance  Officer, 
Major  East-Coast  Construction  Firm,  by  F.  B.  Green,  September  15, 
1978. 

40.  Lecture  Notes  on  "Industrial  Economy"  by  Dr.  Franz  A.  P.  Frisch, 
VPI  &  SU,  1977,  Note  Nos.  163-190. 

41.  Lecture  Notes  on  "Industrial  Economy”  by  Dr.  Franz  A.  P.  Frisch, 
VPI  6  SU,  1977,  Note  Nos.  121-162. 

42.  An  Overview  of  the  Shipbuilding  Labor  Market,  by  Sheldon  E.  Haber 
and  John  C.  Martin,  George  Washington  University,  Institute  for 
Management  Science  and  Engineering,  Serial  T-331,  March  24,  1976. 

43.  The  Conglomerate  Firm,  edited  by  Arnold  C.  Max  and  Zenon  S. 
Zannetos,  Technical  Report  No.  5,  Sloan  School  of  Management,  MIT, 
June  1978. 

Section  Pages  Topic _ 

A  76  to  104  Conglomerates  in  the  shipbuilding 

industry 

B  185  to  209  The  shipbuilding  company  in  a  multi¬ 

divisional  corporate  environment 

44.  Annual  Report  of  the  Maritime  Administration  [not  on  board]. 


45.  Shipyard  Employment  Report,  monthly,  Office  of  Maritime  Manpower, 
Maritime  Administration  [not  on  board]. 

46.  Private  Shipbuilding  and  Conversion  Shipyards  -  Facilities  Data, 
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bility  of  yards 

Size,  capacity,  and  facilities  of 
shipyards 


4 7 .  Development  and  Application  of  a  Computer  Controlled  Ship's  Frame 
Bender  in  the  Automated  Shipyard,  by  Acton  et  a].,  paper  presented 
at  the  Annua]  Meeting  of  the  Society  of  Naval  Architects  and 
Marine  Engineers,  November  16-18,  1978. 


48.  Ways  to  Increase  Shipbuilding  Productivity,  Report  to  the  Congress 
by  the  Comptroller  General  of  the  United  States,  September  23,  1976. 


36 


49.  Brevu  -  Bremer  Vulcan,  Schiffbau  und  Maschinef abrik  (Shipbuilders 
and  Marine  Engin-ers)  ,  brochure  provided  by  Berner  Vulcan,  Weser- 
strabe,  Germany,  1978. 

5 0 .  The  Profitability  of  the  U.S.  Shipbuilding  Industry,  1947-19 76  by 
Edward  M.  Kaitz,  report  prepared  under  Acquisition  RM)  Program  of 
the  Office  of  Naval  Research,  June  20,  1978. 

51.  Productivity,  extracted  from  internal  report  prepared  at  MARAD  in 
December  1977,  provided  by  E.  Karlson. 

52 .  Proceedings  for  Shipbuilding/Industrial  Production  Engineering 
Workshop ,  U.S.  Dept,  of  Commerce,  Maritime  Administration,  in 
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24,  1978. 

Section 

Pages 

Topic 

A 

I 1-1-6 

Production  planning,  scheduling  and 
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Methods  and  standards 
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53 .  Aspects  of  the  National  Shipbuilding  Research  Program  which  Impact 
on  Owners,  Designers,  Regulators,  and  Suppliers,  by  L.  D.  Chirillo, 

R  &  D  Program  Manager,  Todd  Shipyards  Corporation,  Seattle  Division, 
1  October  1977. 

54.  Outfitting  and  Production  Aids  by  Lou  Chirillo  and  others,  March  7, 
1979. 
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"The  Future  of  the  Navy,"  U.S.  Naval  Institute  Proceedings,  June 
1977,  pp.  83-84;  "Attack  Submarine  Development-Recent  Trends  and 
Projected  Needs,"  U.S.  Naval  Institute  Proceedings,  November  1978. 

56.  Questionnaire:  National  Shipbuilding  Research  Program  Product 
Oriented  Work  Breakdown  Structure  by  L.  D.  Chirillo,  R  &  D  Program 
Manager,  Todd  Pacific  Shipyards  Corporation,  Seattle,  Wash., 

December  1,  1978. 

57.  "Industrial  Dynamics:  A  Major  Breakthrough  for  Decision  Makers," 
Harvard  Business  Review,  July-August,  1958. 

58.  Study  no.  T-383,  The  United  States  Shipbuilding  Labor  Market  1960- 
1970,  extract:  author  unknown. 
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59.  The  National  Shipbuilding  Research  Program,  "Study  for  the  Improve¬ 
ment  of  Motivation  in  the  Shipbuilding  Industry,"  by  the  U.S.  Dept, 
of  Commerce  Maritime  Administration  in  cooperation  with  Newport 
News  Shipbuilding  and  Dry  Dock  Company,  June  1976. 

60.  Study  no.  163,  Shipbuilding  Manpower  Study  Executive  Summary,  by 
Mark  Battle  Associates,  Inc.,  Washington,  D.C.,  March  1974. 

61.  Study  no.  163,  Shipbuilding  Manpower  Study,  by  Mark  Battle  Associates, 
Inc.,  March  1974. 

62.  Study  no.  163,  Shipbuilding  Manpower  Study  Appendices,  by  Mark 
Battle  Associates,  Inc.,  March  1974. 

63.  Rational  Approaches  to  Raising  Productivity,  by  Mitchell  Fein, 

PE  Consultant,  WM  &  ME  Monograph  Series  No.  5,  1974. 

64.  Motivation  for  Work,  by  Mitchell  Fein,  WM  &  ME  Monograph  Series  No. 

4,  1974. 

65.  Improved  Planning  and  Production  Control,  Prepared  by  Bath  Iron 
Works  Corporation  under  Marad  Task  0-2  of  the  Ship  Producibi lity 
Research  Program,  August  1977. 

66.  A  Manual  on  Planning  and  Production  Control  for  Shipyard  Use, 

Prepared  by  Corporate-Tech  Planning,  Inc.  for  Bath  Iron  Works 
Corporation  under  MarAd  Task  0-2  of  the  Ship  Producibi! ity  Research 
Program,  September  1978. 

67.  Wage  Incentive  Plans,  by  Mitchell  Fein,  Publication  No.  2,  American 
Institute  of  Industrial  Engineers,  May  5,  1970. 

68.  Potential  Ship  Structure  Misalignment  Tolerance  Guidelines,  Final 
Technical  Report  (Draft),  Part  No.  1  of  Standard  Structural 
Arrangements  Task  S-ll  of  Ship  Producibility  Program,  General 
Dynamics,  Quincy  Shipbuilding  Division. 

69.  Gert  Modeling  and  Simulation:  Fundamentals  and  Applications,  by 
L.  Moore  and  E.  Clayton,  New  York:  Petrocelli-Char ter ,  1976. 

70.  Modeling  and  Analysis  Using  A-GERT  Networks,  by  A.A.B.  Pritsker, 
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SHIPYARD  INTERVIEW  QUESTIONNAIRE 


Introduction: 

My  name  is  _  and  I  am  calling  from  the  Industrial 

Engineering  and  Operations  Research  Department  of  Virginia  Polytechnic 
Institute.  We  are  currently  engaged  in  a  research  project  sponsored  by  the 
Office  of  Naval  Research  entitled  "Analysis  of  the  Cost  of  Variable  Workloads 
on  Shipbuilding".  The  project  is  part  of  a  basic  research  program  in  the 
area  of  ship  Acquisition. 

I  have  been  asked  to  contact  people  at  various  shipyards  in  order  to 
gather  information  on  some  of  the  areas  in  which  our  research  will  be  directed. 
I  have  been  assured  that  all  data  will  be  treated  confidentiality,  and  no  names 
or  titles  will  be  used  without  permission.  If  you  wish  to  verify  authorization 
for  this  project,  I  can  give  you  the  ONR  Project  Number. 

Project  Number  is  230-11-022-106-352965-1. 


Questions; 

In  your  position,  what  phase  of  ship  construction  are  you  most  concerned  with? 
(  )  Hull  fabrication 

(  )  Propulsion  machinery  installation 
(  )  Auxilliary  machinery  installation 
(  )  Outfitting 

(  )  Other  _ 

What  trades  are  directly  involved  in  the  type  of  work  you  supervise? 


At  what  stage  of  ship  construction  do  you  experience  the  most  activity? 

(  )  Pre-construction,  design 

(  )  Early  construction  (up  to  33%  complete) 

(  )  Mid-construction  (33%  to  67%  complete) 

(  )  Final  stages  (67%  to  100%  complete) 

(  )  Post-construction,  retro-fitting  and  overhaul. 

Which  area  of  the  shipyard  (or  shipboard  space  or  compartment)  involves  the  most 
activity  that  you  experience?  (Where  in  the  shipyard  is  the  most  activity 
experienced?)  * 


Which  activity  or  group  of  activities  in  that  area  or  space  involve  the  highest 
concentration  of  personnel? 


What  are  the  trades  involved  in  this  activity  or  group  of  activities? 


We  recognize  that  most  shipyards  do  a  very  good  job  of  meeting  schedules  and 
filling  orders;  but  variations  in  productivity  do  occur,  and  we  are  trying  to 
focus  our  attention  on  the  reasons  for  these  variations.  If  you  could  jot  these 
things  down,  I  am  going  to  give  you  a  list  of  items  which  may  cause  variation  in 
job  completion  time  and  ask  you  to  rate  the  top  three  items  in  order  of  importance 

_  Unavailability  of  labor  at  the  time  needed  (low  labor  progress) 

_  Insufficient  skilled  personnel  (low  manning  level) 

_  Improper  order  or  sequence  of  assigned  labor 

_  Interference  by  different  crafts  (incompatibility) 

_  Too  many  personnel  working  at  one  time  (work  density) 

_  Accelerated  schedules  (work  speedup) 

Other  causes  _  _ * _ 
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(Ask  the  following  questions  if  interference  or  work  density  is  identified  as  a 
cause;  otherwise,  go  to  the  next-to-the-last  question.) 

Where  does  interference  between  workers  or  work  crews  occur  in  the  shipbuilding 
process? 


When  does  this  happen? 


l!ow  do  you  handle  it  if  it  is  a  problem? 


Which  types  of  crews  or  trades  are  most  compatible? 


Which  types  of  crews  or  trades  are  least  compatible? 


Is  this  considered  a  serious  problem? 
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What  type  of  organizational  structure  does  your  shipyard  have? 


(  )  Heirarchical 

(  )  Project  manager 

(  )  Area  coordinator 

(  )  Matrix 

(  )  Other  _ 


Who  in  the  shipyard  should  we  talk  to  if  more  details  are  needed?  (Get  name 
and  phone  number.) 


Who  should  we  contact  for  authorization  to  visit  the  yard? 


CONCLUSION:  (Thank  the  respondent  for  his  participation) 
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APPENDIX  D 

SHIPYARDS  INTERVIEWED 


1 

1^ 

Shipyard 

Size 

Merchant 

Naval 

1 

Bay  Shipbuilding 

Major 

X 

1 

Sun  Shipbuilding 

Ma  j  or 

X 

Bethlehem  Steel 

Major 

X 

1 

Ingall's  Shipbuilding 

Major 

X 

1 

Todd  Shipbuilding 

Major 

X 

X 

1 

Jeff  Boat 

Int 'm 

X 

1 

National  Steel 

Major 

X 

X 

1 

American  Shipbuilding 

Ma  j  or 

X 

1 

Lockheed  Shipbuilding 

Major 

X 

1 

Dravo  Corporation 

Int'm 

X 

1 

Bath  Iron  Works 

Major 

X 

V 

s\ 

1 

Equitable  Shipyards 

Int 'm 

X 

w 

Levingston  Shipbuilding 

Int  'm 

X 

1 

Alabama  Shipbuilding 

Int  'm 

X 

1 

Zigler  Shipyards 

Small 

X 

1 

1  ' 

Avondale 

Major 

X 

X 

■ 

McDermott  Shipyard 

Small 

X 

Nashville  Bridge  Company 

Small 

X 

1 

Peterson  Builders 

Small 

X 

1 

Norfolk  Shipbuilding 

Int  'm 

X 

1 

Newport  News  Shipbuilding 

Major 

X 

• 

X 

1 

*• 

V 

l 
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PREFACE 


The  Industrial  Engineering  and  Operations  Research  Department's 
activities  were  divided  into  three  phases.  The  first  phase,  which  is 
reported  in  "A  Diagnosis  of  the  Workload  Variation  Problem  in  Ship¬ 
building",  was  designed  to  identify  specific  problems  for  further 
analysis.  This  analysis  revealed  the  need  for  the  following: 

(1)  a  system  for  analyzing  how  variations  in  shipyard  workloads 
affect  shipbuilding  cost; 

(2)  a  shipyard  planning  system  which  could  be  utilized  by  ship¬ 
yards  to  minimize  the  cost  of  adjusting  to  workload  variations. 

Phase  2,  which  is  the  subject  of  this  report,  focused  on  developing 
the  framework  of  a  system  for  developing  transfer  function  models  which 
can  be  used  to  develop  the  difference  equations  which  describe  how 
variations  in  shipyard  workloads  affect  shipbuilding  cost.  Phase  3, 
which  is  described  in  "A  Shipyard  Planning  System",  focused  on  develop¬ 
ing  the  framework  of  a  shipyard  planning  system  which  shipyards  could 
use  to  minimize  the  cost  of  adjusting  to  workload  variations. 
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I. 


INTRODUCTION 


The  purpose  of  this  paper  is  to  develop  a  framework  for  analyzing 
how  variation  in  shipyard  workloads  impact  shipbuilding  cost.  Such 
analysis  could  be  used  by  those  responsible  for  establishing  Navy  ship¬ 
building  budgets  to  analyze  how  alternative  funding  programs  influence 
both  procurement  cost  and  mobilization  potential. 

In  order  to  analyze  how  variations  in  workload  affect  shipbuilding 
cost,  it  was  necessary  to  develop  a  means  for  quantifying  the  terms: 

(1)  workload,  and  (2)  shipbuilding  cost.  It  was  also  necessary  to 
establish  a  model  building  process  for  developing  the  difference 
equations  which  describe  how  workload  variations  dynamically  affect 
shipbuilding  cost. 


II.  OPERATIONAL  DEFINITIONS 

The  first  step  towards  quantifying  workload,  and  (2)  shipbuilding 
cost  was  to  develop  operational  definitions  for  these  terms.  The 
process  for  developing  an  operational  definition  of  a  term  was  to  first 
develop  a  precise  phenomenological  definition  of  the  term.  This 
definition  will  be  used  to  construct  an  idealized  measurement  procedure 
where  the  only  constraint  would  be  whether  it  was  physically  feasible — 
for  example,  the  information  systems  cost  will  not  be  considered  at 
this  stage.  This  idealized  measurement  procedure  will  then  be  used  as 
a  criterion  for  evaluating  the  relative  merits  of  alternative  measures 
which  can  be  derived  from  currently  available  sources  of  data. 

The  term  "workload"  is  defined  by  Urdang  and  Flexner  (1969)  as  "the 
amount  of  work  that  a  machine,  employee,  or  group  of  employees  can 
be  or  is  expected  to  perform."  In  practice  the  amount  of  work  which  a 
man  and/or  machine  system  is  expected  to  perform  under  a  given  set  of 
conditions  is  specified  by  a  production  (or  construction)  schedule. 

This  suggests  that  workload  be  redefined  as  the  amount  of  work  which  a 
man  and/or  machine  system  is  scheduled  to  perform  during  a  specified 
time  period  under  specified  conditions. 

In  theory  methods  engineering  can  be  used  to  determine  the  capacity 
of  a  production  facility  operating  in  a  construction  type  environment 
similar  to  those  encountered  in  shipbuilding  (Parker  and  Oglesbuy,  1972). 
The  major  problem  is  that  methods  engineering  has  not  been  used  to 
determine  shipyard  production  capacities  in  the  past  (Bath  Iron  Works 
Corporation,  1977).  Therefore,  such  information  could  not  be  obtained 
from  historical  data  even  if  cost  were  of  no  concern.  Another  problem 
is  that  space  and  facilities  are  limited  in  most  shipyards.  This 
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creates  problems  in  which  workers  in  the  same  craft  and  in  different 


crafts  can  interfere  with  the  productivity  of  one  another.  Presently 
available  methods  engineering  techniques  are  not  capable  of  predicting 
a  priori  when  and  the  extent  to  which  such  interference  will  result. 
However,  time-lapse  motion  pictures  have  been  used  to  obtain  such 
information  on  an  empirical  basis  for  specific  situations  in  civil 
construction  (Parker  and  Oglesby,  1972). 

In  most  organizations  scheduling  is  performed  by  a  staff  group 
such  as  production  control,  production  planning,  or  production 
scheduling,  who  produce  formal  documents  known  as  production  schedules. 
Once  the  schedule  work  has  been  accomplished,  the  production  schedule 
is  of  little  or  no  value.  When  the  value  of  the  information  contained 
in  a  document  is  less  than  the  cost  of  storing  it  the  document  should 
be  destroyed.  Therefore,  the  amount  of  historical  data  on  manufacturing 
facility  workloads  will  be  limited. 

The  term  "shipbuilding  cost"  was  decomposed  into  its  constituent 
parts.  Shipbuilding  is  an  umbrella  term  which  refers  to  the  activities 
involved  in  constructing  a  naval  vessel.  However,  the  set  of  activities 
necessary  for  constructing  one  ship  may  differ  greatly  from  that  of 
another  ship.  Execution  times  for  ship  constructions  are  likely  to 
decrease  as  the  shipyard  becomes  more  experienced  due  to  the  "learning" 
effect.  This  is  especially  pronounced  with  the  construction  of 
identical  (or  near  identical)  vessels.  This  often  leads.to  method 
improvements  in  the  construction  of  subsequent  ships.  Difficulties 
arise  in  predicting  costs  when  the  same  work  on  two  identical  ships  can 
be  performed  by  different  production  facilities  which  utilize  vastly 


different  work  methods.  Also,  weather  conditions  can  also  necessitate 
utilizing  different  production  facilities  and  work  methods  for  accomplish¬ 
ing  the  same  work  package  on  a  particular  ship. 

The  term  "cost"  is  defined  by  Urdang  and  Flexner  (1969)  as  "the 
price  paid  to  acquire,  produce,  accomplish,  or  maintain  anything". 

Note  that  this  definition  gives  the  reader  the  choice  between  the 
following  words:  (1)  acquire;  (2)  produce;  (3)  accomplish;  or  (4) 
maintain.  In  choosing  the  most  appropriate  word  it  will  be  helpful  to 
consider  the  wide  variety  of  activities  which  must  be  performed  in 
building  a  ship.  Raw  materials  and  purchased  parts  components  and 
systems  are  acquired  from  outside  resources.  The  raw  materials  are 
first  converted  into  manufactured  parts.  If  a  large  number  of  identical 
or  highly  similar  parts  are  produced  successively  on  a  periodic  basis, 
then  the  conversion  process  will  be  referred  to  as  a  continuous  produc¬ 
tion  process.  If  only  one  or  a  limited  number  of  copies  of  a  part  are 
made  at  one  time  the  conversion  process  will  be  referred  to  as  a 
discrete  production  process.  A  ship  is  built  by  assembling  manufactured 
purchases  (or  government  supplied)  parts,  components,  and  systems 
according  to  a  specified  plan.  The  process  of  assembling  is  referred  to 
as  the  construction  process.  In  order  to  perform  these  activities 
efficiently  it  is  necessary  for  the  shipyard  to  maintain  its  facilities. 

The  above  definition  of  cost  is  too  limited  to  account  for  the  wide 
variety  of  activities  which  must  be  performed  in  building  a  ship.  In 
order  for  the  definition  to  encompass  those  activities  the  definition 
should  be  modified  so  that  the  reader  is  given  a  choice  between  the 
following  words:  (1)  acquire;  (2)  convert;  (3)  construct;  and  (4) 
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maintain.  Therefore,  the  above  definition  of  cost  should  be  modified 
so  as  to  account  for  each  of  the  above  activities. 

The  price  associated  with  accomplishing  each  of  the  above  activities 
will  typically  be  measured  in  monetary  units.  This  poS.es  a  problem  since 
the  value  of  all  monetary  units  have  been  erroded  by  inflation  during 
recent  years.  Fortunately,  various  indices  are  available  which  can  be 
used  to  adjust  for  inflation.  Another  problem  arises  due  to  the  fact 
that  payments  made  for  accomplishing  the  above  activities  are  typically 
spread  out  over  a  three  to  five  year  time  period.  Each  payment  represents 
the  loss  of  the  opportunity  to  utilize  the  funds  in  alternative  ways. 

The  price  paid  to  accomplish  a  given  activity  can  be  classified 
according  to  whether  the  outlay  was  made  during  the  present  period  or 
in  some  previous  period.  For  example,  wages  will  be  paid  in  the  current 
period,  while  the  cost  of  acquiring  a  capital  asset  may  have  been 
incurred  during  some  previous  period. 

The  cost  of  acquiring  a  capital  asset  equipment  consists  of  the 
purchase  cost,  transportation  cost  and  installation  cost.  From  a  legal 
standpoint  the  costs  are  incurred  at  the  time  title  changes  hands.  The 
mere  fact  that  a  machine  is  purchased  on  credit  with  payments  being 
spread  out  over  time  does  not  alter  the  fact  that  the  cost  of  acquiring 
a  capital  asset  is  incurred  at  the  time  the  title  changes  hands.  How¬ 
ever,  such  costs,  typically  incurred  in  the  past,  complicates  the  task 
of  determining  the  cost  associated  with  accomplishing  the  activities 
performed  during  the  present  period. 

The  problem  revolves  around  determining  the  value  of  the  capital 
assets  which  were  consumed  as  a  result  of  performing  the  current  periods 


activities.  All  solutions  are  arbitrary  to  a  certain  extent.  The 
"generally  accepted"  accounting  practice  is  to  record  the  value  of  the 
asset  at  its  original  acquisition  cost.  The  value  of  the  asset  is  then 
spread  over  its  useful  life  according  to  some  depreciation  method. 
Depreciation  methods  can  be  classified  as  either  straight  line  or 
accelerated  methods.  The  straight  line  method  changes  an  equal  amount 
of  the  assets  original  value  to  each  year  of  its  useful  life.  In 
relation  to  this  straight  line  method  the  accelerated  methods  charge 
higher  amounts  during  the  early  years  and  lower  amounts  during  the 
latter  years  of  the  assets  life. 

The  "generally  accepted"  accounting  practice  of  recording  the  value 
of  an  asset  at  its  original  acquisition  cost  is  undesirable  from  the 
standpoint  that  the  total  amount  of  the  depreciation  charges  will  not 
be  sufficient  to  cover  the  cost  of  acquiring  an  identical  machine  at  • 
the  time  the  asset  has  been  depreciated  to  zero.  Therefore,  it  appears 
that  current  accounting  practices  are  underestimating  the  value  of  the 
capital  assets  which  are  consumed  in  the  manufacturing  process.  An 
alternative  method  of  calculating  depreciation  charges  has  been  proposed 
but  at  present  has  not  found  wide  acceptance.  This  method  involves 
basing  depreciation  charges  on  the  replacement  value  of  the  asset  rather 
than  on  its  original  cost. 

The  analysis  of  the  terms  shipbuilding  and  cost  provided  a  basis 
for  synthesizing  the  following  phenomenological  definition  of  shipbuild¬ 
ing  cost. 

Shipbuilding  cost  is  the  amount  paid  for  performing  the  following 
activity  involved  in  building  a  ship: 
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(1)  acquiring  raw  materials,  and  purchased  parts,  components, 
and  systems  from  outside  sources; 

(2)  converting  raw  materials  into  manufactured  components; 

(3)  assembling  and  joining  the  various  manufactured  parts,  and 
purchased  and/or  government  supplied  parts,  components,  and 
systems  according  to  a  specified  plan; 

(4)  maintaining  manufacturing  facilities  where  the  price  is 
adjusted  for  the  effects  of  inflation,  and  the  time  value 
of  money,  and  where  depreciation  charges  which  are  included 
in  the  price  are  based  on  the  cost  of  replacing  the  capital 
assets . 

In  theory  shipbuilding  cost  as  defined  above  could  be  calculated  from 
presently  available  data.  However,  in  practice,  it  would  be  difficult 
to  determine  the  replacement  cost  for  assets  during  prior  time  periods. 
Thus,  the  amount  of  historical  data  on  shipbuilding  cost  will  be  so 
limited . 

One  strategy  for  coping  with  the  limited  data  problem  is  to  iden¬ 
tify  historically  available  variables  which  logically  would  be  expected 
to  behave  in  a  similar  way  as  manufacturing  facility  workload  and  ship¬ 
building  costs.  Such  variables  will  be  referred  to  as  "proxy  variables". 
A  proxy  variable  for  manufacturing  facility  workload  is  suggested  by  the 
fact  that  this  quantity  was  defined  in  terms  of  the  amount  of  work  that 
a  man  and/or  machine  system  is  scheduled  to  perform  during  a  specified 
time  period.  In  a  labor  intensive  industry  such  as  shipbuilding  it  is 
reasonable  to  assume  that  the  amount  of  work  that  management  expects  to 
be  accomplished  will  be  roughly  proportional  to  the  number  of  workers 
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assigned  to  the  manufacturing  facility.  Monthly  data  on  the  number  of 
production  workers  in  the  shipbuilding  industry  can  be  obtained  from  the 
Bureau  of  Labor  Statistics.  This  data  series  was  started  in  1947. 

Similar  data  is  also  available  at  the  shipyard  level.  However,  special 
permission  will  be  needed  to  gain  access  to  this  data-. 

A  proxy  variable  for  shipbuilding  cost  is  suggested  by  the  fact  that 
approximately  77  percent  of  the  value  added  by  the  shipbuilding  industry 
is  due  to  labor  (Martin,  1977) .  Furthermore,  this  quantity  has  remained 
remarkably  stable  since  1954.  Since  labor  cost  dominates  value  added, 
this  suggests  that  constant  dollar  labor  cost  per  ton  might  be  used  as  a 
proxy  variable  for  shipbuilding  cost.  Shipbuilding  labor  cost  in 
current  dollars  are  reported  on  monthly  basis  by  the  Bureau  of  Labor 
Statistics.  This  data  can  be  converted  to  constant  dollar  labor  cost  by 
a  price  index  (such  as  the  Consumer  Price  Index)  as  a  deflator.  The  • 
labor  cost  is  available  on  a  monthly  basis.  The  labor  cost  data  series 
was  started  in  1948.  However,  tonnage  data  for  shipbuilding  are  not 
available  on  a  monthly  basis. 

The  failure  to  find  adequate  proxy  variables  for  workload  and  ship¬ 
building  cost  which  could  be  derived  from  historical  data  signaled  the 
need  for  a  fresh  approach.  This  led  to  a  qualitative  analysis  of  the 
ways  in  which  workload  variations  could  influence  shipbuilding  cost. 

The  results  of  this  analysis  are  described  in  the  following  section. 
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III.  A  QUALITATIVE  ANALYSIS  OF  THE  IMPACT  OF  WORKLOAD  VARIATIONS  ON 

SHIPBUILDING  COST 

In  the  preceeding  section  operational  definitions  for  workload 
and  shipbuilding  cost  were  developed.  The  definitions  were  then  used 
as  criteria  for  evaluating  the  relative  merits  of  alternative  proxy 
measures  which  could  be  derived  from  historical  data.  The  failure  to 
find  adequate  proxy  variables  which  could  be  derived  from  historical 
data  indicated  the  need  for  a  fresh  approach.  This  section  describes 
the  approach  which  was  adopted.  This  approach  was  .o  analyze  in  a 
qualitative  fashion  the  ways  in  which  workload  variations  could  influence 
shipbuilding  cost. 

In  the  previous  section  workload  was  defined  as  "the  amount  of 
work  which  a  man  and/or  machine  system  is  scheduled  to  perform  during 
a  specified  time  period  under  specified  conditions".  This  suggests 
that  workload  variation  is  a  change  in  the  scheduled  amount  of  work. 

The  fact  that  production/construction  schedules  are  established  by 
shipyard  management  prompted  a  search  for  the  factors  which  could 
cause  production/construction  schedules  to  be  altered.  This  analysis 
revealed  that  changes  in  the  demand  for  ships  is  the  prime  factor  which 
causes  schedule  alterations.  Shortages  of  labor  and  materials  were 
also  found  to  be  important  factors. 

The  above  analysis  indicated  that  the  demand  for  ships  should  be 
regarded  as  a  causal  variable  and  that  workload  should  be  treated  as 
an  intervening  variable.  The  next  step  was  to  analyze  how  changes  in 
the  demand  for  ships  impacts  shipbuilding  cost.  This  analysis  that  a 
change  in  demand  could  cause  reaction  by  shipyard  management  and  ship¬ 
yard  workers  which  could  affect  shipbuilding  cost.  It  was  found  that 
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shipyard  managers  could  respond  to  a  change  in  demand  by  utilizing  one 
or  more  of  the  following: 

(1)  layoff  surplus  workers; 

(2)  hire  new  workers; 

(3)  rehire  workers  previously  laid  off; 

(4)  work  overtime; 

(5)  work  less  than  40  hours  per  week; 

(6)  subcontract  work; 

(7)  use  surplus  production'  workers  to  do  maintenance  work  on 
shipyard  facilities. 

Certain  costs  will  be  associated  with  utilizing  each  of  the  strategies. 
Increasing  the  size  of  the  work  force  can  involve  the  following 
activities:  (1)  recruiting;  (2)  interviewing;  (3)  testing;  (4)  perform¬ 

ing  medical  exams;  (5)  placing;  and  (6)  training.  The  following  often 
result  when  the  size  of  the  workforce  is  decreased:  (1)  exit  interviews; 
(2)  separation  payments;  (3)  increased  unemployment  insurance  premiums; 
and  (4)  a  bad  image  in  the  local  labor  market.  Hancock  (1971,  p.  7-113) 
stated  that  the  cost  for  "the  leaving  of  one  person  and  the  hiring  of  a 
new  one  are  from  5600  to  $2,000".  These  costs  do  not  account  for  the 
hyper  inflation  which  has  transpired  since  1971.  If  an  adjustment  for 
this  is  made  then  the  current  dollar  equivalent  of  these  costs  would  be 
approximately  $1,050  to  $3,500.  Working  overtime  will  require  the  pay¬ 
ment  of  overtime  premiums.  Furthermore,  efficiency  is  lfkely  to  drop. 
Working  less  than  40  hours  per  week  will  tend  to  increase  turnover  cost. 
Subcontracting  work  can  involve  additional  cost  since  the  subcontractor 
must  earn  a  profit.  In  addition  the  subcontractor  is  likely  to  be  less 
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efficient  than  the  shipyard.  Utilizing  skilled  construction  workers  to  do 
maintenance  work  on  shipyard  facilities  will  increase  cost  if  such  work  could 
be  done  equally  well  by  less  skilled  workers.  Obviously,  a  shipyard  would 
like  to  identify  the  combination  of  the  above  options  which  will  minimize 
total  relevant  cost.  A  system  for  doing  this  is  described  in  Terry  (1979). 

It  was  also  found  that  a  change  in  demand  could  motivate  present 
shipyard  workers  to  act  so  as  to  protect  their  interest.  For  example, 
if  workers  suspect  that  they  are  working  themselves  out  of  a  job  they 
will  tend  to  slow  down  to  protect  their  jobs.  In  addition,  some  of  the 
workers  will  probably  quit  to  accept  more  secure  jobs.  This  could  create 
disruptions  which  increase  cost.  In  addition  memories  of  previous  layoffs 
will  make  it  more  difficult  for  a  shipyard  to  recruit  new  workers  once 
demand  increases.  Furthermore,  the  output  rate  of  the  shipyard  will  not 
immediately  increase  once  new  employees  are  hired.  Instead,  it  will  tend 
to  gradually  increase  as  the  new  employees  learn  their  jobs  and  are 
integrated  into  the  workforce. 

The  above  analysis  suggested  the  systems  model  depicted  in  Figure  1 
on  the  following  page.  This  model  contains  three  black  boxes  which 
represent  shipyard  management,  present  shipyard  workers  and  potential 
shipyard  workers.  The  demand  for  ships  serves  as  an  input  to  these  black 
boxes.  Each  of  these  black  boxes  transforms  the  input  into  a  variety  of 
outputs.  The  output  for  shipyard  management  can  be  any  combination  of 
the  following:  (1)  layoff  unneeded  workers;  (2)  hire  new  workers;  (3) 
rehire  workers  previously  laidoff;  (4)  work  overtime;  (5)  work  short  week; 

(6)  subcontract;  and  (7)  use  surplus  workers  to  do  maintenance  work  on 
shipyard  facilities.  The  output  for  present  shipyard  workers  consists  of 
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Rehires 


Figure  1. 


Systems  model  of  the  way  changes  in  the 
demand  for  ships  impacts  shipbuilding  cost. 


decision  to:  (1)  slow  down  to  stretch  out  work;  and  (2)  quit  to  accept 
more  stable  employment.  The  output  for  future  shipyard  workers  consists 
of  the  formation  of  opinions  regarding  the  relative  desirability  of 
working  for  a  shipyard.  These  opinions,  if  negative,  could  make  it 
more  difficult  to  recruit  additional  workers. 

The  above  analysis  suggested  that  order  backlog  could  be  used  as  a 
proxy  variable  for  the  demand  for  shipyard  services.  Order  backlog  data 
for  the  shipbuilding  and  repair  industry  is  maintained  on  a  monthly  basis 
by  the  Census  Bureau  in  their  Current  Industrial  Reports  (Census  Bureau) . 

The  above  analysis  also  prompted  a  search  for  data  from  which 
quantitative  measures  of  the  outputs  of  the  black  boxes  in  Figure  1  could 
be  obtained.  This  search  revealed  that  the  following  variables  are 
recorded  on  a  monthly  basis  by  the  Bureau  of  Labor  Statistics  in  their 
Employment  and  Earnings  Series  (Bureau  of  Labor  Statistics) : 

(1)  production-worker  average  weekly  hours; 

(2)  production-worker  average  weekly  overtime  hours; 

(3)  accession  (permanent  and  temporary  additions  including  both 
new  and  rehired  employees)  per  100  employees; 

(4)  new  hires  (permanent  and  temporary  additions  of  new  employees) 
per  100  employees; 

(5)  separations  (terminations  of  employment  initiated  by  either 
employer  or  employee)  per  100  employees; 

(6)  quits  (terminations  of  employment  initiated  by  employees)  per 
100  employees; 

(7)  layoffs  (suspensions  without  pay  lasting  or  expected  to  last 
more  than  7  consecutive  calendar  days,  initiated  by  the  employer 


without  prejudice  to  the  worker)  per  100  employees. 

It  also  provides  a  basis  for  calculating  rehires: 

rehires  =  accession  -  new  hires. 

Furthermore,  the  extent  to  which  short  work  weeks  are  being  utilized 
can  be  measured  in  terms  of  the  amount  by  which  the  production-workers 
average  weekly  hours  falls  short  of  40.  However,  no  means  was  found  for 
quantifying:  (1)  the  extent  to  which  present  shipyard  workers  slow  down 

in  order  to  protect  their  job;  (2)  the  extent  to  which  potential  ship¬ 
yard  workers  are  easy  or  difficult  to  recruit;  and  (3)  the  extent  to 
which  shipyards  use  more  or  less  subcontracting  than  normal. 


IV.  MODEL  BUILDING  PROCESS 


The  purpose  of  this  section  is  to  describe  a  process  which  can  be 
used  to  develop  models  for  analyzing  how  both  shipyard  management  and 
actual  and  potential  shipyard  workers  react  to  changes  in  shipyard 
order  backlogs.  Management  can  react  by  laying  off  unneeded  workers, 
hiring  new  workers  or  recalling  old  workers,  working  overtime,  or 
reducing  the  length  of  the  workweek.  Actual  shipyard  workers  can  react 
by  quitting  or  slowing  dowm  in  order  to  stretch  their  work  out.  Potential 
shipyard  workers  can  either  respond  or  not  respond  to  shipyard  recruitment 
efforts.  This  suggests  the  need  for  models  to  predict  the  behavior  of 

(1)  shipyard  workers  and  (2)  shipyard  management. 

The  strategy  will  be  to  develop  separate  models  to  explain  how  each 
of  the  following  proxy  variables  for  shipbuilding  cost: 

(1)  production  worker-average  weekly  hours; 

(2)  accessions; 

(3)  new  hires; 

(4)  separations; 

(5)  quits; 

(6)  layoffs; 

are  influence  by  order  backlog.  In  developing  these  models  each  of  the 
above  proxy  variables  for  shipbuilding  cost  will  be  regarded  as  an  output 
of  a  black  box.  Each  of  these  black  boxes  will  have  the  same  input: 
order  backlog.  These  black  boxes  represent  the  process  whereby  changes 
in  the  input  variable  (order  backlog)  are  transformed  into  changes  in 
the  output  variable  (one  of  the  proxy  variables  for  the  component  of 
shipbuilding  cost  which  can  be  influenced  by  variations  in  the  demand 
for  shipyard  services. 
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The  task  of  developing  a  mathematical  model  which  adequately 
describes  the  transformation  is  complicated  by  the  following  factors: 

(1)  A  change  in  the  input  will  not  necessarily  impact  the  output 
immediately.  Instead,  the  reaction  can  be  delayed.  This 
delayed  reaction  effect  is  illustrated  by  Figure  2. 

(2)  The  impact  of  a  one  time  change  in  input  on  the  output  does 
not  have  to  be  concentrated  at  a  single  point  in  time.  It  is 
quite  possible  for  a  change  in  the  input  at  a  single  point  in 
time  to  be  spread  over  several  time  periods.  This  partial 
adjustment  process  is  illustrated  by  Figure  3. 

(3)  The  point  at  which  a  change  in  input  creates  the  maximum  impact 
on  the  output  is  not  known.  This  uncertainty  regarding  the 
location  of  the  point  of  maximum  impact  is  illustrated  by 
Figure  4. 

(4)  The  buildup  to  the  point  of  maximum  impact  can  be  either  slow 
and  gradual  or  fast  and  abrupt  or  some  combination  of  these 
two  extremes.  The  varability  in  the  rate  of  buildup  to  the 
point  of  maximum  impact  is  illustrated  by  Figure  5. 

(5)  The  decline  from  the  point  of  maximum  impact  can  likewise  be 
either  slow  and  gradual  or  fast  and  abrupt  or  some  combination 
of  these  two  extremes.  The  variability  in  the  rate  of  decay 
from  the  point  of  maximum  impact  is  illustrated  by  Figure  6. 

In  order  to  cope  with  this  situation  it  will  be  necessary  to  have  a  model 
building  process  which  can  respond  to  any  mixture  of  the  above  types  of 
behavior.  Such  a  model  building  has  been  developed  by  Box  and  Jenkins 
(1970)  and  is  referred  to  as  the  transfer  function  model  building  pro¬ 
cess.  Since  this  process  utilizes  the  Box-Jenkins  univariate  model 
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Figure  2.  Delayed  reaction  of  a  change  in 
output  to  a  change  in  input 
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Figure  3.  Reaction  of  output  to  a  change  in 
input  is  spread  over  several  time 
periods 


Fast  rate  of 
decay 


Figure  5.  Comparison  of  two  different  rates 
of  decay  from  the  point  of  maximum 
impact 
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building  process  this  process  will  be  described  first. 

An  Overview  of  the  Box-Jenkins  Univariate  Method 

A  single  exogenous  shock  can  affect  the  behavior  of  a  system  in  a 
variety  of  ways.  If  a  shock  alters  the  structure  of  the  system,  then 
its  effect  on  the  systems  behavior  will  be  long  term  or  permanent.  A 
shock  which  does  not  alter  the  structure  of  the  system  will  produce  only 
short  term  changes  in  the  behavior.  These  short  term  effects  on  system 
behavior  can  either  end  abruptly  after  a  specified  period  of  time  or 
gradually  decay.  Figures  7,  8,  and  9  respectively  illustrate  the 
permanent  change,  gradual  decay,  and  abrupt  termination  reaction  patterns. 

In  order  to  predict  the  future  behavior  of  a  system  it  is  not 
necessary  to  determine  what  causes,  or  from  where,  an  exogenous  shock 
came.  It  is  more  important  to  determine  how  the  system  reacted  to  an 
exogenous  shock  regardless  of  its  cause  or  source.  However,  if  the 
cause  of  a  shock  is  known  and  if  the  occurrence  of  the  shock  is  preceded 
by  an  early  warning  signal  then  it  might  be  possible  to  predict  when  the 
shocks  would  occur.  If  this  is  the  case  then  this  knowledge  can  be  used 
to  predict  the  future  behavior  of  the  system  of  interest.  If  a  particu¬ 
lar  type  of  shock  always  causes  the  system  to  react  in  a  particular  way, 
then  the  response  of  the  system  to  that  type  of  shock  can  be  described 
by  a  deterministic  model. 

If  for  every  type  of  shock  the  reaction  of  the  system  can  be  de- 
scribed  in  terms  of  a  deterministic  model,  then  all  aspects  of  the 
behavior  of  the  system  can  be  described  completely  by  the  collection  of 
all  such  models.  However,  such  a  collection  of  deterministic  models 
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Figure  7.  Exogoneous  shock  alters  structure  of 
system  which  cause  permanent  change 
in  system  behavior 
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Figure  8.  Exogoneous  shock  affects  system 

behavior  in  a  manner  which  gradually 
decays  with  the  passage  of  time 
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Figure  9.  Exogoneous  shock  affects  system 
behavior  in  a  manner  which 
persists  at  a  constant  level  for 
a  period  of  time  and  then  abruptly 
returns  to  the  original  level 


will  not  necessarily  enable  one  to  predict  future  behavior  of  the  system 
with  complete  certainty.  When  some  or  all  of  the  shocks  are  not  measur¬ 
able  it  will  not  be  possible  to  predict  the  future  behavior  of  the 
system.  Even  when  all  shocks  can  be  measured,  this  information  is  not 
useful  for  predictive  purposes  if  the  time  required  to  make  a  prediction 
exceeds  the  time  required  for  the  shock  to  influence  the  system.  Further¬ 
more,  if  the  cost  of  measuring  the  shocks  plus  the  cost  of  using  this 
information  to  predict  future  behavior  of  the  system  exceeds  the  value 
derived  from  more  accurate  predictions,  then  a  deterministic  model  of 
the  systems  behavior  will  not  be  useful. 

Useful  prediction  models  can  often  be  developed  for  situations  in 
which  either  or  both  of  the  following  conditions  occur:  (a)  it  is 
either  too  expensive  or  impossible  to  measure  the  shocks;  or  (b)  the 
system's  reaction  to  a  shock  can  not  be  described  adequately  by  a  deter¬ 
ministic  model. 

The  strategy  for  developing  predictive  models  in  the  above  situa¬ 
tions  exploits  the  fact  that  the  effects  of  a  shock  on  the  system's 
behavior  tend  to  be  spread  out  over  some  interval  of  time  rather  than 
being  concentrated  at  a  particular  point.  This  suggests  that  the  current 
value  of  a  time  series  which  describes  the  behavior  of  a  system  is  a 
function  of  the  current  and  past  shocks.  Therefore,  it  follows  that  the 
sequence  of  values  which  describes  the  behavior  of  the  system  contains 
information  which  can  be  used  to  establish  a  stochastic  grocess  model 
which  tells  how  the  time  series  which  describes  the  behavior  of  the 
system  was  generated. 
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Box  and  Jenkins  (1970 )  developed  a  systematic  procedure  for  estab¬ 
lishing  such  a  model  for  any  given  time  series.  In  its  simplest  form 
the  Box-Jenkins  method  considers  the  case  in  which  there  is  no  quantita¬ 
tive  information  available  regarding  the  exogenous  forces  which  affect 
the  behavior  of  the  system.  In  any  given  time  period'  it  is  quite 
possible  for  a  number  of  these  exogenous  forces  to  simultaneously  act  on 
the  system.  The  net  effect  of  ail  the  forces  acting  on  the  system  dur¬ 
ing  any  given  time  period  will  be  denoted  by  a^.  The  time  series  at, 
a^_1,  at-2’  •  •  •  assumed  to  be  a  sequence  of  independent,  identically 

O 

distributed  Normal  random  variables  with  mean  zero  and  variance  a  . 

a 

Such  a  sequence  of  random  variables  is  often  referred  to  as  "white 
noise."  The  time  series  which  describes  the  behavior  of  the  system  will 
be  denoted  by  z-^,  z^,  .  .  .  z^  and  will  be  regarded  as  a  realization  of 
a  jointly  distributed  random  variable. 

The  Box-Jenkins  method  is  based  on  the  idea  that  a  time  series  in 
which  successive  values  are  highly  dependent  can  be  regarded  as  the 
result  of  passing  white  r.oise  through  a  linear  filter.  This  linear 
filtering  operation  considers  the  behavior  of  the  system  to  be  a  linear 
combination  of  the  present  and  past  shocks, 

zt  =  y  +  at  +  Vt-i  +  Vt-2  +  ’  1  *  +  Vt-N  (1) 

where  y  is  the  level  of  the  process,  are  weights  applied  to  the 

previous  shocks,  the  sum  of  which  must  converge,  and  N  is  the  number  of 

observations  in  the  time  series.  If  this  model  were  used  to  represent 

the  behavior  of  a  time  series,  it  would  be  necessary  to  estimate  N  +  2 

2 

parameters  (y,  ...  o  ),  This  exceeds  the  number  of  observa¬ 

tions.  However  it  is  possible  to  simplify  the  above  model  when  it  can 


92 


be  assumed  that  the  mean  and  variance/covariance  matrix  (of  the  joint 
probability  distribution  of  the  stochastic  process  which  generated  the 
time  series)  exist  and  are  invariant  with  respect  to  time.  A  time 
series  generated  by  a  stochastic  process  for  which  the  mean  and  variance/ 
covariance  matrix  exists  and  is  invariant  with  respect  to  time  is  re¬ 
ferred  to  as  a  stationary  time  series. 

In  order  to  explain  why  it  is  necessary  to  require  that  the  time 
series  be  stationary,  it  will  be  helpful  to  first  briefly  describe  the 
process  used  to  establish  a  simplified  model.  Each  of  the  steps 
involved  in  this  process  as  well  as  methods  for  identifying  and  coping 
with  nonstationary  time  series  will  be  discussed  in  greater  detail  in 
later  sections. 

The  Box-Jenkins  process  for  simplifying  the  model  in  Equation  1 
makes  use  of  the  autocorrelation  function  (acf)  and  the  partial  auto-, 
correlation  function  (pacf)  as  the  means  for  tentatively  identifying 
which  of  a  large  family  of  stochastic  models  best  accounts  for  the 
behavior  of  a  given  time  series.  The  acf  measures  the  degree  of 
association  or  mutual  dependence  between  values  of  the  same  time  series 

separated  by  lags  1,  2,  -  time  periods.  The  £th  term  in  the  acf 

represents  the  product  moment  correlation  coefficient  between  values  of 
the  time  series  separated  by  i  time  periods.  The  pacf  measures  the 
degree  of  association  between  values  of  the  time  series  separated  by 

lags  1,  2,  -  time  periods  where  the  effects  of  all  lags  less  than 

the  lag  of  interest  are  removed.  The  £th  term  in  the  pacf  represents 
the  product  moment  correlation  coefficient  between  values  of  the  time 
separated  by  i  time  period  with  the  effects  of  the  correlation  between 
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values  separated  by  lags  1,  2,  - ,  £-1  removed.  The  acf  and  pacf 

describe  how  well  a  system  remembers  its  past. 

The  model  simplification  process  for  doing  this  first  involves 
calculating  the  acf  and  pacf  for  the  observed  time  series.  The  shapes 
of  the  sample  estimates  of  the  acf  and  the  pacf  are  then  compared  with 
the  shapes  of  the  theoretical  acf's  and  the  pacf's  derived  from  the 
family  of  stochastic  processes.  The  parameters  of  the  linear  filter 
whose  theoretical  acf  and  pacf  best  match  the  sample  acf  and  pacf  are 
estimated  by  a  nonlinear  least  squares  technique.  The  adequacy  of  the 
fitted  model  is  then  checked  by  a  procedure  which  necessitates  calcula¬ 
ting  the  autocorrelation  function  of  the  residuals.  If  the  fit  is  not 
adequate  the  autocorrelation  function  is  used  to  determine  how  the  model 
should  be  modified.  Once  an  adequate  model  has  been  found,  it  is  used 
to  forecast  future  values  of  the  time  series. 

In  order  for  the  above  approach  to  work,  it  is  necessary  for  the 
first  and  second  moments  of  the  Joint  probability  density  function  of 
a  sequence  of  observations  generated  by  the  stochastic  process  under 
consideration  to  be  invariant  with  respect  to  time.  The  requirement 
that  the  second  moment  be  invariant  results  from  the  fact  that  the 
autocorrelation  function  is  defined  in  terms  of  the  autocovariances 
and  the  variance: 

_  Yi 

Pi  ~  y—  for  i  =  1,  2,  .  .  . 

T0 
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where 


Yq  =  variance  of  process 

Yi  =  autocovariance  between  two  observations  separated  by  i  time 
periods 

=  autocorrelation  between  two  observations  separated  by  i  time 
units . 

Since  both  yq  and  Yi  are  defined  in  terms  of  the  mean  it  follows 

that  the  first  moment  of  the  time  series  must  also  be  invariant.  A 
process  for  which  the  first  two  moments  of  its  joint  probability  density 
function  are  invariant  with  respect  to  time  will  be  referred  to  as  a 
stationary  process. 

Nonstationari ty  will  always  result  when  the  mean  of  a  process 
changes  over  time.  This  can  occur  when  either  the  level  of  the  process 
abruptly  changes  or  the  level  of  the  process  is  subject  to  a  trend. 
Figure  10  illustrates  a  situation  where  the  level  of  the  process 
suddenly  shifts.  Figure  11  illustrates  a  situation  where  the  level 
of  the  process  is  subject  to  a  trend.  When  this  occurs  the  variance 
and  autocovariances  of  the  process  are  not  meaningful  since  they  are 
defined  in  terms  of  the  mean.  A  nonstable  mean  is  not  the  only  way 
nonstationarity  can  manifest  itself.  Even  when  the  mean  is  stable  it 
is  still  possible  for  the  variance  of  the  process  to  change  over  time 
as  illustrated  by  Figure  12. 
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Figure  10.  Nonstationaity  is  caused  by 

the  sudden  shift  in  the  level 
of  the  process 
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Figure  11.  Nonstationarity  is  caused  by  a 
trend  in  the  process  level. 
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Figure  12.  Example  of  a  Time  Series  which  is  Nonstationary 
due  to  the  Fact  that  the  Variance  is  Heterogeneous 

Nonstationarity  also  manifests  itself  when  the  autocovariance 
structure  of  the  process  changes.  This  situation  is  illustrated  by 
Figure  13. 

If  the  stochastic  process  generating  the  time  series  is  nonstation¬ 
ary,  then  it  will  not  be  possible  to  calculate  the  autocorrelation 
function.  When  this  is  the  case  the  Box-Jenkins  method  for  developing 
a  model  for  describing  the  process  is  not  applicable  since  it  uses  the 
acf  and  pacf  to  identify  a  tentative  model  and  to  check  the  adequacy  of 
the  fitted  model. 


Figure  13.  Example  of  a  Time  Series  which  is 
Nonstationary  due  to  a  .  hange  in  the  Autocovariance 
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All  models  considered  by  the  Box-Jenkins  method  must  be  stationary. 
The  task  of  verifying  that  each  of  these  models  is  stationary  can  be 
simplified  by  exploiting  the  fact  that  all  models  considered  by  the 
Box-Jenkins  method  are  special  cases  of  the  linear  filter  process.  Thus, 
if  the  linear  filter  process  can  be  shown  to  be  stationary,  then  all  of 
its  special  cases  will  also  be  stationary.  To  establish  that  the  linear 
filter  process  is  stationary  it  will  be  necessary  to  show  that  its  mean 
and  variance /covariance  matrix  are  independent  of  time. 

The  mean  of  the  linear  filter  process  described  in  Equation  1 
given  by: 

E(zt)  =  y  +  E(at  +  ^a^  +  *2*t_2  *  •  •  •  ) 

Since  the  a  '  s  represent  independent  random  drawings  from  a  Normal  dis¬ 
tribution  with  mean  zero  and  variance  (J2  it  follows  that 

a 

E(zt)  =  y  (2) 

The  variance  of  the  process  is  defined  as: 

y0  =  E  (zt  -  E(zt)}2  =  E  {at  +  +  .  .  .}2 

=  E  {a2  +  \Jj2  a^  ^2+  .  .  .}  +  E  (cross  product  terms) 

Since  the  a^'s  are  independent  random  variables  the  expected  values  of 
the  cross  product  terms  will  be  zero.  Therefore, 

Y0  =  E{a2  +  ^at_12  +  .  .  .} 

Since  the  a  ' s  are  random  selections  from  a  probability  distribution 

+  • 

with  zero  mean  it  follows  that 
E(at2)  =aa2 

Therefore,  the  variance  of  the  linear  filter  process  described  in 
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Equation  1 


can  be  expressed  as 


Y0 

where 


(3) 


^0  “  1 


The  covariance  between  z  and  z  .  can  be  derived  in  a  similar  fashion 

t  t-j 


Yj  =  E{zt  "  E(zt)}  {zt-j  -  E(zt-j)} 


where 


-  E«\  +  Vt-i  *  •  •  •>  (Vj  *  Vt-j-i  * 
'  E<<Vt-jS)  *  (Vj*iVw!)  *  •  •  •> +  E 

=  a  2  .1  i*> 

a  i=0  i  i+J 


.  .  .)} 

(cross  product  terms) 

(A) 


V 


0 


1 


From  Equations  2,  3,  and  4  it  can  be  seen  that  the  mean,  variance, 
and  covariances  are  independent  of  time.  Therefore,  it  follows  that  the 
linear  filter  process  is  stationary.  The  following  discussion  will 
utilize  the  back  shift  operator  B  defined  such  that 


z 


t-k 


to  simplify  notation. 

A  special  case  of  the  linear  filter  model  is  the  moving  average  (MA) 
model  which  postulates  that  the  current  value  of  a  time  series  is  a 
linear  function  of  its  previous  disturbances.  The  general  form  of  an 

m 

MA  model  of  order  q  is: 

Zt  =  at  -  0 -  0 2at-2  "  •  •  •  "  6qat-q  +  V  ^ 

where 
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y  =  mean  of  MA  process. 

MA(q)  will  be  used  to  denote  a  q~h  order  MA  model.  The  backshift  oper¬ 
ator  B  can  be  used  to  obtain  a  more  compact  expression  for  MA(q) 

z.  =  (i  -  e.B  -  e^B2  -  e  Bq)a  +  y  (f>) 

t  1  d.  q  t 

or 

z^  =  0(B)at  +  y 

where 

0(B)  =  (1  -  0  B  -  0  B2  -  ...  -  0  Bq) 

12  q 

Another  special  case  of  the  linear  filter  model  is  the  autoregres¬ 
sive  model  (AR)  which  postulates  that  the  current  value  of  a  time  series 
is  a  linear  function  of  all  its  past  values.  The  general  form  of  an  AR 
model  of  order  p  is: 

zt  =  ^lzt-l  +  ^2Zt-2  +  •  •  •  +  ^pzt_p  +  at  +  5  (7) 

where 

6  =  mean  of  AR  process. 

A  more  compact  expression  for  AR(p)  can  be  obtained  by  utilizing  the 
backshift  operator.  The  symbol  AR(p)  will  be  used  to  denote  a  p^h 
order  AR  model 

(1  -  (Jj-jB  -  4>2B2  -  ...  -  <t>pBp)zt  =  at  +  6  <8> 

or 

d>(B)z  =  a,  +  6 
Y  t  t 

where 

♦  (B)  =  (1  -  <p,B  -  <J>  B2  BP) 

1  2  p 
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For  every  finite  order  AR  model  it  can  be  shown  that  there  exists 
a  corresponding  infinite  order  MA  model  (Chatfield,  1975).  The  strategy 


It  can  also  be  shown  that  for  every  finite  order  MA  model  there  exists 
a  corresponding  infinite  order  AR  model  (Nelson,  1973).  The  strategy 
for  doing  this  will  be  illustrated  for  the  case  of  anMA(l)  model 

zt  =  at  -  e1at_1  +  y  U2) 

The  term  a^  ^  can  be  eliminated  from  an  MA(l)  model  by  substituting 

zt-l  =  at-l  "  0lat-2  +  y 
in  Equation  k.22  to  obtain 

\  =  at  ~  0]_zt_i  -  9i2at-2  +  V^1  +  ei^ 

Next  the  term  a  0  can  be  eliminated  from  an  MA(l)  model  by  substituting 

zt-2  =  at-2  "  0lat-3  +  y 
to  obtain 

zt  =  at  -  e1zt_1  -  0122t_2  ~  ei3&t_3  +  y(!  +  0!  +  e.^2) 

Proceeding  in  this  fashion  yields 


00  i  u 

zt  =  at  “  til  ei  zt-i  + 


which  is  an  AR(°°)  model  in  which 

6i  =  \i 


(13) 
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6  = 


1  -  61 

Thus,  the  MA  model  and  the  AR  model  represent  alternative  ways  for 
describing  a  time  series.  In  a  particular  situation  the  model  which 
explains  the  situation  with  the  smallest  number  of  parameters  will  be 
preferred.  The  AR  model  will  be  better  for  situations  in  which  the 
effect  of  a  disturbance  on  future  observations  gradually  dies  out.  The 
MA  model  will  be  better  for  situations  in  which  the  effect  of  a  distur¬ 
bance  ends  abruptly  after  a  specified  number  of  observations.  In  situ¬ 
ations  between  these  two  extremes  a  mixed  auto-regressive-moving  average 
(ARMA)  model  will  provide  the  simplest  explanation. 

The  general  form  of  the  ARMA  process  is  given  by 

+  6  +  a.  -  0na 


=  +  •  •  •  +  4>pZt_p  ■  -  ■  -  "i^t-l 


-  6  a+  (14) 

q  t-q 


where 


=  value  of  the  observed  time  series 


zi  =  1 

a  =  i^*1  value  of  the  disturbance 
6  =  mean  of  the  ARMA  process 

ARMA(p,q)  will  be  used  to  denote  an  ARMA  model  with  p  autoregres¬ 
sive  and  q  moving  overage  terms.  The  backshift  operator  can  be  used  to 
obtain  a  more  compact  expression  for  the  ARMA  model 

<j>(B)z  =  0(B)a  +  <5  <15> 

T> 

or 


=  _fiil 


<(>  (B 


at  +  6 


where 
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0(B)  -  (1  -  4)1B  -  02B2  ....  -  <J^BP) 

6(B)  =  (1  -  9. B  -  60B2  -  .  .  ,  ■  0  B?) 

12  q 

Through  the  use  of  an  ARMA  model  it  is  often  possible  to  describe 
a  time  series  vith  fewer  parameters  than  would  be  the  case  if  either  an 
AR  model  or  an  MA  model  were  used  alone.  However,  it  also  includes  both 
the  AR  and  MA  models  as  special  cases.  It,  therefore,  represents  the 
most  general  model  for  describing  a  nonseasonal  stationary  time  series. 
However,  most  economic  time  series  are  both  seasonal  and  nonstationary. 

The  strategy  for  dealing  with  nonstationary  time  series  is  to 
transform  the  original  data  in  a  manner  which  will  remove  the  nonstat- 
ionarity.  Differencing  the  original  data  represents  the  primary  means 
for  removing  the  aspect  of  nonstationarity  caused  by  changes  in  the 
mean  level  of  the  process.  In  practice  it  is  seldom  necessary  to  go 
beyond  the  second  difference.  If  w^  is  defined  as  the  sequence  of  first 
differences 

vt  =  zt  "  zt-l  (16) 

then  the  ARMA  model  applied  to  the  differenced  data  becomes: 

\  *  Vt-i +  •  •  •  *  Vt-p  *  “t  ♦  -  •  •  •  -  Vt..,  <17) 

Note  that  the  process  of  taking  the  first  difference  eliminates  the 
constant  term  6  from  the  ARMA  model  given  by  Equation  14  .  The  back¬ 
ward  shift  operator  B  can  be  used  to  obtain  a  more  compapt  expression 


for  Equation  17 

♦(B)  (1  -  B)dz  =  6 ( B ) a 

L>  vi 
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where  d  represents  the  level  of  differencing  required  to  remove  the  non- 
stationarity  and  the  remaining  terms  have  been  defined  previously.  The 
model  in  Equation  18  describes  the  change  in  the  level  of  the  time 
series.  The  level  at  any  given  time  will  be  equal  to.  the  sum  of  all 
past  changes.  Since  summation  is  a  discrete  counterpart  of  integration 
the  model  in  Equation  18  is  referred  to  as  integrated  autoregressive- 
moving  average  (ARIMA)  model.  An  ARIMA  model  in  which  an  ARMA(p,q)  model 
is  used  to  describe  the  d'*'*1  difference  of  the  data  is  referred  to  by 
ARIMA (p,d,q) . 

In  economic  time  series,  nonstationarity  often  arises  from  a 
tendency  of  the  variance  of  a  time  series  to  increase  over  time.  This 
source  of  nonstationarity  can  often  be  removed  by  working  with  logs  of 
the  original  data. 

The  seasonality  encountered  in  economic  time  series  is  usually 
multiplicative  in  nature.  To  cope  with  multiplicative  seasonality  Box 
and  Jenkins  recommended  using  the  following  model: 

4>  ( B )  ( 1  -  B)d$(B)(l  -  BS)Dz.  =  0  (B)0(Bs  )a  (19) 

L  t 

where 

s  =  period  of  seasonality 

<J>(B)  =  (l  -  4>  B  -  4>„B2  -  .  .  .  -  <f>  Bp)  =  nonseasonal  AR  operator 

1  d  P 

$>(B)  =  (l  -  ^B5-  s  -  ...  -  4>pB^)  =  seasonal  AR  operator 

0(B)  =  (l  -  0..B  -  0^B2  -  .  .  .  -  0  B*l)  =  nonseasonal  MA  operator 

12  q  • 

0(B)  =  (1  -  O^B3  -  0oB2S  -  ...  -  ®qB^)  =  seasonal  MA  operator 

d  =  order  of  nonseasonal  differencing 

D  =  order  of  seasonal  differencing 
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(l  -  B)d  =  nonseasonal  difference  operator 

(l  -  Bs)d  =  seasonal  difference  operator. 

The  symbol  ARIMA(p,P,d,D,q,Q)  will  be  used  to  represent  the  model 
described  by  Equation  19. 

The  ARIMA(p,P,d,D,q,Q)  model  encompasses  most  forecasting  models 
as  -  special  cases.  These  special  cases  include  the  fol¬ 
lowing  forecasting  methods:  (a)  single  exponential  smoothing;  (b) 
double  exponential  smoothing;  (c)  triple  exponential  smoothing;  (d) 
Winters  method;  (e)  adaptive  smoothing;  and  (f)  X-ll  variant  of  the 
Census  Method  II  Seasonal  Adjustment  Program  developed  by  the  United 
States  Bureau  of  the  Census. 

The  flow  chart  in  Figure  9  provides  an  overview  of  the  process  used 
to  determine  which  model  within  the  ARIMA(p,P,d,D,q,Q)  family  is  best 
for  forecasting  the  behavior  of  a  particular  time  series. 

Identification.  The  identification  process  is  concerned  with 
tentatively  selecting  a  model  from  the  ARIMA(p ,P ,d  ,D  ,q  ,Q)  family  of 
models.  This  involves  specifying  values  for  p,  P,  d,  D,  q,  and  Q.  The 
first  step  is  to  determine  whether  or  not  the  series  is  stationary.  A 
rough  idea  as  to  whether  or  not  a  particular  series  is  stationary  can  be 
obtained  by  examining  a  plot  of  the  raw  data.  Nonstationarity  will  be 
evidenced  by  changes  in  the  mean  and/or  variance  of  the  process  over 
time.  Trends  will  also  indicate  nonstationarity.  The  apf  can  also  be 
used  to  determine  whether  or  not  a  time  series  is  stationary.  If  a 
series  is  nonstationary  the  autocorrelations  will  die  out  slowly  (Box  & 
Jenkins,  1970,  p.  17^-175;  200-201).  The  nonstationarity  which  exists 
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Figure  14.  Flow  Chart  of  the  Process  for  Determining  the  Specific 
Model  within  the  ARIMA(p,P,d,D,q,Q)  Family  of  Models 
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in  a  series  can  be  either  seasonal  or  nonseasonal  or  both,  ie: 

1.  Nonseasonal  nonstationarity  is  evidenced  when  the  autocorrela¬ 
tions  at  successive  lags  die  out  slowly  as  in  Figure  15. 

2.  Seasonal  nonstationarity  is  evidenced  when  there  is  a  cyclic 

'  ) 
pattern  in  the  autocorrelations  as  in  Figure  16. 

3.  Seasonal  nonstationarity  is  also  evidenced  when  the  autocor- 

* 

relations  separated  by  the  seasonal  period  die  out  slowly  as 
in  Figure  17. 

4.  Both  seasonal  and  nonseasonal  nonstationarity  are  evidenced 
when  1  and  either  2  or  3  occur. 

Nonstationarity  can  be  reduced  by  differencing  the  original  data. 
There  are  two  basic  types  of  differences:  (a)  nonseasonal  differences 
which  are  obtained  by  subtracting  observations  separated  by  one  time 
period;  and  (b)  seasonal  differences  which  are  obtained  by  subtracting 
observations  separated  by  the  length  of  the  seasonal  period. 

It  is  very  important  to  correctly  identify  the  orders  of  the  non¬ 
seasonal  differences  and  the  seasonal  differences.  Undifferencing 
leaves  dependencies  in  the  data  which  should  be  removed,  while  over¬ 
differencing  introduces  dependencies  in  the  data  which  should  not  be 
there.  The  following  can  be  used  as  guidelines  for  identifying  the 
proper  orders  of  nonseasonal  and  seasonal  differences: 

1.  The  variance  for  a  correctly  differenced  series  will  be  lower 
than  the  variances  for  the  series  which  have  be^n  either 
overdifferenced  or  underdifferenced. 

2.  The  value  of  the  Box-Pierce  (1970)  test  statistic 

k  2 

Q  -  T  I  r/  (a) 
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where 


T  =  total  number  of  observations  minus  the  maximum  order  of 
differencing 

K  =  number  of  values  of  sample  acf  that  have  been  calculated 

2  ' 

r^~(a)  =  sample  acf  of  the  difference  time  series  will  be  smaller 
for  the  correctly  differenced  series  than  the  values 
associated  with  series  that  have  been  either  overdifferenced 
or  underdifferenced. 

The  next  step  in  the  identification  process  after  obtaining  a  sta¬ 
tionary  time  series  is  to  determine  the  orders  of  the  nonseasonal  and 
seasonal  autoregressive  and  moving  average  operators.  This  will  be 
accomplished  by  comparing  the  shapes  of  the  sample  acf  and  pacf  with  the 
shapes  of  the  theoretical  acf's  and  pacf's  corresponding  to  the  various 
models  in  the  ARIMA(p,  P,  d,  D,  q,  Q)  family  of  models.  The  model  whose 
theoretical  acf  and  pacf  best  match  the  shapes  of  the  sample  acf  and 
pacf  will  be  tentatively  selected  as  the  model  to  be  fitted  to  the 
observed  data. 

Since  in  practice  it  is  seldom  necessary  for  any  of  the  operators 
in  the  ARIMA(p,P,d,D,q,Q)  model  to  be  higher  than  the  second  order.  Box 
and  Jenkins  (1970,  Chapter  3)  have  analyzed  the  shapes  of  the  theoreti¬ 
cal  acf's  and  pacf's  for  these  models.  This  analysis  provided  the  basis 


111 


for  constructing  the  charts  in  Appendix  A.  These  charts  can  be  used  to 
facilitate  the  process  of  determining  the  tentative  orders  for  the  non- 
seasonal  and  seasonal  autoregressive  and  moving  average  operators. 

These  charts  should  be  used  as  follows: 

1.  For  nonseasonal  models  compare  the  shape  of  the  estimated  acf 
of  the  stationary  series  with  the  shapes  of  the  theoretical  acf's  shown 
in  the  charts  and  select  the  model  (or  models)  whose  theoretical  acf 
best  matches  the  shape  of  the  sample  acf. 

2.  For  a  model  which  is  completely  seasonal  ignore  the  values  of 
the  sample  acf  at  lags  which  do  not  correspond  with  the  seasonal  period 
and  then  determine  the  model  whose  theoretical  acf  best  matches  the 
shape  of  the  relevant  portions  of  the  sample  acf. 

3.  For  a  model  which  has  both  seasonal  and  nonseasonal  operators: 
(a)  Disregard  the  significant  seasonal  correlations  in  the  sample  acf 
and  apply  step  1  to  the  relevant  portion  of  the  sample  acf.  (b)  Dis¬ 
regard  the  significant  nonseasonal  correlations  in  the  sample  acf  and 
apply  step  2  to  the  relevant  portion  of  the  sample  acf. 

The  process  for  determining  the  tentative  orders  for  the  seasonal 
and  nonseasonal  autoregressive  and  moving  average  operators  places 
primary  emphasis  on  matching  the  sample  acf  with  the  theoretical  acf 
corresponding  to  one  of  the  ARIMA(p ,P,d,D,q,Q)  models.  If  acf's  of  two 
or  more  of  these  models  appear  to  match  a  sample  acf  equally  well,  then 
the  pacf's  should  be  compared  to  determine  which  is  most  appropriate. 
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The  reason  for  placing  primary  emphasis  on  the  acf  is  that  the  procedure 
for  estimating  the  pacf  "becomes  very  sensitive  to  rounding  errors  and 
should  not  be  used  if  the  values  of  the  parameters  are  close  to  the 
nonstationary  boundaries"  (Box  and  Jenkins,  1970,  p.  65). 

The  autocorrelations  calculated  from  the  stationary  time  series 
are  only  estimates  of  the  actual  correlations  and  are  therefore  subject 
to  sampling  error.  To  test  the  significance  of  an  individual  autocor¬ 
relation  coefficient  Box  and  Jenkins  (1970,  p.  3^—3 5 )  recommended  using 
an  approximate  expression  for  the  variance  of  the  autocorrelation  coef¬ 
ficient  developed  by  Bartlett  (19^+6), 

V(r  )  =  JL.  U  +  2  .1  p.2} 

J  II  i=l  l 

where 

N  =  number  of  observations  in  stationary  time  series 

rj  =  sample  estimate  of  autocorrelation  coefficient  for  a  pair  of 
observations  separated  by  J  periods 

=  theoretical  autocorrelation  coefficient  for  a  pair  of  obser¬ 
vations  separated  by  i  periods 

V(r  )  =  variance  of  r 
J  J 

q  =  separation  beyond  which  theoretical  autocorrelation  coefficients 
are  zero 

If  the  mean  of  a  stationary  time  series  is  significantly  different 
from  zero,  then  it  will  be  necessary  to  add  a  constant  term  to  the 
ARIMA(p,P,d,D,q,Q) .  This  situation  will  arise  when  there  is  a  persistent 
trend  in  the  undifferenced  time  series. 
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Estimation.  Once  the  tentative  form  of  the  model  has  been  decided 
upon  the  next  step  is  to  estimate  the  parameter  values  in  the  model . 

The  least  squares  criterion  is  used  to  determine  the  best  estimates 
of  the  parameters  of  the  ARIMA(p,P,d,D,q,Q)  and  their  standard  devia¬ 
tions.  When  the  a^'s  are  Normally  distributed  the  least  squares 
estimators  of 


4>  =  (<t1#  4>2, 
$  =  *2, 
e  =  (61,  e2, 
0  =  (ei,  e2. 


in  the  ARIMA(p,P,d,D,q,Q)  model 

<J>(B)(1  -  B)d$  (B)  (1  -  Bs )Dz  =  0(B)0  (Bs)a  (20) 

L>  X> 

where 


s  =  period  of  seasonality 

<)>(B)  =  (l  -  <f>  B  -  4>  B2  —  .  .  .  —  4>  B^)  =  nonseasonal  AB  operator 

1  2  p 

$(B)  =  (l  -  $^BS-  <t>2BJ  -  ...  -  4>pB^)  =  seasonal  AR  operator 

0(B)  =  (l  -  S^B  -  SgB2  -  ...  -  =  nonseasonal  MA  operator 

0(B)  =  (l  -  0^BS  -  02B2S  -  ...  -  0q®^)  =  seasonal  MA  operator 

d  =  order  of  nonseasonal  differencing 


D  =  order  of  seasonal  differencing 

(l-B)^  =  nonseasonal  difference  operator 
(l-Bs)D=  seasonal  difference  operator 


will  closely  approximate  the  maximum  likelihood  estimates  of  these  para¬ 
meters  . 
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The  expression  for  is  obtained  by  rewriting  Equation  20  as 
follows : 

at  =  e_1  (B)0-1  Cb)4>(B)4>(B)  (1  -  B)d(l  -  Bs)Dzt  ''‘lj 

Note  that  a^  is  a  nonlinear  function  of  the  parameters  of  the  ARIMA(p,P, 
d,D,q,Q)  model.  Hill  climbing  techniques  are  used  to  determine  the 
parameter  values  which  minimize  the  sum  of  squares. 

The  computational  effort  required  to  solve  this  nonlinear  least 
squares  problem  can  be  reduced  by  starting  the  search  procedure  at 
parameter  values  which  are  close  to  their  true  values.  Table  1  summa¬ 
rizes  the  procedure  for  obtaining  preliminary  estimates  of  the  para¬ 
meters  for  nonseasonal  models . 

Diagnostic  testing.  The  adequacy  of  a  model  will  be  indicated  when 
the  following  conditions  are  satisfied:  (a)  the  residuals  are  not 
significantly  different  from  white  noise;  and  (b)  all  of  the  parameters 
of  the  model  are  significantly  different  from  zero. 

Box  and  Pierce  (1970)  have  developed  a  test  for  determining  whether 
or  not  a  sample  acf  is  significantly  different  from  white  noise.  This 
test  utilizes  the  following  test  statistic 

Q  =  T  E  r  2(a) 
j-1  J 

where 

T  =  total  number  of  observations  minus  the  maximum  order  of 

differencing 

K  =  number  of  values  of  sample  acf  that  have  been  calculated 

Tj2(a)  =  sample  acf  of  the  estimated  residual  series. 
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TABLE  1 


Summary  of  Procedures  for  Obtaining  Preliminary  Estimates 
of  the  Parameters  of  Nonseasonal  Models 


p 

d 

<1 

Preliminary  estimate  of  parameter 

1 

any 

positive 

integer 

0 

4^  =  p.^  subject  to  -1  <  <f>  <  1 

0 

any 

positive 

integer 

1 

See  Table  A,  p.  517-518,  in  Box 
and  Jenkins  (1970)  to  obtain 
estimate  of  0^. 

2 

any 

positive 

integer 

0 

See  Chart  B,  p.  518,  in  Box  and 
Jenkins  (1970)  to  obtain  estimates 
of  4>1  and  <t>2. 

0 

any 

positive 

integer 

2 

See  Chart  C,  p.  519,  in  Box  and 
Jenkins  (1970)  to  obtain  estimates 
of  0^  and 

1 

any 

positive 

integer 

i 

See  Chart  D,  p.  520,  in  Box  and 
Jenkins  (1970)  to  obtain  estimates 
of  4>1  and  0^, 

any 

positive 

integer 

any 

positive 

integer 

any 

positive 

integer 

See  Appendix  A6.2,  p.  201-205,  in 
Box  and  Jenkins  (1970). 

Q  is  approximately  chi-square  distributed  with  K-p-q  degrees  of  freedom 
where  p  and  q  represent  the  number  of  autoregressive  and  moving  average 
parameters  in  the  model  under  consideration. 

If  the  above  test  indicates  that  the  residual  series  is  signifi¬ 
cantly  different  from  white  noise,  then  the  sample  acf  of  the  residuals 
should  be  analyzed  to  determine  how  the  model  should  be  improved 

Another  means  of  assessing  the  adequacy  of  the  fitted  model  is  to 
test  the  significance  of  the  individual  terms  in  the  sample  acf.  This 
test  is  based  on  the  fact  that  the  distribution  of  the  estimated  auto¬ 
correlations  of  residuals  is  approximately  Normal  with  mean  zero  and 
variance  equal  to  the  reciprical  of  the  sample  size  (Box  &  Jenkins, 

1970,  p.  290). 

Using  more  parameters  in  a  model  than  necessary  is  referred  to  a-s 
overfitting.  This  should  be  avoided  since  the  presence  of  the  unneces¬ 
sary  parameters  will  probably  bias  the  estimates  of  the  essential  para¬ 
meters.  Unnecessary  parameters  will  be  indicated  by  the  fact  that  the 
1  -  a  confidence  limits  for  the  parameter  estimate  include  the  value 
zero  (The  procedure  for  establishing  these  confidence  limits  is 
described  in  Box  and  Jenkins  (1970,  p.  22U-231).). 

An  Overview  of  the  Box-Jenkins  Transfer  Function  Model 
The  Box-Jenkins  univariate  method  utilizes  only  the  historical 
values  of  a  time  series  in  forecasting  future  behavior.  Often  there 
will  exist  one  or  more  variables  which  influence  the  behavior  of  the 
variable  being  forecast.  These  causal  variables  have  the  potential  for 
increasing  forecast  accuracy.  To  exploit  this  potential  Box  and  Jenkins 
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(1970)  developed  a  process  for  building  forecasting  models  which  utilizes 
both  the  historical  values  of  the  series  being  forecast  and  leading 
variables  to  forecast  future  behavior.  This  model  building  process  is 
referred  to  as  the  transfer  function  method. 

The  transfer  function  model  visualizes  the  leading  variables  as 
inputs  to  a  process  which  produce  the  variable  being  forecast  as  an 
output.  The  objective  of  the  transfer  function  model  is  not  to  discover 
a  model  which  describes  in  great  detail  the  causal  mechanism  whereby  the 
input  is  transformed  into  the  output.  Instead,  the  transfer  function 
model  regards  the  transformation  as  a  black  box  and  seeks  only  to  dis¬ 
cover  which  of  a  large  family  of  mathematical  models  best  describe  how 
a  change  in  the  input  affects  the  output  over  time.  The*  function  which 
describes  how  the  input,  during  a  particular  period,  affects  the  output 
during  future  time  periods  is  referred  to  as  the  impulse  response 
function. 

The  impulse  response  functions  found  in  practice  take  on  a  wide 
variety  of  shapes.  Some  will  follow  an  oscillatory  pattern  while  others 
will  be  oscillation-free.  For  those  following  an  oscillatory  pattern 

some  will  follow  a  damped  sine  wave  pattern  while  others  will  follow  a 
pattern  in  which  the  effects  in  successive  periods  will  alternate  in 
sign.  For  those  that  are  oscillation-free  some  will  be  unimodal  while 
others  will  be  monotone.  In  some  situations  the  effects- will  decay 
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gradually  over  time.  In  other  situations  the  effects  will  persist  for 
a  period  of  time  and  then  end  abruptly.  The  delay  between  the  change 
in  the  input  and  its  affect  on  the  output  will  vary  from  one  situation 


to  another.  The  variety  of  possible  shapes  which  the  input  response 
function  can  take  on  in  practice  necessitates  that  the  family  of  math¬ 
ematical  models  from  which  the  theoretical  impulse  response  is  selected 
must  be  extremely  flexible.  The  required  flexibility  can  be  found  in  a 
model  of  the  type: 

fi(B)y  =  w(B)Bbxt  (22) 

where 

x^  =  input 

yt  =  output 

B  =  backshift  operator 

b  =  delay 

6(B)  =  (l  -  6-^B  -  ...  -  6^Br)  =  output  lag  operator 
U)(B)  =  (cl)0  -  io^B  -  ...  -  wgBs)  =  input  lag  operator 

This  model  can  be  rewritten  as 


=  6(B)xt-b 


where 


UQ  +  6  B  +  62  B  + 


119 


The  polynomial  operator  v (B)  represents  the  transfer  function 
which  summarizes  the  dynamic  structure  of  the  effect  transferred  from 
the  input  sequence  to  the  output  sequence.  The  values  vq’  vi’  V2’ 
represent  the  changes  in  output  expected  from  a  one-time-only  unit 
change  in  input  and  are  referred  to  as  the  impulse  response  coefficients. 
The  only  restriction  on  v (B)  is  that  if  the  input  is  held  at  a  fixed 
level  Xf  then  the  output  eventually  should  settle  at  an  equilibrium 
level  y  .  This  requirement  often  is  not  fulfilled  by  the  data  in 
the  original  form  because  of  trends  and  seasonality.  Fortunately,  such 
problems  can  often  be  removed  by  appropriately  differencing  the  data. 

The  input  lag  operator  w(B)  describes  the  magnitude  of  the  more 
immediate  effects  of  the  input.  The  output  lag  operator  describes  the 
duration  and  pattern  of  decay  for  these  effects. 

Identification.  The  process  for  tentatively  identifying  the  trans¬ 
fer  function  makes  use  of  the  cross  correlation  function  (ccf)  which 
measures  the  degree  of  association  between  the  present  value  of  a  given 
time  series  variable  and  past,  present,  and  future  values  of  another 
time-series  variable.  The  ccf  is  defined  as 

.  Wk) 

a  a 

x  y 

=  cross  correlation  between  x  and  y  at  lag  k 
=  cross  covariance  between  x  and  y  at  lag  k 
=  standard  deviation  of  x 
=  standard  deviation  of  y 


p  (k) 
Mxy 


where 


Vk) 

YXy(k) 
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The  procedure  is  to,  first,  estimate  the  sample  ccf  for  the  observed 
data  and  then  to  find  a  theoretical  transfer  function  model  whose  ccf 
closely  matches  the  shape  of  the  sample  ccf.  However,  the  sample  ccf 
can  not  be  estimated  directly  from  the  input  and  output  series  if  the 
input  series  is  highly  autocorrelated.  This  results  from  the  fact  that 
autocorrelation  in  the  input  series  induces  autocorrelation  in  the  ccf. 
This  will  distort  the  shape  of  the  sample  ccf  which  will  make  it  more 
difficult  to  correctly  identify  the  transfer  function.  Fortunately, 
this  problem  can  be  circumvented  by  reducing  the  input  series  to  white 
noise.  This  can  be  accomplished  by  utilizing  the  Box-Jenkins  univariate 
model  described  in  the  preceding  section  to  fit  an  autoregressive  inte¬ 
grated  moving  average  (ARIMA)  model  to  the  input  series.  The  process  of 
reducing  the  input  series  to  white  noise  is  referred  to  as  prewhitening 
the  input.  In  order  to  preserve  the  relationship  between  input  and  out¬ 
put  it  will  be  necessary  to  apply  the  same  model  to  the  output  series. 

The  result  of  this  operation  is  referred  to  as  the  prewhitened  output 
series.  Since  the  same  transformation  is  applied  to  both  the  input  and 
output  series,  the  transfer  function  model  which  fits  the  prewhitened 
series  will  also  fit  the  original  series. 

Box  and  Jenkins  (1970,  p.  3*47-353)  recommend  the  following  procedure 
for  using  the  sample  ccf  tc  tentatively  identify  the  transfer  function 
model : 
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1.  The  value  of  b,  the  delay,  will  be  equal  to  the  number  of  lags 
before  the  first  non-zero  cross  correlation  coefficient  is  encountered. 


r  = 


2.  The  value  of  r  will  be  determined  by  the  type  pattern  followed 
by  the  remaining  values 

0  if  there  are  no  values  following  a  pattern 

1  if  the  remaining  values  decay  exponentially 

2  if  the  remaining  values  follow  a  damped  sine  wave. 

•3.  The  value  of  s  can  be  obtained  as  follows: 

bp  -  (b  +  l)  for  np  >  1  for  r  >  1 

np  -  1  r  =  0 


s  = 


where 


bp  =  number  of  values  up  to  and  including  the  value  at  which  the 
pattern  begins 

np  =  number  of  values  not  following  a  pattern 

In  practice  the  input  will  almost  always  be  able  to  explain  only  a 
portion  of  the  variation  in  output.  The  inability  of  the  input  to  fully 
account  for  all  of  the  variation  in  output  results  from  the  presence  of 
measurement  errors  and  the  influence  of  other  causal  variables.  The 
residual  series 


nt  =  yt  -  V(B)xt_b 


where 

A 

V(B)  =  preliminary  estimate  of  transfer  function 
is  used  to  estimate  the  noise  series.  The  Box-Jenkins  univariate  model 
is  then  used  to  develop  an  ABIMA  model  for  the  noise  series 


4>(B)(1  -  B)dnt  =  0(B)at 


or 


e(B) 

(1  -  B)d4>(B) 


Thus  the  combined  transfer  function-noise  model  is 

v  =  a(B)  Bbx  + _ all} _ . 

*  5(B)  t  (2  _  B)d4,(B)  t 

The  process  for  identifying  a  transfer  function  noise  model  is 
shown  in  the  flow  chart  in  Figure  14 . 

Estimation.  Once  the  tentative  forms  of  the  transfer  function  and 
the  noise  model  have  been  determined,  the  next  step  is  to  simultaneously 
estimate  all  of  the  parameters  in  the  combined  model 


=  —IB) 


B  x.  + 
5(B)  ^ 


(1  -  B)d<(>(B) 


where 


6(B)  =  (1  -  61B  - 


-  6rBr)  =  output  lag  operator 


U)(B)  =  (wQ  -  co^B  -  ...  -  WsBS)  =  input  lag  operator 

<£(B)  =  (l  -  4^B  -  (Jj^B2  -  ...  -  4>pBP)  =  nonseasonal  AR  operator 

9(B)  =  (l  -  0^B  -  0OB2  -  ...  -  BqB^)  =  nonseasonal  MA  operator 

d  =  order  of  nonseasonal  differencing 
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(l-B)^  =  nonseasonal  difference  operator 
b  =  delay 

The  least  squares  criterion  is  used  to  determine  the  best  estimates 
of  the  parameters  of  the  combined  model  and  their  standard  deviations. 
For  the  case  when  the  at's  are  Normally  distributed,  the  least  squares 
estimators  of 

(jj  ~  ((a)q  *  *  •  •  •  *  ) 

. i) 

1  2  r 

4>  =  4>2>  •  •  •  >  ♦p) 

6  =  (e1,  02,  ....  0  ) 

in  the  combined  model  will  closely  approximate  the  maximum  likelihood 
estimators  of  these  parameters  (Box  &  Jenkins,  1970,  p.  388). 

The  values  of  the  a  are  estimated  by  the  following  three  stage  ' 
process: 

1.  The  estimated  output,  y  ,  from  the  transfer  function  model  is 
obtained  by  solving  Equation  22  for  y  and  then  substituting  the 
actual  input,  x^,  into  the  resulting  expression 

9t  =  6"1(B)cj(B)Bbxt 

2.  The  noise  series,  ,  represents  the  difference  between  the 
estimated  output  and  the  actual  output 


3.  The  values  of  the  a  series  are  obtained  by  solving  the  noise 
model  for  a  and  then  passing  the  n  series  through  the  resulting 
expression 
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Note  that  a^  is  a  nonlinear  function  of  the  parameters  of  the  transfer 
function  and  noise  models. 

Hill  climbing  techniques  are  used  to  determine  the  values  of  the 
parameters  in  the  combined  transfer  function  noise  model  which  minimize 
the  sum  of  squares.  The  computational  effort  required  to  solve  this  non¬ 
linear  least  squares  problem  can  be  reduced  by  starting  the  search 
procedure  at  parameter  values  which  are  close  to  their  true  values.  Box 
and  Jenkins  (1970,  p.  378-380)  developed  a  procedure  for  obtaining  pre¬ 
liminary  estimates  of  the  parameters  of  the  transfer  function.  These 
should  be  used  in  conjunction  with  their  procedures  for  obtaining  pre¬ 
liminary  estimates  of  the  parameters  of  the  noise  model  which  are 
summarized 

Diagnostic  testing.  The  process  for  assessing  the  adequacy  of  the 
model  utilizes  both  the  acf  of  the  estimated  residuals  and  the  ccf  where 
the  prewhitened  input  series  is  lagged  relative  to  the  estimated  residual 
series.  When  the  acf  of  the  estimated  residuals  differs  significantly 
from  that  of  white  noise,  a  problem  in  the  noise  model  is  indicated.  A 
ccf  which  differs  significantly  from  that  which  would  be  obtained  when 
two  white  noise  series  are  cross  correlated  indicates  that  the  transfer 
function  model  is  deficient. 

The  procedure  for  determining  whether  or  not  the  acf  differs 
significantly  from  white  noise  has  already  been  described 

in  connection  with  diagnostic  testing  of  a  univariate  model. 


Box  and  Pierce  (1970,  p.  395)  developed  a  test  for  determining 
whether  or  not  a  ccf  differs  significantly  from  that  of  two  white  noise 
series.  This  procedure  utilizes  the  test  statistic 

Q  =  m  £  r2^(k) 

k=0  aa 

where 

m  =  n-u-p 

n  =  number  of  observations  available  for  analysis 

u  =  max  {r,  s  +  b} 

p  =  number  of  moving  average  terms  in  noise  model 

b  =  delay 

r  =  number  of  output  lag  parameters 

s  -  number  of  input  lag  parameters  minus  one. 

Q  is  approximately  chi-square  distributed  with  K+l-(r+s+l)  degrees 
of  freedom.  If  the  above  test  indicates  that  the  residual  series  is 
significantly  cross  correlated  with  the  lagged  values  of  the  prewhitened 
input  series,  then  the  sample  ccf  should  be  used  to  determine  how  the 
model  should  be  improved. 

Box  and  Jenkins  (1970,  p.  39^-395)  also  recommended  testing  the 
significance  of  the  individual  terms  in  the  sample  ccf  as  a  means  for 
determining  whether  or  not  a  ccf  differs  significantly  from  that  of  two 
white  noise  series.  This  test  is  based  on  the  assumption  that  the  dis¬ 
tribution  of  the  estimated  cross  correlations  is  approximately  Normal 
with  mean  zero  and  variance  equal  to  the  reciprocal  of  the  sample  size. 
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Using  more  parameters  in  a  model  than  necessary  is  referred  to  as 
overfitting.  This  should  be  avoided  since  the  presence  of  the  unneces¬ 
sary  parameters  will  probably  bias  the  estimates  of  the  essential  param¬ 
eters.  Unnecessary  parameters  will  be  indicated  when  the  1  -  a  confi¬ 
dence  limits  for  a  parameter  estimate  includes  the  value  zero  (the 
procedure  for  establishing  these  confidence  limits  is  described  in  Box 
and  Jenkins,  1970,  p.  391). 
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APPENDIX  A 

Theoretical  Shapes  of  ACF's  and  PACF's 
for  Selected  ARMA(p,q)  Models 
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APPENDIX  A 

THEORETICAL  SHAPES  OF  ACF's  AND  PACF'd 
FOR  SELECTED  ARMA ( p , q )  MODELS 


Theoretical  Autocorrelation  Functions 

(p=0,  q=l) 


dj  >  o  e2  <  o 
-1  0+1  -10  +1 


(p=0,  q*2) 

e2  >  0  0-l  >  o  <  o  e1  <  o 

e2  >  o  e2  <  o  e2  >  o  e2  <  o 


-1.0+1  -10  +1  -1  0+1  -10  +1 


(p=l,  q=0) 


4>2  >  0  <f>  <  o 

-1  0+1  -10  +1 


i 


(p=2,  q=0) 


CD 

M 

A 

O 

V 

r—i 

CD 

0 

e-L  >  o 

e-L  >  o 

O 

V 

C\J 

CD 

e2  > 

0 

o2  <  o 

CD 

ro 

V 

o 

■1  0  +1 

-1  0 

+1 

-1  0  +1 

-1  0  +1 

(p=0,  q=l) 


<  0  4>i  >  0 

-1  O+l  -10  +1 


(p=0,  q=2) 


<  0  <  0 

4>2  <  0  <j>2  >  0 

-1  0+1  -10  +1 


<t>!  >  0 
4*2  <  o 
-1  0  +1 


4>i  >  o 
4>2  >  o 
-1  0  +1 


133 


(P=l,  q=l) 


e1>(j)1>o  4>1>e1>o 

and 

<t>T<  0<61 
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PREFACE 


The  Industrial  Engineering  and  Operations  Research  Department's 
activities  were  divided  into  three  phases.  The  first  phase,  which  is 
reported  in  "A  Diagnosis  of  the  Workload  Variation  Problem  in  Ship¬ 
building,"  was  designed  to  identify  specific  problems  for  further 
analysis.  This  analysis  revealed  the  need  for  the  following: 

(1)  a  system  for  analyzing  how  variations  in  shipyard  workloads  affect 
shipbuilding  cost; 

(2)  a  shipyard  planning  system  which  could  be  utilized  by  shipyards  to 
minimize  the  cost  of  adjusting  to  workload  variations. 

Pht’se  2,  which  is  described  in  "A  Framework  for  Analyzing  How  Variations 
in  Shipyard  Workloads  Impact  Shipbuilding  Cost,"  focused  on  developing 
the  framework  of  a  system  for  developing  transfer  function  models  which 
can  be  used  to  develop  the  difference  equations  which  describe  how 
variations  in  shipyard  workloads  affect  shipbuilding  cost.  Phase  3, 
which  is  the  subject  of  this  report,  focused  on  developing  the  framework 
of  a  shipyard  planning  system  which  shipyards  could  use  to  minimize  the 
cost  of  adjusting  to  workload  variations. 
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I .  INTRODUCTION 


Building  a  ship  is  a  large  scale  endeavor  which  typically  requires 
three  or  four  years  to  complete.  During  this  time  the  services  of  the 
construction  trades  listed  in  Table  1  will  be  used. 

These  trades  will  be  responsible  for  producing,  erecting,  installing, 
assembling,  and  testing  an  enormous  number  of  parts  as  indicated  in  Table 
2. 

Frisch  (1976)  observed  that  the  most  appropriate  method  for  producing 
any  of  the  above  parts  or  assembling  any  number  of  them  into  subsystems 
will  depend  on  the  following  factors:  (1)  the  quantity  of  the  item  to  be 
produced  In  a  single  run;  (2)  the  weight  of  the  items  being  produced  or 
assembled;  and  (3)  the  complexity  of  the  items  or  subsystem  which  is  a 
function  of  the  tolerances  which  must  be  maintained  in  production  and/or 
assembly.  He  then  defined  seven  different  levels  of  complexity  which  are 
illustrated  in  Table  3.  He  then  constructed  a  table  (Table  4) 
to  contrast  the  complexity  of  a  ship  with  that  of  a  watch  and  a  car. 

A  wide  variety  of  equipment  will  be  used  in  performing  the  numerous 
tasks  involved  in  constructing  a  ship.  These  tasks  are  typically  classified 
into  three  categories:  (1)  hull  construction;  (2)  outfitting;  and  (3) 
machinery  installation.  Table  5  shows  the  major  processes  involved  in 
hull  construction  and  the  major  types  of  equipment  used  during  each  process. 

In  addition  there  are  numerous  major  items  which  are  installed  during 
outfitting.  These  are  listed  in  Table  6. 

The  major  items  of  equipment  installed  during  machinery  installation 
are  listed  in  Table  7.  Numerous  types  of  specialized  equipment  will  be 
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Table  1 


CONSTRUCTION  TRADES  IN  SHIPBUILDING 


Automotive  Operator  (all  types)  Machine  Operator 

Blacksmith  Machinist  (inside-outside) 

Boilermaker  Maintenance  Man 


Calker  and  Chipper 
Carpenter 

Crane  Operator  (all  types) 
Diesel  Mechanic 
Electrician  (inside-outside) 
Electronics  Mechanic 


Material  Chaser 
Nuclear  (all  types) 

Painter  (all  types) 

Pipe  Coverer  and  Insulator 
Pipefitter 

Quality  Control  Technician 


Engine  Locomotive  Operator 
(all  types) 

Engine  and  Pump  Operator 

Flame  Cutter  (Burner) 

Heat  Treating  Mechanic 

(Bulkhead  Straightman) 

Hooker-on 

Hydraulic  Mechanic 

Instrument  Calibrator 

Joiner 

Laborer,  Shipyard 
Lead  Burner 
Lof tsman 


Rigger 
Sandblaster 
Sheetmetal  Worker 
Shipf itter 
Tractor  Operator 
Tool  Clerk 
Tool  Repairman 
Ventilation  Man 
Warehouse  Man 
Welder 

X-Ray  Radiographer 


Source:  Mark  Battle  Associates,  Inc.,  Shipbuilding  Manpower  Stud^ 
Executive  Summary,  March,  1974. 
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Table  2 


NUMBER  OF  UNITS  OR  PIECES  IN  BUILDING  A  SHIP 


Number 

or 

of  Units 

Pieces 

10,000  ton 

dwt  cargo  60,000  ton 
ship  dwt  tanker 

Number  of  Sizes 

or 

Classes  of  Items 

Structural  steel  (cut 
pieces) 

25,000 

40,000 

NA 

Sheet  steel  (cut  pieces) 

1,500 

2,500 

NA 

Sections  of  pipe 

2,500 

3,500 

NA 

Electric  cable 
(various) 

100,000 

200,000 

NA 

Propulsion  engines, 
other  machinery 
units,  motors, 
pumps,  fans,  etc. 

200 

350 

NA 

Various  purchased 
material  items 

22,500 

60,000 

NA 

Stock  and  supply  items 
by  sizes  and  types 
(bolts,  nuts,  stock 
fittings,  pipe 
flanges,  etc.) 

1,500-2,500 

Source:  Mack-Forlist,  Daniel  M.  and  Arthur  Newman,  The  Conversion  of 
Shipbuilding  From  Military  To  Civilian  Markets,  Praeger 
Publishers,  New  York,  1970,  p.  36. 
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j  COMPLEXITY  ASSOCIATED  TOLERANCE 
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Table  A 


COMPLEXITY  COMPOSITION 


EXAMPLES 

COMPLEXITY 

Watch1 

Car2 

Ship3 

GROUP 

Weight 

7o 

Weight 

% 

Weight 

1 

1 

(EXTRA  FINE) 

92 

1 

1 

2 

(FINE) 

8 

5 

2 

3 

(MEDIUM  FINE) 

— 

A5 

13 

A 

(MEDIUM) 

— 

A 

20 

5 

(MEDIUM  COARSE) 

— 

A5 

15 

6 

(COARSE) 

— 

— 

A5 

7 

(EXTRA  COARSE) 

— 

— 

A 

HlECHANICAL  WRISTWATCHES 
2passenger  CARS 


3GENERAL  CARGO  SHIPS  -  COMMERCIAL 


Source:  Frisch,  F.  A.  P.,  "Production  and  Construction: 

A  Comparison  of  Concepts  in  Shipbuilding  and 
Other  Industries",  Paper  presented  at  the  Fifth 
Annual  DOD  Procurement  Research  Conference, 
November  17-18,  1976,  p.  62. 


146 


HULL  CONSTRUCTION 
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trades  are  incompatiable  with  one  another.  For  example,  none  of  the 
other  trades  can  work  in  the  area  below  welders.  Another  example,  no 

trade  can  work  in  the  same  area  as  painters.  Even  if  there  are  no  prece¬ 
dence  relationships  between  two  or  more  work  packages  and  work  space  is 
not  a  limiting  factor,  these  work  packages  cannot  be  scheduled  simulta¬ 
neously  if  there  is  insufficient  manpower. 

The  rates  at  which  the  various  trades  are  used  in  constructing  a 
given  ship  is  not  constant  over  time.  During  the  early  stages  of  con¬ 
struction,  trades  involved  in  constructing  the  hull,  such  as  welders, 
shipfitters,  riggers,  loftsmen,  hooker-on  men,  and  flame  cutters,  will 
dominate  the  work  force  assigned  to  a  given  ship.  During  the  latter 
stages  of  construction,  trades  involved  with  outfitting,  such  as  elec¬ 
tricians,  pipefitters,  sheetmetal  workers,  and  carpenters,  will  dominate 
the  workforce.  The  number  of  workers  in  each  trade  that  will  be  needed 
during  each  scheduling  period  is  difficult  to  determine.  This  results 
from  the  fact  that  "in  ship  construction,  the  model  is  the  product  and 
everything  which  cannot  be  predetermined  with  paper  studies  must  be  tried 
out  during  production  in  the  form  of  rework,  change  orders,  or  something 
else"  (Frisch,  1976,  p.  40).  This  ad  hoc  approach  to  ship  construction 
makes  it  virtually  impossible  to  determine  priori  all  of  the  activities 
which  will  be  necessary  to  construct  a  ship.  Consequently,  shipbuilders 
are  not  able  to  establish  standard  methods  to  the  extent  that  this  is 
done  for  products  that  are  produced  continuously  or  in  large  lot  sizes. 
The  lack  of  standard  methods  makes  it  difficult  to  develop  accurate  time 
standards  which  form  the  basis  for  estimating  manpower  requirements. 
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Nevertheless,  not  having  the  right  number  of  workers  in  each  trade 


during  each  time  period  can  be  extremely  costly.  If  too  few  workers  are 
assigned  to  a  bottleneck  activity,  then  the  ship  will  not  be  completed 
on  time.  The  penalty  for  not  completing  a  ship  on  time  consists  of: 

(1)  penalty  payments  for  late  delivery  explicitly  specified  in  the  con¬ 
tract;  and  (2)  opportunity  cost.  Typically,  contracts  for  Navy  ships  do 
not  explicitly  contain  late  delivery  penalty  clauses.  The  opportunity 
cost  for  a  50  million  dollar  ship  have  been  estimated  to  be  at  least 
$20,000  per  day  (Corporate-Tech  Planning,  Inc.,  1978).  A  shipyard  will  typically 
take  one  or  more  of  the  following  actions  when  a  planned  delivery  date 
is  in  jeopardy:  (1)  schedule  overtime;  (2)  increase  crew  sizes;  (3)  post¬ 
pone  work  until  after  launch;  and  (4)  increase  management  involvement  in 
day-to-day  operation  (Corporate-Tech  Planning,  Inc.  ,  1978).  On  the  other  hand 
having  too  many  workers  would  create  idle  time  and  increase  labor  cost 
unnecessarily . 

Periodically  wage  rates  for  the  various  trades  and  cost  of  the  vari¬ 
ous  raw  materials  increase  as  a  result  of  inflation.  However,  these  rate 
and  cost  increases  do  not  necessarily  occur  at  the  same  time.  The  times 
at  which  these  increases  occur  can  be  determined  from  union  contracts  for 
the  various  trades  and  by  examining  the  times  in  the  past  at  which  pre¬ 
vious  price  increases  have  occurred  for  the  various  raw  materials.  Note 
that  the  rate  of  increase  for  each  of  the  trades  and  for  each  of  the  raw 
materials  will  not  necessarily  be  the  same.  Under  a  fixed  cost  contract, 
inflationary  trends  provide  an  incentive  to  the  shipyard  to  complete  a  ship 
as  soon  as  possible.  However,  attempting  to  do  this  can  result  in  inter¬ 
trade  and  intratrade  interferences. 
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TABLE  6 

MAJOR  ITEMS  INSTALLED  DURING  OUTFITTING 


1.  Doors 

,  a.  Watertight 
b.  Non-watertight 

2.  Airports,  port  lights  and  windows 

3.  Manholes 

4.  Tonnage  openings  (optional) 

5.  Cargo  hatch  covers 

6.  Bulwarks 

7.  Cargo  handling  systems 

8.  Surface  grating 

9.  Ladders 

10.  Anchors  and  anchor  handling 

11.  Mooring  lines 

12.  Lifting  equipment 


13.  Sparrings 

14.  Deck  covering 

15.  Insulation 

-  Fire  insulation 

-  Sound  insulation 

-  Refrigeration  insulation 

-  Comfort  insulation 

16.  Support  wiring 

17.  Electrical  wiring 

18.  Mooring  and  towing  gears 

19.  Rudders 

20.  Steering  gear  and  electrical  equipment 
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TABLE  6  (continued) 


21.  Piping  systems 

a.  Bilge  and  ballast 

b.  Fuel  oil  transfer 

c.  Fresh  water  system 

d.  Fire  systems 

-  sprinklers 

-  CO2  system,  foam 

e.  Cargo  system 

f.  Sanitary  system 

22.  Ventilation  systems 

a.  Cargo  hold  ventilation  and  dehumidification 

b.  Living  and  working  area  ventilation  and  dehumidification 

23.  Heating  and  air  conditioning 

a.  Cargo  heating  (typically  for  tankers) 

b.  Living  and  working  areas 

24.  Life  saving  equipment 

a.  Life  boats,  rafts  and  floats 

b.  Life  boat  mechanism  and  winches 

c.  Life  jackets 

25.  Navigation  equipment 

a.  Signal  equipment 

b.  Radio  direction  finder. 

c.  Navigation  lights 

d.  Sounding  gear 

e.  Compasses 

f.  Radar 

g.  Loran 

26.  Weapon  systems  (for  combatant  ships) 
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TABLE  7 


MAJOR  ITEMS  OF  EQUIPMENT  INSTALLED  DURING  MACHINERY  INSTALLATION 

PRESERVATION  AND  MAINTENANCE 

1.  Painting 

a.  Anti  fouling  marine  coating 

b.  Anti  corrosion 

2.  Sanding 

3.  Blasting 

PROPULSION  MACHINERY 

1.  Steam  turbine 

2.  Diesel 

3.  Reduction  gears 
A.  Thrust  bearing 

5.  Propeller  shaft 

6.  Propeller 

7.  Crankshaft 

AUXILIARY  MACHINERY 

1 .  Pumps 

2.  Boilers 

3.  Electric  diesels  and/or  turbines 

4.  Evaporators 
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used  by  each  of  the  many  trades  involved  in  outfitting  and  machinery 
installation. 

From  the  above  it  is  obvious  that  constructing  a  ship  is  an  enor¬ 
mously  complicated  process  which  involves  a  staggering  number  of  diverse 
activities.  The  activities  for  each  trade  are  grouped  in  categories  on 
the  basis  of  geographic  proximity.  These  categories  are  known  as  work 
packages.  Typically  the  hull  erection  process  and  the  outfitting  process 
will  be  divided  into  approximately  2,500  work  packages  for  each  process. 

A  typical  work  package  requires  approximately  three  months  and  500  man¬ 
hours  to  complete  (Corporate-Tech  Planning,  Inc.,  1978). 

These  work  packages  cannot  be  scheduled  independently  of  one  another. 
In  many  instances  there  will  be  either  a  preferred  sequence  or  a  manda¬ 
tory  sequence  in  which  certain  of  the  work  packages  should  be  performed. 
When  a  preferred  sequence  is  not  followed  unnecessary  work  will  result. 

For  example,  suppose  that  the  crane  available  at  the  erection  site  has 
the  capacity  to  lift  the  main  engine, whereas  the  one  at  the  outfitting 
berth  does  not.  The  main  engine  can  be  lifted  aboard  as  an  assembled 
unit  prior  to  launching.  After  launching  it  will  be  necessary  to  partial¬ 
ly  disassemble  the  main  engine,  lift  the  components  separately,  and  re¬ 
assemble  the  main  engine  on  board.  A  mandatory  sequence  results  from  the 
fact  that  it  is  physically  impossible  to  perform  certain  activities  prior 
to  others.  For  example,  the  electrical  cable  in  a  compartment  cannot  be 
installed  before  the  compartment  has  been  built.  Another  factor  which 
precludes  performing  several  work  packages  in  parallel  is  limited  work 
space.  Interference  results  when  too  many  workers  from  either  the  same 
or  different  trades  are  assigned  to  a  confined  workspace.  Also  certain 
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Certain  costs  will  be  incurred  when  the  size  of  the  workforce  is 
either  increased  or  decreased.  Increasing  the  size  of  the  work  force  can 
involve  the  following  activities:  (1)  recruiting;  (2)  interviewing;  (3) 
testing;  (4)  performing  medical  exams;  (5)  placing;  and  (6)  training. 

The  following  often  result  when  the  size  of  the  workforce  is  decreased: 

(1)  exit  interviews; (2)  separation  payments;  (3)  increased 
unemployement  insurance  premiums;  and  (4)  a  bad  image  in  the  local  labor 
market.  Hancock  (1971  p.  7-113)  stated  that  the  cost  for  "the  leaving 
of  one  person  and  tne  hiring  of  a  new  one  are  from  $600  to  $2,000".  These 
costs  do  not  account  for  the  hyper  inflation  which  has  transpired  since 
1971.  If  an  adjustment  for  this  is  made  then  the  current  dollar  equivalent  of 
these  costs  would  be  approximately  $1,050  to  $3,500. 

The  typical  shipyard  has  several  ships  in  varying  stages  of  comple¬ 
tion  at  any  given  time.  These  ships  each  compete  for  a  share  of  the  ship¬ 
yard's  manpower  and  facilities.  This  competition  takes  place  during  each 
and  every  period  in  the  planning  horizon.  Management's  task  is  to  deter¬ 
mine  how  the  size  of  each  trade's  workforce  should  be  changed  from  period 
to  period  so  as  to  minimize  the  sum  of  the  following:  (a)  lateness 
penalty  cost;  (b)  cost  of  changing  the  size  of  the  workforce;  (c)  labor 
cost;  and  (d)  materials  cost.  This  type  of  problem  is  referred  to  as  the 
multi-resource/multi-project  planning  problem. 

The  multi-project/multi-resource  planning  problem  is  a  very  real  one 
in  shipyards.  This  was  evidenced  by  a  study  in  which  manufacturing  execu¬ 
tives  in  the  major  U.  S.  shipyards  were  asked  to  rank  the  causes  of  delays 
and  cost  overruns  in  shipbuilding  (Terry,  Green,  and  Magnuson,  1979). 


1 54 


The  most  important  finding  of  this  survey  was  that  62.5  percent  of  the 
shipyards  engaged  in  Naval  construction  ranked  the  unavailability  of 
skilled  labor  as  the  most  important  cause  of  delays  and  cost  overruns; 

37.5  percent  of  the  shipyards  ranked  improper  sequence  of  assigned  labor 
as  the  most  important  cause. 

It  is  obvious  from  the  above  that  the  multi-resource/multi-project 
(mrmp)  planning  problems  which  shipyards  must  face  are  extremely  complex. 
Furthermore,  the  cost  of  making  wrong  decisions  can  be  extremely  high. 

This  would  seem  to  suggest  that  shipyards  should  be  using  the  mrmp  plan¬ 
ning  models.  In  fact  the  first  such  model  was  described  in  a  paper  en¬ 
titled  "Multi-Ship,  Multi-Shop,  Workload  Smoothing  Program"  (Levy, 
Thompson,  and  Wiest,  1963).  However,  a  survey  (Terry,  Green  and  Magnuson, 
1979)  to  determine  the  types  of  planning  tools  being  used  by  shipbuilders 
revealed  that  not  one  out  of  a  sample  of  21  U.  S.  shipbuilders  were  using 
any  type  of  mrmp  planning  model  (Levy,  Thompson,  and  Wiest,  1963; 
Lambourne,  1963;  Wiest,  1967;  Pritsker,  Watters,  and  Wolfe,  1969;  Dar-El, 
Behmoaram,  and  Tur,  1978).  Furthermore,  none  were  using  any  of  the 
multi-resource/single  project  models  that  had  been  developed.  Davis 
(1973)  reviewed  the  state-of-the-art  for  such  models  up  through  mid-1972. 
In  this  paper  Davis  referred  to  138  papers  which  he  deemed  to  be  relevant 
to  the  problem  of  project  scheduling  under  resource  constraints.  Since 
this  review  the  following  papers  have  appeared  in  the  open  literature: 
Davis  and  Patterson,  1973;  Dar-El  and  Tur,  1977;  Dar-El -and  Tur,  1978; 
Dar-El,  Behmoaram,  and  Tur,  1978;  Horroelen,  1973;  Patterson,  1973;  and 
Patterson  and  Huber,  1974.  The  survey  described  by  Terry,  Green,  and 
Magnuson  also  revealed  only  half  of  the  shipyards  were  using  PERT  or  CPM. 
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None  were  using  the  probabilistic  features  of  PERT  or  the  time-cost 
features  of  CPM. 

The  irony  of  this  situation  compels  the  question:  why  are  ship¬ 
yards  not  using  project  planning  tools  which  prima  facie  evidence 
suggests  that  they  desperately  need?  The  factors  responsible  for  this 
unexpected  result  are  explored  in  the  following  section. 
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II.  FACTORS  RESPONSIBLE  FOR  LACK  OF  USE  OF  MULTI-RESOURCE/MULTI-PROJECT 

PLANNING  MODELS 

Several  factors  were  identified  which  were  suspected  to  be  respon¬ 
sible  for  the  failure  of  U.  S.  shipbuilders  to  use  multi-resource/multi- 
project  (mrmp)  planning  models.  They  are:  (1)  shipyard  executives  are 
not  aware  of  the  existence  of  such  models;  (2)  the  models  are  too  com¬ 
plicated  for  shipyard  personnel;  (3)  the  models  require  data  that  could 
not  be  obtained  at  a  reasonable  cost;  and  (4)  the  models  fail  to  recog¬ 
nize  the  fact  that  decisions  concerning  enterprise  goals  and  the  acquisi¬ 
tion  of  resources  are  typically  made  at  a  higher  management  level  than 
decisions  concerned  with  the  allocation  of  existing  resources. 

It  is  quite  likely  that  many  shipyard  executives  are  not  aware  of 
the  existence  of  mrmp  planning  models.  These  models  are  not  usually 
covered  in  introductory  operations  research/management  science  courses. 

Most  shipyard  executives  are  not  likely  to  be  familiar  with  the  journals 
in  which  these  models  have  been  described.  Furthermore,  most  and  per¬ 
haps  all  of  the  articles  on  mrmp  planning  have  not  been  written  in  a 
style  which  would  appeal  to  and  be  understood  by  most  shipyard  executives. 

However,  the  biggest  potential  barrier  to  the  use  of  such  models 
arises  from  their  data  requirements.  Current  mrmp  planning  models  re¬ 
quire  that  the  user  specify  the  precedence  relationships  between  the 
various  activities  involved  in  building  a  ship.  The  size  and  scope  of  a 
ship's  construction  would  result  in  a  combinational  problem  too  large  and 
too  expensive  for  any  single  heuristic  mrmp  planning  model  to  handle 
(Levy,  Thompson,  and  Wiest,  1963;  Lambourne,  1963;  Wiest,  1967;  and  Dar-El, 
Behmoaram,  and  Tur;  1978). 
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Another  barrier  to  implementing  such  models  results  from  the  fact 


that  constructing  a  ship  is  an  evolutionary  process  (Frisch,  1976)  . 

Under  such  conditions  it  is  not  possible  to  identify  a  priori  all  of  the 
important  activities  which  will  be  involved  in  constructing  a  ship.  In 
some  situations  it  will  not  be  possible  to  specify,  at  planning  time, 
which  of  several  mutually  exclusive  approaches  for  accomplishing  a  given 
job  will  be  optimal.  Circumstances  can  arise  which  can  make  it  necessary 
to  perform  certain  activities  which  were  in  no  way  expected  at  planning 
time.  Conversely  activities  which  were  believed  to  be  absolutely  neces¬ 
sary  at  planning  time  can  be  found  to  be  unnecessary  as  the  project  pro¬ 
gresses.  These  difficulties  are  compounded  by  the  fact  that  "missing  the 
date  on  one  part  of  the  construction  cycle  can  cause  severe  cascade  effects 
on  other  parts"  (Corporate-Tech  Planning,  Inc.,  1978,  p.  D-9).  For  ex¬ 
ample,  failure  to  get  the  main  engine  installed  prior  to  launch  can  neces¬ 
sitate  partially  disassembling  it  so  that  its  components  can  be  lifted  by 
the  crane  at  the  outfit  pier.  This  creates  the  necessity  for  assembling 
it.  on  the  ship  where  work  space  is  limited  (Frisch,  1979).  The  problem  of 
specifying  precedence  relationships  is  further  complicated  by  the  inter¬ 
action  between  ships  as  evidenced  by  the  following: 

Since  all  major  shipyards  have  more  than  one  ship 
under  construction  at  one  time,  the  planning  of  successive 
ships  is  constrained  by  the  status  of  ships  already  under 
construction.  Build  methods  are  influenced  by  other  ship 
work,  whether  planned  or  already  started.  (Corporate-Tech 
Planning  Inc.,  1978,  p.  D-3) 

This  situation  clearly  suggests  that  the  precedence  relationships  tend  to 
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evolve  over  time  as  work  on  the  various  ships  in  the  yard  progresses. 

This  indicates  that  frequent  reruns  of  the  mrmp  planning  model  will  be 
necessary  which  will  magnify  the  difficulty  of  an  already  difficult  prohit 
Nevertheless  it  is  clear  that  shipyards  need  a  more  effective  method 
for  coping  with  the  mrmp  planning  problem.  Evidence  to  support  this 
claim  is  provided  by  the  following: 

"Introduction  of  new  or  improved  management  procedures, 
particularly  those  which  improve  ship  construction 
scheduling  and  thus  minimize  the  present  queueing  pro¬ 
blems  that  occur  when  certain  shipyard  occupations  are 
used  intensely  for  sometime,  only  to  stand  idle  later. 

The  lessening  of  insecurity  about  this  layoff /rehire 
pattern  could  attract  more  workers",  (Mark  Battle 
Associates,  Inc.,  1974,  p.  8). 

This  need  is  particularly  great  in  the  outfitting  area  as  noted  by 
Millen  (1978,  p.  8) 

"Few  shipyards  can  claim  to  have  an  effective  production 
planning  organization  for  the  steelworking  operation  which 
itself  is  responsible  for  at  least  half  of  the  work  content 
of  each  vessel.  Still  fewer  shipyards  achieve  any  effective 
planning  and  coordination  of  the  outfitting  departments  each 
of  which,  because  of  its  relatively  smaller  contribution  and 
through  concern  for  overhead  costs,  is  usually  encouraged  to 
carry  out  its  planning  on  an  informal  basis". 

The  survey  of  shipyard  executives  and  engineers  (Terry,  Green  and  Magnuson 
1979)  revealed  that  most  delays  occurred  in  outfitting.  This  is  not 
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suprising  in  view  of  the  informal  approach  to  outfit  planning.  Frisch 
(1979)  observed: 

"Wien  outfitting  problems  arise  there  is  hardly  ever 
sufficient  time  to  acquire  additional  resources.  The 
problems  typically  have  to  be  solved  with  the  resources 
at  hand". 

The  need  for  a  more  effective  method  for  solving  the  mrmp  planning 
problem  is  further  emphasized  by  the  fact  that  a  sample  of  296  Navy 
ships  built  since  the  early  1950' s  were  finished,  on  the  average,  10.5 
months  late.^  This  represents  a  great  economic  loss  to  both  the  ship¬ 
builders  and  the  nation. 

The  following  section  describes  a  theoretical  framework  which 
provides  the  basis  for  developing  models  for  solving  the  mrmp  planning 
problems  of  the  size  and  difficulty  faced  by  shipbuilders. 


1.  The  basis  for  this  sample  was  ship  progress  curves  maintained  by  the 
Naval  Sea  Systems  Command. 


III.  THEORETICAL  FRAMEWORK  FOR  DECOMPOSING  MULTI-RESOURCE/MULTI-PROJECT 
PROBLEMS  INTO  SOLVABLE  PARTS 

As  seen  in  the  previous  section  the  multi-resource/multi-project 
(mrmp)  planning  problems  which  shipyards  encounter  are  far  too  large  to 
be  solved  by  existing  mrmp  planning  models.  This  section  describes  a 
two  stage  approach  to  decompose  this  problem  into  solvable  parts:  (1) 
decompose  the  total  problem  into  strategic  and  tactical  components;  and 
(2)  decompose  activity  network  into  sub-networks  that  can  be  handled 
separately. 

Decomposition  into  Strategic  and  Tactical  Components 

Anthony  (1965)  provided  a  taxonomy  which  classifies  planning 
decisions  into  two  hierarchical  levels:  strategic  planning  and  tactical 
planning.  This  taxonomy  provides  the  basis  for  the  first  stage  of 
decomposing  the  mrmp  planning  problem  into  manageable  parts. 

Strategic  Planning.  The  purpose  of  strategic  planning  is  to  identify 
the  goals  of  the  enterprise,  to  evaluate  the  internal  strengths  and 
weaknesses  of  the  enterprise,  and  to  monitor  the  environment  to  identify 
problems/opportunities  associated  with  accomplishing  these  goals.  The 
information  obtained  from  monitoring  the  environment  and  evaluating 
internal  strengths  and  weaknesses  typically  give  rise  to  two  types  of 
questions:  (1)  ’lew  to  solve  the  problems  which  the  enterprise  is  likely 

to  face  in  the  future?  (2)  How  to  exploit  the  opportunities  which  are 
likely  to  occur?  Answers  to  these  questions  enable  the  organization  to 
identify  the  resources  needed  to  achieve  its  objectives. 

The  resource  needs  of  the  organization  are  then  translated  into 
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"broad  brush"  plans  for  acquiring  the  needed  resources  and  disposing  of 
those  not  required.  This  plan  focuses  attention  on  those  resources  that 
are  costly  and/or  hove  long  procurement  lead  times.  This  plan  should 
specify  what  the  levels  of  these  major  resources  should  be  from  period 
to  period.  In  order  to  do  this  it  will  be  necessary  to  determine  how 
these  resources  should  be  deployed.  However,  detailed  information  on  how 
the  major  resources  should  be  deployed  is  not  needed  for  the  purposes  of 
formulating  this  plan.  The  appropriate  level  of  detail  on  the  intended 
utilization  of  such  resources  should  correspond  to  that  needed  to  make 
an  intelligent  decision  regarding  whether  or  not  the  resources  should  be 
acquired . 

The  results  of  the  strategic  planning  process  creates  the  environ¬ 
ment  in  which  tactical  decisions  must  be  made.  The  goal  of  strategic 
planning  is  to  create  an  environment  for  tactical  planning  which  is  in 
some  sense  optimal.  However,  strategic  decisions  are  made  in  the  face 
of  great  uncertainty.  As  the  future  unfolds  uncertainty  is  reduced. 
Therefore,  strategic  plans  need  to  be  updated  periodically  to  reflect 
the  impact  of  intervening  events. 

Tactical  Planning.  Tactical  planning  is  concerned  with  developing 
action  plans  for  effectively  utilizing  available  resources  to  accomplish 
short  range  objectives  which  have  been  derived  from  the  strategic  plan. 
Another  function  is  to  plan  for  the  acquisition  and  disposal  of  resources 
as  specified  in  the  strategic  plan.  If  the  tactical  planner  is  unable  to 
meet  the  broad  general  objectives  specified  by  the  strategic  plan,  then 
the  tactical  planner  should  immediately  advise  the  strategic  planners 
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of  this  situation.  There  are  certain  items  of  detailed  information  that 
are  considered  in  the  tactical  planning  process  that  are  either  ignored 
or  dealt  with  in  an  aggregate  form  in  the  strategic  planning  process. 

Examples  of  this  type  of  information  are  detailed  precedence  relation¬ 
ships  and  intercraft  and  intracraft  interference. 

Differences  Between  Strategic  and  Tactical  Planning.  Any  type  of 
planning  must  be  based  on  a  sample  of  what  is  likely  to  happen  in  the 
future.  However,  strategic  planning  and  tactical  planning  differ  in 
terms  of  the  type  of  forecast  required.  Strategic  planning  requires  a 
longer  range  forecast  than  tactical  planning.  Tactical  planning  requires 
a  more  detailed  forecast  than  strategic  planning.  This  suggests  that  the 
data  base  which  supports  tactical  planning  must  be  more  detailed  and 
updated  more  frequently  than  the  one  for  strategic  planning.  The  differ¬ 
ence  in  data  base  requirements  for  strategic  and  tactical  planning 
suggests  that  using  a  monolithic  model  which  simultaneously  solves  both 
the  strategic  planning  and  the  tactical  planning  problems  will  increase 
the  size  of  the  problem  unnecessarily  and  largely  explain  why  current 
algorithms  are  not  being  used.  Furthermore,  it  overloads  both  the 
strategic  and  tactical  planners  with  irrelevant  information. 

The  strategic  planner  is  given  too  much  detail  while  the  planning  horizon 
would  be • too  'long  f or  the  tactical  planner. 

The  strategic  planner  has  two  degrees  of  freedom  which  are  not 
available  to  the  tactical  planner:  (1)  freedom  to  shift* target  completion  dates 
and  (2)  freedom  to  vary  the  levels  of  the  various  resources  from  period  to 
period.  Making  such  decisions  at  the  tactical  or  operational  level  is 
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undesirable  due  to  the  enormous  potential  cost  of  a  wrong  decision  and 
the  limited  purview  of  the  tactical  planner.  However,  information 
indicating  a  need  to  revise  the  strategic  plan  will  not  always  surface 
at  the  strategic  management  level.  When  it  surfaces  at  a  lower  level 
it  should  be  immediately  transmitted  to  the  strategic  level. 

The  inputs  to  the  strategic  planning  process  occur  on  a  more  random 
and  less  predictable  basis  than  do  inputs  to  the  tactical  planning 
process.  This  suggests  that  strategic  planning  should  be  done  on  an 
exception  basis  while  tactical  planning  should  be  done  on  a  periodic 
basis.  The  rational  for  this  practice  is  that  the  time  and  cost  associated 
with  generating  a  strategic  plan  are  too  great  to  be  squandered  unnecess¬ 
arily.  What  is  needed  is  a  more  efficient  and  cost  effective  means  of 
doing  strategic  planning,  one  which  is  more  responsive  to  important 
changes  in  both  ''internal"  and  ’’external”  environments.  Less  important 
changes  should  be  handled  informally  if  the  cost  of  revising  the  strategic 
plan  exceeds  the  benefits  of  doing  so. 

The  mrmp  planning  models  reported  in  the  literature  attempt  to  deal 
with  both  the  strategic  and  tactical  planning  problems  simultaneously, 
but  these  cannot  cope  with  the  enormous  planning  problems  faced  by 
shipbuilders.  These  models  also  fail  to  recognize  that  all  of  the  infor¬ 
mation  used  in  tactical  planning  is  not  necessary  for  strategic  planning. 

Decompose  Activity  Network  into  Sub-networks  which  can  be  Solved 
Independently.  Approximately  5,000  work  packages  (roughly  split  50-50 
between  hull  erection  and  outfitting)  are  involved  in  constructing  the 
typical  ship  (Corporate-Tech  Planning,  Inc.,  1978)  which  utilizes  about 
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40  trades  for  its  completion  (Mark  Battle  Associates,  Inc.,  1974).  The 
average  work  package  requires  about  three  months  to  complete,  but  these 
cannot  be  performed  in  an  arbitrary  sequence;  for  example,  either  a 
mandatory  sequence  or  a  preferred  sequence  may  exist.  A  mandatory 
sequence  exists  when  it  is  physically  impossible  to  perform  certain 
activities  prior  to  certain  prerequisite  activities.  For  example,  the 
main  engine  cannot  be  installed  until  it  has  been  built  and  the  engine 
mounts  have  been  fabricated.  A  preferred  sequence  occurs  when  it  is 
more  cost  effective  to  perform  activities  in  a  certain  sequence,  but 
nevertheless  possible  to  perform  them  in  a  different  sequence.  For 
example,  heavy  piping  may  be  installed  at  a  lower  cost  when  the  blocks 
are  being  assembled  and  work  space  is  less  restricted  than  anytime 
thereafter.  Space  limitations  can  also  prevent  certain  activities  from 
being  performed  simultaneously.  Such  activities  are  said  to  be  incom¬ 
patible.  For  example,  floor  coverers  and  electricians  cannot  work 
effectively  in  the  same  small  compartment. 

An  enormously  large  and  complicated  activity  network  would  be 
necessary  for  describing  the  numerous  precedence  and  incompatability 
relationships  which  exist  between  the  5,000  or  so  work  packages  involved 
in  constructing  a  ship. 

The  typical  shipyard  will  have  several  ships  under  construction  at 
any  one  time.  An  enormously  large  activity  network  will  be  required  for 
each  of  these  ships.  The  number  of  work  packages  necessary  for  the  ships 
in  a  given  shipyard  at  any  given  time  could  run  into  tens  of  thousands. 
Furthermore,  these  activities  will  be  competing  for  the  limited  number 
of  workers  in  each  trade's  workforce. 
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The  tactical  shipyard  planning  problem  is  concerned  with  allocating 
each  trade's  limited  workforce  to  the  various  activities  competing  for 
their  services.  Shipyard  management  would  obviously  like  to  make  this 
allocation  in  the  most  effective  manner.  However,  formulating  a  mathe¬ 
matical  model  of  this  problem  could  result  in  a  combinatorial  problem  of 
immense  proportions. 

The  problem  of  determining  how  to  allocate  each  trade's  workforce  to 
the  various  activities  is  complicated  by  the  great  amount  of  uncertainty 
in  the  construction  process.  This  uncertainty  manifests  itself  in  two 
forms:  (1)  The  activity  duration  times  are  influenced  by  .a  large  number 

of  factors  which  cannot  be  controlled  due  to  the  somewhat  exploratory 
nature  of  shipbuilding  (Frisch,  1976).  This  makes  it  difficult  to 
accurately  estimate  the  durations  for  the  various  activities.  (2)  The 
set  of  activities  necessary  for  constructing  a  ship  cannot  be  accurately 
foreseen.  Factors  beyond  the  control  of  shipyard  management  can  cause 
the  set  of  activities  to  change.  Frisch  (1976,  p.  40)  gives  the  follow¬ 
ing  examples  of  such  activities. 

"a  subcontractor  may  not  be  able  to  deliver;  new  weapon 
developments  must  be  accommodated;  weather  in  open 
building  areas  influences  worker  efficiency". 

Changes  in  the  set  of  activities  necessary  to  construct  a  ship  can  be 
caused  by  decisions  made  by  shipyard  management.  For  example,  the 
original  plan  might  call  for  the  main  engine  to  be  installed  prior  to 
launching.  However,  if  building  the  main  engine  is  delayed  shipyard 
management  might  decide  to  install  it  after  launching. 

Random  variations  in  activity  duration  times  can  be  handled  by 


the  probability  features  of  PERT.  However,  MacCrimmon  and  Ryavec  (1964) 
have  carefully  analyzed  the  assumptions  underlying  this  method.  Their 
analysis  revealed  a  number  of  fundamental  problems  which  are  summarized 
below. 

1.  Beta  distribution  might  be  appropriate. 

2.  When  there  are  several  paths  from  the  origin  node  to  the 
terminal  node,  then  the  duration  of  the  project  is  the  maximum 
of  the  durations  of  each  path.  Typically  the  paths  will  share 
common  resources.  Consequently,  they  will  not  be  independent. 
However,  even  if  they  are  independent  the  duration  of  the 
project  will  not  be  normally  distributed  since  the  maximum  of 
the  two  or  more  normally  distributed  variables  is  not  normally 
distributed  (Clark,  1961).  Furthermore,  even  if  the 
distribution  of  the  duration  of  the  project  is  assumed  to  be  . 
normally  distributed  the  PERT  procedure  yields  biased  estimates 
of  both  the  mean  and  the  variance  (Tippett,  1925).  Specifically, 
the  PERT  procedure  underestimates  the  mean  and  overestimates 

the  variance  of  the  project  duration.  Furthermore,  the  degree 
of  bias  increases  as  the  number  of  paths  through  the  network 
increases. 

3.  The  PERT  formulas  for  calculating  the  mean  and  variance  restrict 
the  values  of  the  shape  parameters  of  the  Beta  distribution  to 
the  extent  that  the  coefficient  skewness  of  the  probability 
distribution  of  activity  duration  can  take  only  three  values 
(-2.5,  2.5,  or  0).  Thus  the  contention  that  the  PERT  procedure 

is  sufficiently  flexible  to  represent  any  probability  distribution 
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of  activity  duration  times  is  misleading. 

As  a  general  rule  the  probability  features  of  PERT  are  not  being  used. 
This  is  evidenced  by  the  fact  that  the  National  Aeronautics  and  Space  Admin¬ 
istration  no  longer  includes  the  probability  of  meeting  schedules  in  the 
output  of  its  NASA-PERT  system  (Moder  and  Phillips,  1970).  Furthermore, 
none  of  the  shipyards  surveyed  in  Terry,  Green,  and  Magnuson  (1979)  re¬ 
ported  using  the  probability  features  of  PERT.  The  most  likely  reason 
for  this  is  that  most  people  have  difficulty  estimating  the  optimistic 
and  the  pessimistic  activity  durations. 

Changes  in  the  set  of  activities  necessary  to  construct  a  ship  caused 
by  factors  beyond  the  control  of  shipyard  management  can  be  modeled  by 
simulation  languages  such  as  GERT  (Pritsker,  1974)  and  TRANSIM  (McMichael 
Orleans,  1975)  which  permit  one  of  several  mutually  exclusive  arcs  leading 
from  a  node  to  be  determined  probabilistically.  However,  such  a  simula¬ 
tion  model  could  not  be  used  to  specify  how  to  allocate  each  crafts  work¬ 
force  to  the  various  work  packages.  Changes  in  the  set  of  activities 
necessary  to  construct  a  ship  caused  by  shipyard  management  decisions  can 
be  modeled  by  the  use  of  Decision  CPM  (Crowston,  W.  0.  and  G.  L.  Thompson) . 
This  technique  utilizes  total  project  cost  as  the  basis  for  selecting  a 
job  method  from  a  set  of  mutually  exclusive  job  methods.  Unfortunately 
this  technique  is  not  able  to  handle  uncertainty  which  results  from 
either  random  activity  durations  or  from  activities  which  cannot  be 
foreseen. 

On  the  basis  of  the  above,  it  would  appear  that  the  uncertainty  in 
the  ship  construction  process  would  make  a  problem  which  is  already 
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hopelessly  difficult  even  more  difficult.  Fortunately,  this  need  not  be 
the  case.  This  uncertainty  can  actually  be  used  to  breakup  the  activity 
networks  for  each  ship  into  a  number  of  subnetworks.  The  rational  for 
doing  this  rests  on  the  belief  that  detailed  plans  should  be  based  on 
hard  facts,  indicating  that  detailed  plans  should  be  developed  for  only 
those  parts  of  the  activity  network  for  which  all  necessary  activities 
of  a  major  nature  can  be  clearly  foreseen. 

The  horizon  for  tactical  planning  should  not  extend  past  the  point 
in  the  activity  network  at  which  the  first  uncertain  activity  occurs. 

For  example,  suppose  there  is  a  question  as  to  whether  or  not  the  main 
engine  can  be  built  in  time  to  be  installed  prior  to  launch.  If  the 
engine  were  ready  prior  to  launch,  it  can  be  installed  as  a  unit.  If 
not  ready  prior  to  launch, then  the  crane  capacity  at  the  outfit  pier 
would  dictate  that  it  be  partially  disassembled ,  lifted  aboard  piece  by 
piece  and  reassembled  in  situ  which  itself  represents  an  additional 
activity.  Furthermore,  workspace  will  be  more  limited  after  launch. 

Thus,  installing  the  main  engine  post- launch  might  necessitate  using  a 
vastly  different  work  method  than  required  for  a  pre-launch  installation. 

In  addition  freeing  up  the  necessary  space  for  assembling  the  main 
engine  on  board  might  necessitate  drastic  changes  in  the  sequence  in  which 
certain  outfitting  activities  are  performed.  Thus,  it  should  be  obvious 
that  the  occurrence  of  uncertain  events  could  drastically  alter  the  nature 
of  the  remainder  of  the  project.  Thus,  it  does  not  seem. prudent  to  develop 
detailed  plans  beyond  these  points. 

Another  method  for  reducing  the  size  of  the  tactical  planning  problem 
is  to  "schedule  Oie  larger  items  of  work  and  expect  the  support  effort-like 
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scaffolding,  material  transfer,  lighting,  preparation  of  welding  and 
cutting  machines  to  follow  along"  (Corporate-Tech  Planning,  Inc.,  1978, 


p.  C-4). 
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IV.  STRATEGIC  SHIPYARD  PLANNING 
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Model  Formulation 

At  the  strategic  level,  shipyard  management  is  concerned  with:  (1) 
determining  target  completion  dates  for  the  various  ships  under  construc¬ 
tion;  and  (2)  specifying  how  the  workforce  levels  for  each  trade  should 
be  changed  from  period  to  period. 

Strategic  planning  is  also  concerned  with  determining  the  most 
appropriate:  (1)  financial  structure;  (2)  manufacturing  technology; 

(3)  portfolio  of  assets;  and  (4)  organization  structure.  The  results 
of  these  decisions  determine  the  environment  in  which  the  strategic 
shipyard  plan  must  be  formulated.  As  such  they  can  be  regarded  as 
super-strategic  decisions.  For  the  purposes  of  the  present  discussion 
it  will  be  assumed  that  the  above  types  of  super-strategic  decisions  . 
have  already  been  made. 

None  of  the  multi-resource/multi-project  planning  models  reported 
in  the  open  literature  simultaneously  considers  both  workforce  levels 
and  target  completion  dates  as  decision  variables.  Fendley 
(1968)  assumed  that  resource  levels  were  fixed  and 
considered  target  project  completion  dates  as  decision  variables.  A 
number  of  authors  (Levy,  Thompson,  and  Wiest,  1963;  McGee  and  Markarian, 
1963;  Lambourne,  1963;  Moshman,  Johnson  and  Laresen,  1963;  Wiest,  1967; 
and  Patterson,  1973)  have  developed  multi-resource/multi-project  planning 
models  based  on  the  assumption  that  due  dates  for  the  various  projects 
are  fixed.  Dar-El  and  Tur  (1978)  developed  a  model  which  determines 
resource  levels  for  a  single  project  so  as  to  minimize  the  "weighted 
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deviation  squares  between  future  resource  loads  and  currently  assigned 
levels".  In  this  model  they  assumed  that  the  resources  required  to 
perform  each  of  the  activities  was  fixed.  Penalties  for  splitting 
activities  and  for  not  scheduling  critical  activities  can  be  added  to 
the  weighted  deviation  squares.  These  penalty  charges  can  be  varied 
so  as  to  generate  a  set  of  solutions  corresponding  to  different  project 
durations.  Bit  level  storage  was  used  to  generate  the  set  of  feasible 
activity  combinations  for  each  day.  In  FORTRAN  the  unit  of  storage  is 
a  computer  word  Consisting  of  32  bits  for  an  IBM  360.  Thus,  the  use  of 
bit  level  storage  makes  it  possible  to  reduce  storage  requirements  by  a 
factor  of  32  for  an  IBM  360.  It  will  also  result  in  a  substantial 
reduction  in  computer  effort  since  logic  operations  require  less  time 
at  the  bit  level  than  at  word  level.  Nevertheless  these  computational 
refinements  will  not  be  sufficient  to  allow  their  algorithm  to  be 
generalized  to  handle  multi-resource/multi-project  problems  in  which 
both  target  completion  dates  and  resource  levels  are  decision  variables. 

Clearly,  the  state-of-the-art  as  reported  in  the  literature  indicate 
that  a  new  model  be  developed  for  solving  mrmp  planning  problems.  A 
prerequisite  for  developing  a  model  for  the  strategic  shipyard  planning 
problem  is  a  clear  understanding  of  how  workforce  levels  and  project 
completion  dates  impact  cost. 

Note  that  the  total  cost  of  the  portfolio  of  ships  under  construction 
is  not  influenced  by  the  completion  date  of  a  given  ship  up  to  the  time 
at  which  the  completion  date  exceeds  the  due  date.  After  the  due  date, 
cost  increases  due  to  the  lateness  penalty  charges  for  the  ship  which  is 
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late.  However,  being  late  on  one  project  can  create  a  series  of  events 
which  can  cause  the  cost  of  other  projects  to  increase.  For  example, 
being  late  on  one  ship  might  delay  the  start  of  a  number  of  other  ships. 
This  could  increase  costs  in  several  ways.  Lateness  penalty  charges 
could  be  incurred  on  subsequent  ships  which  would  not  have  been  incurred 
otherwise.  However,  it  might  be  possible  to  reduce  these  lateness  penalty 
charges  by  temporarily  increasing  the  size  of  the  workforce,  but  this 
will  cause  certain  additional  costs  to  be  incurred.  When  the  workforce 
is  expanded  it  will  be  necessary  to:  (1)  recruit;  (2)  interview;  (3) 
test;  (4)  examine;  (5)  place;  and  (6)  train  the  new  employees.  Reducing 
the  size  of  the  workforce  can  involve:  (1)  exit  interviews;  (2) 
separation  payments;  (3)  increased  premiums  for  employment  insurance; 
and  (4)  a  bad  image  in  the  local  labor  market.  The  manner  in  which 
the  size  of  the  workforce  influences  total  cost  for  a  portfolio  when  all 
other  decision  variables  are  held  constant,  is  illustrated  by  Figure  1  . 

This  figure  shows  that  total  costs  increase  more  sharply  as  the  size 
of  the  workforce  is  decreased  than  for  the  case  in  which  it  is  increased. 
However,  this  will  not  always  be  the  case. 

Wage  rates  for  each  of  the  trades  and  unit  costs  for  the  various 
materials  will  tend  to  increase  abruptly  at  discrete  points  in  time  due 
to  inflation.  This  is  illustrated  by  the  graph  in  Figure  2  which  shows 
how  the  wage  rate  for  a  given  trade  changes  with  the  passage  of  time. 

Both  the  points  in  time  at  which  such  increases  occur  and  their 
amounts  can  often  be  determined  from  the  multi-year  union  contracts.  A 
similar  type  step  function  will  typically  exist  for  each  of  the  trades 
and  each  of  the  raw  materials.  For  many  raw  materials  price  increases 
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Relevant 

Cost 


Change  in  size  of  workforce 


FIGURE 


1.  Influence  of  change  in  workforce  size  to  total  relevant 
cost.  • 
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tend  to  be  announced  at  specified  times  during  the  year.  Forecasting  the 
size  of  such  increases  will  require  an  analysis  of  the  factors  which 
influence  both  the  supply  and  demand  for  the  material  in  question. 
Inflation  in  wage  rates  for  the  various  trades  and  unit  cost  of  raw 
materials  cause  total  cost  to  increase  as  the  duration  of  a  particular 
ship  construction  project  increases.  This  is  depicted  by  Figure  3. 

The  time  required  to  complete  a  given  ship  will  depend  on  the  number 
of  workers  for  each  trade  assigned  to  work  on  it  during  each  time  period. 
This  is  illustrated  by  Figure  4  which  shows  two  hypothetical  percentage 
completion  curves  for  constructing  a  given  ship. 

The  most  obvious  way  in  which  completion  dates  influence  cost  is 
through  late  delivery  penalty  charges.  The  penalty  cost  for  late 
delivery  consists  of:  (1)  penalty  payments  for  late  delivery  explicitly 
specified  in  the  contract;  and  (2)  opportunity  costs  of  the  shipyard 
facilities  tied  up  by  the  late  delivery.  In  theory,  bonuses  for  early 
delivery  are  possible.  However,  they  are  rarely  encountered  in  practice 
and  will  not  be  considered  as  part  of  the  model.  Figure  5  illustrates 
how  the  total  cost  of  a  portfolio  of  projects  behaves  as  the  completion 
date  of  a  single  ship  is  varied  given  that  all  other  variables  which 
influence  total  costs  are  held  constant. 

The  strategic  shipyard  planning  problem  is  to  determine  the  target 
completion  dates  and  workforce  levels  for  each  trade  during  each  period 
in  the  planning  horizon  such  that  total  relevant  cost  will  be  minimized. 
The  following  notation  will  be  used  in  formulating  a  mathematical  model 
of  this  problem: 
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FIGURE  3. 


Duration 
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Inflation  in  wage  rates  and  cost  of  raw  materials  increases 
total  cost  of  a  portfolio  of  ship  construction  projects  as 
the  duration  of  a  project  increases. 
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FIGURE  4  . 
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Start  on  time 


Impact  of  start  date  and  labor  application  rates  on 
completion  time.  • 
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FIGURE  5  . 


Completion  Date 


Impact  of  completion  date  of  a  given  ship  on  total  cost 
of  a  portfolio  ship  construction  project  given  that  all 
other  decision  variables  which  influence  total  cost  are 
held  constant. 
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i  =  index  for  trades 


I 

j 

J 

labor 


material 


hiring 


laying  off 


index  set  of  trades,  i.e.,  I 


{i} 


index  for  ships 


index  set  of  ships. 


i.e. ,  J  =  '  {j} 


labor  hours  of  trade  i  assigned  to  ship  j 
during  period  t 


=  labor  hours  of  trade  i  available  during 
period  t 

=  total  number  of  estimated  labor  hours  of 
trade  i  required  on  ship  j 


unit  cost  per  labor  hour  of  trade  i  for 
period  t 


the  amount  of  material  required  per  labor 
hour  of  trade  i  assigned  to  ship  j 

unit  cost  of  material  associated  with 
craft  i  for  period  t 


=  labor  hours  of  craft  i  hired  during  time 
period  t 

=  upper  limit  on  the  number  of  labor  hours  of 
craft  i  hired  during  any  period 

=  cost  of  hiring  a  labor  hour  of  craft  i 


=  labor  hours  of  trade  i  layed  off  during 
time  period  t 

=  upper  limit  on  the  number  of  labor  hours  of 
trade  i  layed  off  in  any  period 

=  cost  of  laying  off  a  labor  hour  of  trade  i 
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=  start  date  for  ship  j 
Dj  =  due  date  for  ship  j 


T.  =  delivery  date  for  ship  j 

timeliness/ 

lateness 


=  cost  of  lateness  for  ship  j  in  period  t 


P  =  number  of  periods  in  the  planning  horizon 


0  if  construction  on  ship  j  is  not  permitted  during 
period  t 

1  if  construction  on  ship  j  is  permitted  during 
period  t 


In  formulating  this  model  it  will  be  assumed  that: 

(1)  all  hiring  and  laying  off  is  done  at  the  beginning  of  each  time 

period; 

(2)  shipyard  facilities  are  adequate  and  that  lack  of  such  facilities 
will  not  constitute  a  major  source  of  delay; 

(3)  both  intertrade  and  intratrade  interference  are  negligible; 

(4)  overtime  is  not  permitted; 

(5)  strict  precedence  relationships,  which  specify  that  activity  B 

cannot  be  started  until  activity  A  has  been  completed,  are  not 

necessary  in  view  of  the  resulting  long  time  periods  used  in 
strategic  planning. 

The  last  three  assumptions  are  made  to  reduce  the  size  and  complexity  of 

the  problem  so  that  it  will  be  solvable.  However,  these ^factors  will  be 

explicitly  considered  in  the  tactical  shipyard  planning  model. 

Total  relevant  costs  consists  of  the  sum  of:  (1)  cost  of  late 


deliveries;  (2)  labor  costs;  (3)  materials  costs;  (4)  hiring  costs;  and 
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(5)  layoff  costs. 

The  cost  of  late  deliveries  is  given  by 


1  1  L-<-  z-. 
t=D.+1  J 


If  ship  j  is  not  completed  until  period  D^+n  then  penalty  costs  wi] 1  be 
incurred  in  periods  D^+l  t0  D^+n.  This  results  from  the  fact  that  the 
binary  decision  variable  z  was  defined  to  be  equal  to  one  for  the 
case  in  which  construction  on  ship  j  is  permitted  during  period  t. 

Labor  costs  are  given  by  the  following  expression: 


Labor  costs  =  I  E  c.  y. 

_  ,  •  T  it  it 

t=l  id 


Materials  costs  are  given  by 


Materials  costs  =  E  E  E  c.  k. .  x. . 

.  T  •  t  _  it  1J  lit 

iel  icJ  t=s .  J  J 

3 

This  equation  is  based  on  the  assumption  that  trade  i  consumes 


materials  at  the  rate  of  k„  units  of  material  per  labor  hour  working 
on  ship  j . 
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The  cost  of  hiring  additional  workers  is  given  by 


1 

Hiring  costs  =  I  c.  h  (4) 

t=l  1  11 

while  the  cost  of  laying  off  workers  is  given  by 

P  * 

Layoff  costs  =  I  c.  i.  (5) 

t-1  1  1C 

Thus  the  objective  is  to  minimize  the  sum  of  the  quantities  in  equations 
1  through  5  subject  to  a  number  of  constraints  which  must  be  satisfied 
in  order  for  a  solution  to  the  problem  to  be  valid. 

P 

2  x...  =  w. .  for  all  i,  j  (6) 

t=l  iJt  ij 

This  constraint  states  that  the  total  number  of  hours  of  trade  i  assigned 
ship  j  during  the  periods  in  the  planning  horizon  must  be  equal  to  the 
total  number  of  hours  of  trade  i  needed  to  construct  ship  j. 

1  ^4.  zAt  _5  for  all  i,  t  (7) 

jeJ  J  J 

This  constraint  states  that  the  total  number  of  hours  of  trade  i  assigned 
to  each  of  the  J  ships  under  construction  during  time  period  t  must  not 
exceed  the  number  of  hours  available  in  the  trade  i  workforce. 


for  all  i,  t 


(8) 


yit  =  yit-l  +  hit  -  *it 


This  constraint  states  that  the  size  of  craft  i  workforce  at  time  t  is 
equal  to  its  size  during  the  preceeding  period  adjusted  for  the  number 
of  workers  hired  and  laid  off  during  that  period. 


h.  <  u.  for  all  i,  t 

it  —  1 


This  constraint  states  that  there  is  a  limit  to  the  number  of  workers  in 
trade  i  that  can  be  hired  during  time  period  t. 


*it  -<Ui 


for  all  i,  t 


This  constraint  states  that  there  is  a  limit  to  the  number  of  workers  in 
trade  i  that  can  be  laidoff  during  time  period  t. 


T  T 

I  x...  <r  E  x  for  T=l,  2 . p  (11) 

t-i  ^  ^  t-i  gjt 

This  constraint  states  that  the  cumulative  number  of  hours  of  trade  i 
assigned  to  ship  j  cannot  exceed  the  cumulative  number  of  hours  of  trade 

g  assigned  to  ship  j  multiplied  by  a  factor  r.  .  which  represents  the 

J  • 

ratio  of  total  number  of  hours  required  for  trade  i  to  the  total  number 


of  hours  required  for  trade  g  on  ship  j.  This  constraint  is  needed  for 
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only  those  situations  in  which  the  cumulative  progress  of  one  trade 
restricts  the  cumulative  progress  of  another  trade.  For  example,  such 
a  constraint  would  be  needed  to  reflect  the  fact  that  the  number  of  stee] 
plates  welded  into  place  can't  exceed  the  number  of  plates 
that  have  been  cut.  In  addition  to  the  above  constraints  the  decision 
variables  are  restricted  to  the  ranges  defined  below. 


x .  .  >  0,  y .  > 

ijt  —  *it  — 


(0  or  1) 


(12) 


The  model  as  formulated  above  limits  attention  to  a  finite  number 
of  future  time  periods.  The  typical  shipyard  will  continue  to  exist 
well  beyond  this  time.  Thus,  it  is  necessary  to  ensure  that  the  sizes 
of  the  work  forces  for  the  various  trades  be  consistent  with  operations 
beyond  that  point. 

If  periods  beyond  the  planning  horizon  are  ignored,  then  the  model 
as  formulated  could  recommend  dismissing  all  employees  at  the  end  of  the 
recommended  horizon.  Fortunately,  this  problem  can  be  handled  by 
specifying  a  minimum  size  for  the  workforce  for  each  craft  at  the  end 
of  the  planning  horizon. 

It  is  possible  that  not  all  ships  in  the  shipyards  order  backlog 
can  be  completed  by  the  end  of  the  planning  horizon.  Whgn  this  arises 
it  will  be  necessary  to  specify  the  amounts  of  work  which  must  be  done 
on  each  ship  during  the  horizon. 
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The  above  model  Is  a  mixed  integer  quadratic  programming  problem. 


A  procedure  for  solving  this  problem  is  described  in  the  following 
section. 

Solution  Strategy 

To  date  there  have  been  no  algorithms  developed  for  solving  the 
mixed  integer  quadratic  programming  problem  that  was  formulated  in  the 
preceeding  section.  A  strategy  which  can  be  used  to  solve  this 
problem  consists  of  using  the  projection  technique  (Geoffrion,  1970)  to 
transform  the  mixed  integer  quadratic  programming  problem  into  a  linear 
programming  problem  which  will  be  referred  to  as  the  inner  minimization 
problem.  This  will  be  accomplished  by  temporarily  fixing  the  binary 
decision  variables  according  to  a  specified  pattern.  These  patterns 
will  be  determined  in  a  sequential  fashion  according  to  the  method  of- 
steepest  decent.  This  process  will  be  referred  to  as  the  outer 
minimization  problem. 

Commercial  linear  programming  codes  appear  to  be  able  to  handle 
real  world  size  strategic  shipyard  planning  problems.  Lasdon  (1978) 
has  noted  that  such  codes  are  capable  of  handling  problems  with  8,000 
to  16,000  rows  provided  sufficient  core  storage  is  available.  He  also 
notes  that  even  larger  problems  with  a  generalized  upper  bounding  struc¬ 
ture  can  be  solved  and  notes  that  Hirschfeld  (1972)  has  reported  solving 
a  50,000  row  problem  with  generalized  upper  bounding  structure.  This  is 
most  encouraging  since  the  strategic  shipyard  planning  problem  has  a 
generalized  upper  bounding  structure. 


The  above  linear  programming  problem  will  have  to  be  solved  a 
number  of  times  since  it  is  the  inner  maximization  problem.  However, 


the  nature  of  the  strategic  shipyard  planning  problem  is  such  that  com¬ 
position  of  the  optimal  basis  for  the  linear  programming  problem  should 
not  change  appreciably  from  one  search  iteration  to  the  next.  This 
situation  can  be  exploited  to  reduce  computation  time  by  utilizing  the 
capacity  of  commerical  linear  programming  codes  to  save  a  previous  opti¬ 
mal  basis  so  as  to  provide  an  advanced  starting  solution  for  a  subsequent 
problem  (Orchard-Hays,  1974). 
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V.  TACTICAL  SHIPYARD  PLANNING 


Introduction 

The  strategic  shipyard  plan  specifies:  (1)  target  completion 
dates  to  the  various  ships  on  the  shipyard's  order  book  or  currently 
under  construction;  and  (2)  the  size  of  each  trade's  workforce  for 
each  period  during  the  planning  horizon.  These  target  completion  dates 
and  workforce  levels  represent  respectively  the  goals  for  tactical 
management  and  the  resources  which  they  can  use  in  pursuing  these 
goals. 

The  tactical  planners  task  is  to  develop  plans  which  specify  for 
each  period  in  the  tactical  planning  horizon  how  each  trade’s  workforce 
should  be  allocated  to  the  various  activities  necessary  for  constructing 
each  of  the  ships.  In  developing  these  plans  it  will  be  necessary  for 
the  tactical  planner  to  recognize  that:  (1)  both  intertrade  and  intra¬ 
trade  interference  can  cause  substantial  productivity  losses;  (2) 
failure  to  consider  strict  precedence  relationships  can  cause  unnecessary 
delays  and  costs  to  be  incurred;  and  (3)  overtime  can  be  used  when  it  is 
cost  effective  to  do  so.  Recall  that  the  formulation  of  the  strategic 
shipyard  planning  problem  assumed  that  it  was  not  necessary  to  consider 
these  factors  at  the  strategic  planning  level.  This  was  done  to  reduce 
the  size  and  complexity  of  the  problem  so  that  it  could  be  solved. 

The  tactical  planning  process  is  further  complicated  by  random 
variation  in  activity  duration  times  and  uncertainty  regarding  sets  of 
activities  which  will  actually  be  necessary.  These  factors  make  it 
necessary  to  revise  tactical  shipyard  plans  from  time  to  time.  The 
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combined  effect  of  the  above  factors  imply  that  the  tactical  planning 
problem  is  enormously  complicated. 


The  complexity  and  importance  of  the  tactical  shipyard  planning 
problem  suggests  that  it  would  be  desirable  to  have  an  efficient  method 
for  solving  such  problems.  The  purpose  of  this  section  is  to  describe 
the  development  of  such  a  method.  The  development  of  such  a  procedure 
is  accomplished  in  four  phases.  The  first  phase  specifies  a  model 
which  describes  how  the  variables  under  the  control  of  the  tactical 
planner  influence  cost.  The  second  phase  analyzes  and  evaluates  exist¬ 
ing  multi-resource/multi-project  (mrmp)  planning  models.  The  third 
phase  develops  a  procedure  for  solving  the  problem.  The  fourth  phase 
programs  the  solution  procedure. 

Phase  1:  Model  Specification 

In  order  to  formulate  a  model  for  solving  the  tactical  shipyard 
planning  problem  it  is  necessary  to  identify  how  each  of  the  variables 
under  the  control  of  the  tactical  planner  influences  cost. 

If  a  project  is  not  completed  on  time,  then  a  lateness  penalty 
cost  will  be  incurred.  This  cost  consists  of  the  opportunity  cost  of 
having  the  shipyard's  productive  facilities  tied  up  and  additional 
payments  of  late  delivery  explicitly  specified  by  the  contract.  In 
addition,  failure  to  meet  delivery  dates  can  jeopardize  the  shipyard's 
ability  to  get  future  business. 

The  tactical  planner  can  influence  the  amount  of  lateness  penalty 
cost  which  a  shipyard  incurs  through  the  choice  of  the  times  at  which 
activities  which  have  a  high  probability  of  being  on  a  critical  path 
are  started.  For  example,  suppose  that  each  activity  on  the  critical 


path  for  a  given  ship  is  started  at  its  appropriate  time  based  on  the 
assumption  that  each  of  the  activities  will  be  completed  in  the  expected 
time  for  that  activity.  However,  the  actual  time  required  to  complete 
each  activity  will  tend  to  vary  about  its  mean.  This  variation  in  the 
actual  critical  path  activity  duration  time  will  cause  the  project 
duration  to  vary.  As  a  result  of  the  central  limit  theorem  the 
probability  density  function  for  project  duration  times  will  tend  to 
be  normally  distributed  with  mean  and  variance  respectively  equal  to 
the  sum  of  the  means  and  variances  for  each  of  the  activities  on  the 
critical  path.  Since  the  normal  distribution  is  symmetrical  this 
implies  that  starting  each  activity  on  the  critical  path  at  its  late 
start  date  will  result  in  a  probability  of  .5  that  the  project  will  be 
completed  on  time  (this  assumes  that  other  "near  critical"  paths  do  not 
exist).  This  is  illustrated  in  Figure  6  which  shows  that  if  the 
expected  completion  date  for  a  project  is  equal  to  its  due  date,  then 
there  will  be  a  50  percent  chance  that  the  project  will  not  be  completed 
on  time. 

The  tactical  planner  can  reduce  the  odds  that  the  project  will  be 
late  by  starting  some  or  perhaps  all  of  the  activities  on  the  critical 
path  earlier  than  their  respective  start  times  (for  example,  by 
increasing  resource  levels) .  This  will  cause  the  expected  completion 
date  to  be  less  than  the  due  date  which  corresponds  to  shifting  the 
probability  density  function  for  project  completion  times  to  the  left. 
This  will  in  turn  reduce  the  probability  of  being  late. 
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Due  Date  =  E (Completion  date) 


Time 


If  expected  completion  date  equals  due  date,  then 
probability  of  being  late  is  equal  to  0.5. 


I 


Starting  the  critical  path  activities  earlier  than  their  scheduled 
start  times  can  be  visualized  as  insurance  against  lateness.  The  cost 
of  resources  can  be  regarded  as  the  premium  on  the  insurance  policies 
for  the  activities.  The  tactical  planner's  task  is  to  determine:  (1) 
how  much  to  spend  on  insurance;  and  (2)  the  best  portfolio  of  policies 
to  purchase. 

Overtime  can  be  used  to  reduce  the  amount  of  lateness  penalty 
charges  incurred  when  a  project  is  not  completed  on  time,  but  the  use 
of  overtime  involves  overtime  premium  payments.  If  there  are  two  or 
more  activities  on  which  overtime  can  be  used  to  reduce  the  amount  of 
lateness  and  these  activities  both  use  one  or  more  trades  which  are  in 
short  supply,  then  the  problem  of  determining  the  optimal  allocation  of 
overtime, to  the  activities  arise.  The  tradeoff  between  lateness  and 
overtime, assuming  that  all  other  cost  influencing  factors  are  held 
constant  is  depicted  by  Figure  7 . 

In  shipbuilding  there  are  certain  situations  in  which  it  is  not 
possible  to  simultaneously  perform  two  or  more  activities  as  efficiently 
as  the  activities  could  be  performed  separately.  This  can  arise  when  two 
trades  interfere  with  the  work  of  one  another.  This  phenomenon  is 
referred  to  as  intertrade  interference.  For  example,  if  electrical 
cables,  pipes,  and  ventilation  ducts  all  pass  through  a  confined  space 
then  electricians,  pipefitters,  and  vent  installers  could  impede  each 
other's  progress.  It  is  also  possible  for  members  of  the  same  trade  to 
interfere  with  each  other's  work.  This  phenomenon  is  referred  to  as 
intratrade  interference.  For  example,  assigning  too  many  floor  coverers 
to  a  given  area  could  cause  each  to  get  in  the  other's  way. 
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Cost 


FIGURE  7. 


Tradeoff  between  lateness  penalty  cost  and  overtime  cost 
assuming  that  all  other  cost  influencing  facjtors  are  held 
constant. 
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The  duration  of  a  project  can  be  reduced  by  judiciously  utilizing 
intertrade  and/or  intratrade  interference.  This  is  illustrated  by 
Figure  8.  However,  if  the  use  of  intertrade  and/or  intratrade  inter¬ 
ference  is  carried  to  the  extreme,  then  the  project  duration  will 
increase.  This  is  illustrated  by  Figure  9  where  the  minimum  represents 
the  point  where  the  incremental  time  reduction  which  results  from 
performing  certain  activities  in  parallel  and/or  using  more  workers  on 
a  job  is  exactly  offset  by  the  interference  which  results  from  such. 

Some  activities  can  be  scheduled  to  start  at  anytime  while  others 
cannot  be  started  until  a  set  of  prerequisite  activities  have  been 
satisfied.  These  prerequisites  can  be  classified  as  mandatory  or 
preferred.  A  given  activity  cannot  be  accomplished  until  all  of  its 
mandatory  prerequisites  have  been  met.  Failure  to  satisfy  all  preferred 
prerequisites  will  not  block  the  accomplishment  of  the  activity,  but  • 
will  increase  the  cost  of  accomplishing  it.  For  example,  it  might  be 
desirable  to  paint  compartments  prior  to  installing  floor  coverings  and 
fixtures  since  this  will  eliminate  the  need  for  covering  floors  and 
fixtures  with  drop  cloths.  Thus,  another  option  for  reducing  project 
lateness  is  to  selectively  disregard  preferred  prerequisites. 

Random  variations  in  delivery  date  leadtimes  can  cause  material 
shortages  to  arise.  Such  shortages  can  delay  the  start  and/or  completion 
of  certain  of  certain  activities.  When  an  activity  on  the  critical  path 
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(c)  intertrade 

interference 


FIGURE  8 .  Illustration  of  how  project  duration  can  be  reduced  by 
utilizing  intertrade  and  intratrade  interf etence . 
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FIGURE  9  .  Relationship  between  the  project  duration  and  the  degree 
of  interference  allowed. 


is  delayed,  then  project  lateness  can  result.  However,  this  need  not 
necessarily  be  the  case,  since  there  are  a  variety  of  actions  which  can 
be  taken  to  make  up  for  the  time  loss  caused  by  the  delay.  For 
example,  the  following  methods  for  reducing  project  duration:  (1)  over¬ 
time,  (2)  intertrade  interference,  (3)  intratrade  interference,  and  (4) 
activity  splitting;  can  be  used  more  extensively  than  called  for  in  the 
predelay  plan. 

The  risk  of  material  shortages  can  be  reduced  by  allowing  for 
longer  delivery  leadtimes  in  purchasing  raw  materials.  However,  use 
of  longer  delivery  leadtimes  will  result  in  increased  inventory  carry¬ 
ing  cost.  The  tradeoff  between  the  additional  leadtimes  for  the  case 
in  which  all  other  cost  influencing  factors  are  held  constant  is  shown 
in  Figure  10. 

Another  option  for  reducing  the  amount  of  project  lateness  is  to. 
shift  workers  from  activities  whose  completion  can  be  postponed  without 
causing  their  project  to  be  completed  late  to  a  more  urgent  activity. 
This  situation  is  referred  to  as  activity  splitting.  However,  this  will 
typically  cause  set-up  cost  and  transportation  costs  to  increase. 

Figure  11  depicts  the  tradeoff  between  the  number  of  activities 

split  and  project  lateness  for  the  case  in  which  all  other  cost 
influencing  factors  are  held  equal. 

The  above  discussion  suggests  that  the  tactical  planner  can 
influence  the  following  costs:  (1)  cost  of  productivity  lost  as  a 
result  of  splitting  activities;  (2)  cost  of  overtime  premiums;  (3) 
cost  of  productivity  lost  as  a  result  of  intertrade  interference;  (A) 
cost  of  productivity  lost  as  a  result  of  intratrade  interference;  (5) 


Number  of  activities  split 


FIGURE  11  .  Tradeoff  between  the  number  of  activities  split  and 
project  lateness  in  which  all  other  cost  influencing 
factors  are  held  equal. 
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cost  of  project  lateness;  (6)  cost  of  carrying  additional  material 
inventories  to  protect  against  late  deliveries;  and  (7)  cost  of 
disruption  as  a  result  of  late  deliveries  of  key  materials.  Thus 
the  tactical  shipyard  planning  problem  is  to  assign  each  trade's  work¬ 
force  to  the  various  activities  in  a  manner  which  will  minimize  the 
sum  of  the  above  costs  subject  to  the  following  constraints:  (1) 
resource  availabilities  cannot  be  exceeded;  and  (2)  no  activity  can 
be  started  until  all  of  its  mandatory  precedent  requirements  have  been 
met.  However,  an  idle  resource  represents  an  opportunity  loss.  Therefore, 
the  cost  of  idle  time  associated  with  tactical  shipyard  plans  should  be 
calculated  and  reported  to  the  strategic  shipyard  planner.  The 
strategic  planner  will  then  compare  the  cost  of  laying  off  and  rehiring 
of  workers  necessary  to  eliminate  the  idle  time  with  the  savings  in 
labor  cost  and  revise  the  strategic  shipyard  plan  if  necessary. 

Phase  2:  Literature  Review 

The  next  step  in  developing  a  solution  procedure  is  to  review  the 
solution  strategies  used  in  solving  existing  mrmp  planning  models. 

These  procedures  are  found  to  be  either  analytic  or  heuristic. 

Only  one  analytic  approach  has  been  published  to  date.  Pritsker, 
Watters,  and  Wolfe  (1969)  formulated  the  problem  as  a  zero-one  integer 
linear  programming  problem.  Zero-one  variables  are  used  to  indicate 
whether  or  not  an  activity  is  completed  during  each  period  in  the 
planning  horizon.  If  it  has  been  completed,  then  its  successors  may 
start.  If  not,  then  its  successors  cannot  start.  Unfortunately,  such 
an  approach,  if  applied  to  problems  of  the  size  faced  by  shipyards, 
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would  result  in  a  problem  the  size  of  which  would  far  exceed  the 
capacity  of  present  zero-one  integer  linear  programming  algorithms 
(Wiest  and  Levy,  1977).  In  fact,  Elmaghraby  (1977,  p.  217)  reported 
that  Lenstra  (1976)  has  shown  that: 

The  problem  of  scheduling  activities  on  multiple,  resources 
when  the  activities  are  subject  to  precedence  constraints 
and  the  availability  of  the  resources  is  limited  is  known 
to  be  "NP  hard"....  This  implies  that,  in  all  probability, 
there  shall  be  no  "efficient"  algorithm  for  solving  this 
problem  (in  the  sense  of  achieving  optimality) . 

A  number  of  authors  (Fendley,  1968;  Gonguet,  1969;  Pascoe,  1964;  and 
Patterson,  1973)  have  reported  simulation  experiments  that  were  designed 
to  evaluate  the  effectiveness  of  a  number  of  scheduling  rules.  The 
effectiveness  of  these  scheduling  rules  were  evaluated  in  terms  of 
of  the  following  characteristics:  (1)  on  time  completion;  and  (2) 
efficient  resource  utilization.  None  of  the  scheduling  rules  were 
found  to  be  best  for  all  performance  measures.  In  general,  due  date 
oriented  scheduling  rules  performed  better  than  resource  oriented  rules 
when  the  performance  was  measured  in  terms  of  on  time  completion.  The 
reverse  was  true,  in  general,  when  performance  was  measured  in  terms 
of  efficient  resource  utilization. 

None  of  the  heuristic  approaches  to  the  mrmp  scheduling  problems 
with  the  exception  of  a  model  developed  by  Dar-El,  Behmoaram,  and  Tur 
(1978)  utilized  a  cost  based  objective  function.  The  heuristic  used 
in  Levy,  Thompson  and  Wiest  (1963)  was  designed  to  reduce  peak 
resource  requirements  and  smooth  out  period-to-period  assignments 
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subject  to  a  constraint  on  project  duration.  SPAR-1  (Wiest,  1967) 
utilizes  a  heuristic  which  places  primary  emphasis  on  completing  the 
job  as  early  as  possible.  RAMPS  which  was  developed  for  proprietary 
use  utilizes  a  heuristic  which  involves  minimizing  a  weighted  function 
of  variables  such  as  total  slack,  idle  resources,  project  deiay  cost, 
and  number  of  successors  to  a  given  job  (Lambourne,  1963;  Moshman, 
Johnson  and  Laresen,  1963;  Wiest,  1963;  and  Wiest,  1969).  Jenett 
(1970)  lists  a  number  of  other  network  analysis  programs  which  are 
commercially  available  but  for  which  the  scheduling  heuristic  are 
kept  secret. 

The  Dar-El,  Behmoaram,  and  Tur  (DBT)  model  was  found  to  provide 
an  excellent  conceptual  foundation  for  developing  a  model  for  solving 
the  tactical  shipyard  planning  problem.  The  objective  function  of 
this  model  is  equal  to  the  sum  of  the  following  terms:  (1)  total  cost 
of  idle  resources;  (2)  total  activity  splitting  penalties;  (3) 
incremental  total  cost  of  overtime  work;  and  (4)  total  project  lateness 
penalties . 

The  following  notation  is  used  to  specify  the  formulas  for  calcula¬ 
ting  each  of  the  above  cost: 

CD  =  calendar  date 

CF(i)  =  expediting  factor  for  resource  type  i 

CL(j)  “  lateness  penalty 

CS(i)  =  unit  cost  for  resource  type  i 

CSLK(j,k)  =  critical  slack  for  activity  k  on  project  j  (a  user 
fixed  parameter  which  can  be  used  to  specify  the 
priority  associated  with  starting  critical  activities 
prior  to  their  late  start  date) 

DUR(j,k)  =  remaining  duration  of  the  activity  k  on  project  j 

HRA(i)  =  higher  resource  availability  level  for  resource  type  i 
LF(j,k)  =  late  finish  date  of  the  activity  k  on  project  j 
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LRA(i)  =  lower  resource  availability  level  for  resource  type  i 
RA(i)  =  number  of  units  available  of  resource  type  i 
RR(i)  =  number  of  units  required  by  resource  type  i 
SPEN(i)  =  splitting  penalty  for  resource  type  i 
TCIR  =  total  cost  for  idle  resources 
TCOT  =  incremental  total  cost  of  overtime  work 
TASP  =  total  activity  splitting  penalties 
TPLP  =  total  project  lateness  penalties 
The  formulas  for  the  constituent  terms  in  the  objective  function  are  as 
follows: 

(1)  total  cost  of  idle  resources 


TCIR  =  Z  [RA  (i)-RR(i) ]  *  CS(i) 
All  i 


(2)  total  activity  splitting  penalties 


TASP  =  Z 

All  i 


SPEN(i)*RR(i)  if  activity 

was  scheduled  in  immediately 
preceeding  period,  but  not 
completed,  and  not  scheduled 
on  current  day 


0  otherwise 

(3)  total  incremental  cost  of  overtime  work 


tcot  =  z  <; 

All  i 


RR(i)-(HRA(i) )  *CS(i)*CF(i)  if 

LRL(i)  <  RR(i)  <  HRA(i) 


otherwise 
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(4)  total  project  lateness  penalties 


TPLP 


(CD-  LF(j,k)-DUR(j,k)-CSLK(j,k|>'LP(j) 
if  CD>LF(j,k)-DUR(j,k)-CSLK(j,k) 


I  I  ^ 
All  j  All  k 


0  otherwise 


Behmoaram  explained  the  need  for  a  purpose  of  critical  slack  as 
follows: 

"An  activity  whose  start  is  delayed  beyond  its  late  start 
date  will  cause  project  tardiness.  A  lateness  penalty 
should  be  paid  if  an  assignable  (candidate)  activity.... 
is  delayed  beyond  its  late  start  date....  Thus  when 
lateness  penalty  applies  the  project  is  already  in  a 
potential  state  of  being  late  and  since  certain  resources 
are  scarce,  there  is  great  likelihood  that  some  project 
delay  will  occur.  To  increase  the  effectiveness  of  the 
lateness  penalty  and  thereby  avoid  project  tardiness  the 
critical  slack  is  introduced....  the  critical  slack  causes 
the  activity  to  be  considered  late  CSLK  time  units  before 
the  activity  is  really  late..,.  The  critical  slack  therefore 
operates  as  a  safety  factor  to  reduce  the  likelihood  of 
project  tardiness".  (Behmoaram,  1976,  p.  27-30) 

However,  he  failed  to  discuss  two  important  questions:  (1)  How  should 
the  value  of  CSLK  be  determined?;  and  (2)  What  form  should  the  penalty 
function  take?. 

The  problem  is  to  allocate  resources  to  the  various  activities  in 
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a  manner  which  will  minimize  total  relevant  cost  subject  to  the  con¬ 
straints  which  specify  that:  (1)  resource  availabilities  cannot  be 
exceeded;  and  (2)  all  precedent  requirements  must  be  met  before  an 
activity  can  be  started.  This  model  assumes  that  the  level  of  each 
resource  availability  is  constant  for  all  time  periods. 

Dar-El  et  al.  (1978)  realized  that  a  heuristic  approach  would  be 
necessary  for  solving  real  world  sized  mrmp  planning  problems  and 
developed  a  heuristic  approach.  This  model  is  called  SCREAM  (SCarce 
REsource  Allocation  Model)  and  is  programmed  in  FORTRAN.  The  basic 
flow  chart  for  SCREAM  is  shown  in  Figure  12.  Note  that  each  block 
in  this  flow  diagram  has  a  number  in  the  upper  right  hand  corner. 

These  numbers  will  be  used  to  identify  the  blocks  when  they  are  discussed 
in  further  detail  below. 

Block  1  reads  in  the  following  data:  (1)  arrival  and  due  dates  of 
every  project;  (2)  lateness  penalty  and  critical  slack  (a  user  fixed 
parameter  which  can  be  used  to  specify  the  priority  associated  with 
starting  critical  activities  prior  to  their  late  start  dates);  (3)  total 
number  of  activities  in  each  project;  (4)  set  of  immediate  followers 
for  each  activity;  (5)  the  levels  of  the  various  resources  required  to 
accomplish  each  activity;  (6)  unit  cost  of  each  resource;  (7)  amount  of 
each  resource  available  ;  (8)  maximum  amount  of  overtime 

permitted  for  each  resource;  (9)  overtime  cost  for  .each 
resource;  and  (10)  splitting  penalties  associated  with  each  resource. 

Block  2  calculates  the  critical  path  completion  date  for  each  project 
and  the  early  start,  early  finish,  late  start,  late  finish,  and  total 


Figure  12.  Basic  Flow  Chart  for  Scream 


Source:  Dar-El  ,  Behmoaram,  and  Tur  (1978) 
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slack  for  each  activity.  These  items  are  stored  in  matrices  for  further 
use  in  the  scheduling  process. 

Block  3  initializes  the  calendar  to  day  1. 

Block  A  schedules  the  activities  based  on  the  following  logic.  At 
the  start  of  each  period  every  activity  associated  with  any  of  the 
projects  to  be  scheduled  must  be  in  one  of  the  following  mutually 
exclusive  sets:  (1)  activities  which  have  been  completed;  (2)  activities 
which  have  been  started  but  have  not  yet  been  completed;  (3)  activities 
not  yet  started  that  have  no  precedent  activities;  (A)  activities  not 
yet  started  that  have  no  precedent  activities  yet  to  be  completed;  and 
(5)  activities  that  cannot  be  started  due  to  the  fact  that  one  or  more 
precedent  activities  have  not  yet  been  completed. 

In  a  given  period,  activities  in  sets  2,  3,  and  A  are  eligible  for 
the  assignment  of  resources.  Such  activities  are  referred  to  as  assign¬ 
able  activities.  However,  resource  availabilities  will  not  usually  be 
sufficient  to  allow  all  assignable  activities  to  be  assigned  simulta¬ 
neously  in  a  given  period.  This  creates  the  problem  of  determining 
which  activities  should  be  assigned  in  any  given  period. 

The  assignments  made  in  a  given  period  determine  the  set  of 
assignable  activities  for  future  periods.  The  minimum  cost  attainable 
in  any  given  period  will  be  dependent  on  the  set  of  assignable  activities 
for  that  period.  This  implies  that  in  order  to  minimize  total  system 
cost  it  will  be  necessary  to  account  for  the  impact  that  the  assignments 
made  in  a  given  period  have  on  the  costs  of  future  periods.  Unfortunately, 
this  will  result  in  a  combinatorial  problem  which  is  far  too  large  to 
be  solved  in  a  cost  effective  manner. 


In  order  to  avoid  this  problem  Behmoaram  utilized  a  strategy  which 
was  first  utilized  by  Weist  (1963).  This  strategy  consists  of  decom¬ 
posing  the  multi-period  optimization  problem  into  a  series  of  sub-problems, 
one  for  each  period  in  the  planning  horizon,  and  solving  these  sub¬ 
problems  independently.  Admittedly,  such  an  approach  cannot  guarantee  to 
minimize  total  system  cost.  Nevertheless  this  appears  to  be  a  resonable 
approach  for  dealing  with  a  problem  which  would  otherwise  be  unsolvable. 

The  sub-problem  for  each  period  is  solved  as  follows:  (1)  a  "binary 
enumeration  technique"  developed  by  Dar-El  and  Tur  (1975)  is  used  to  gen¬ 
erate  all  possible  combinations  of  assignable  activities.  This  binary 
enumeration  technique  uses  bit  level  storage  in  order  to  -exploit  the  fact 
that  a  computer  word  in  FORTRAN  on  the  IBM  360  contains  32  bits.  In 
dition,  the  time  required  to  perform  logic  operations  is  much  less  at  the 
bit  level  than  at  the  word  level;  (2)  the  resource  requirements  for  each 
assignable  activity  combination  is  computed  and  compared  to  the  amounts 
of  the  various  resources  which  are  available  during  the  period  of  interest. 
These  combinations  which  require  one  or  more  of  the  various  resources 
than  is  available  are  classified  as  infeasible  and  discarded.  The  re¬ 
mainder  of  the  combinations  are  feasible;  (3)  the  value  of  the  objective 
function  is  computed  for  each  feasible  combination  of  assignable  activi¬ 
ties;  and  (4)  the  minimum  value  of  the  objective  function  for  the  single 
period  problem  is  determined  by  comparing  each  value  computed  with  the 
previous  minimum  value.  The  new  value  replaces  the  previous  minimum  when 
the  new  value  is  less  than  the  previous  minimum.  Otherwise,  the  new  value 
is  discarded. 


Block  5  is  an  updating  routine  which  is  called  when  the  scheduling 
process  for  a  given  day  has  been  completed.  It  performs  the  following 
operations:  (1)  reduces  the  duration  for  the  activities  scheduled; 

(2)  decreases  the  total  slack  of  unscheduled  activities  by  one  unit; 

(3)  increases  the  calendar  date  by  one  unit;  and  (A)  removes  activities 
that  have  been  completed  from  the  set  of  available  activities. 

Block  6  checks  whether  or  not  all  activities  have  been  scheduled. 
Blocks  7,  8,  and  9  form  a  subsystem  which  permits  the  model  to  be  used 
in  an  interactive  mode  the  need  for  which  is  explained  below 

The  amount  of  time  required  to  accomplish  a  given  activity  will 
be  dependent  on  the  work  method  used  to  complete  it.  For.  example, 
tamden-arc  automatic  welding  equipment  requires  about  25  percent  less 
time  than  single-arc  automatic  welding  (Mack-Forlist  and  Newman,  1970). 
Another  example,  the  amount  of  time  required  to  fabricate  a  panel  will 
be  dependent  on  the  extent  to  which  assembly  line  methods  are  used, 
the  degree  to  which  the  assembly  line  has  been  balanced,  and  the  priority 
rules  which  determine  the  order  in  which  different  types  of  panels  enter 
the  line.  The  use  of  a  particular  work  method  for  a  given  activity  af¬ 
fects  not  only  the  time  required  to  complete  that  activity  but  the  re¬ 
sources  which  will  be  available  to  the  other  activities.  Therefore,  it 
would  appear  that  the  tactical  shipyard  planner  should  consider  alter¬ 
native  work  methods.  One  approach  for  doing  this  has  been  proposed  by 
Wiest  (1967)  in  his  SPAR-1  model.  This  model  allows  the  user  to  specify 
three  crew  sizes:  (1)  a  normal  crew  size  which  represents  the  number  of 
men  or  other  resources  typically  assigned  to  the  job;  (2)  a  maximum  crew 
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size  which  represents  the  number  of  men  required  for  crashing  the  job; 
and  (3)  a  minimum  crew  size  which  represents  the  smallest  number  of  men 
which  can  be  assigned  to  the  job.  Pritsker,  Watters,  and  Wolfe  (1969) 
recognized  that  certain  resources  are  substitutable  and  note  that  such 
situations  can  be  incorporated  into  their  zero-one  integer  linear  pro¬ 
gramming  model  through  the  use  of  multiple  choice  type  constraints. 
Dar-El  and  Tur  (1976)  permit  the  user  to  specify  a  number  of  alternative 
resource  combinations  (ARC’s)  which  can  be  used  to  accomplish  each  ac¬ 
tivity.  Behmoaram  notes  that: 

"Tf  an  absolute  minimum  for  the  system  cost  is  being  sought, 
then  every  ARC  for  all  activities  should  be  considered  in 
the  scheduling  process.  However,  this  would  enormously 
complicate  an  already  complicated  combinatorial  problem  and 

increase  the  computation  time  and  cost  exponentially . " 

(Behmoaram,  1976,  p.  35) 

However,  he  suggests  that  one  way  to  tackle  this  problem  is  "to  build  a 
monitoring  function  into  the  model  which  sets  off  an  alarm,  requiring 
the  scheduler  to  investigate  the  utility  of  substituting  one  or  more 
activity  ARC's". 

The  purpose  of  Block  7  is  to  perform  this  monitoring  function. 

This  is  accomplished  by  comparing  DELTA,  minimum  cost  value  of  the 
objective  function  to  an  upper  control  limit.  An  alarm  is  set  off  when 
the  value  of  DELTA  exceeds  its  upper  control  limit  as  indicated  by 
Block  8.  When  the  alarm  occurs  the  scheduler  should  attempt  to  identify 
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new  ARC's.  The  model  will  be  rerun  with  each  of  these  new  ARC's  and 
the  one  with  the  minimum  value  of  DELTA  will  be  utilized.  This  is  an 
interactive  process  which  is  indicated  by  Block  9. 

The  strategy  of  generating  ARC's  is  likely  to  require  non-trivial 
amounts  of  engineering  time.  Therefore,  it  is  likely  to  be  expensive. 
Therefore,  postponing  the  generating  of  ARC's  until  an  alarm  occurs  is 
appealing  from  the  standpoint  of  saving  engineering  time.  This 
advantage  is  offset  somewhat  by  the  necessity  of  having  to  make  multiple 
runs  which  will  increase  computation  cost.  In  addition  the  problem  of 
how  to  determine  the  value  of  the  upper  control  limit  has  not  been 
answered  satisfactorily. 

Phase  3:  Specification  of  Solution  Procedure 

The  tactical  shipyard  planning  problem  is  to  minimize  the  sum  of 
the  following  cost: 

(1)  cost  of  productivity  lost  as  a  result  of  splitting  activities; 

(2)  cost  of  overtime  premiums; 

(3)  cost  of  project  lateness; 

(4)  cost  of  productivity  lost  as  a  result  of  intertrade  inter¬ 
ference; 

(5)  cost  of  productivity  lost  as  a  result  of  intratrade  inter¬ 
ference; 

(6)  cost  of  disruption  resulting  from  late  deliveries  of  materials 
subject  to  the  following  constraints: 

1.  resource  availabilities  cannot  be  exceeded; 

2.  no  activity  can  be  started  until  all  of  its  precedent 
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requirements  have  been  met. 


In  theory  this  problem  could  be  formulated  as  a  zero-one  integer  linear 
programming  problem.  However,  Lenstra  (1976)  as  reported  by  Elmaghraby 
(1977,  p.  217)  has  shown  that 

The  problem  of  scheduling  activities  on  multiple  resources 
when  the  activities  are  subject  to  precedence  constraints 
and  the  availability  of  resources  is  limited  is  known  to 
be  "NP  hard"....  This  implies  that,  in  all  probability, 
there  shall  be  no  "efficient"  algorithm  for  solving  this 
problem  (in  the  sense  of  achieving  optimality). 

This  implies  that  the  tactical  shipyard  planning  problem  is  "NP  hard". 
This  indicates  that  a  search  for  an  analytic  procedure  for  solving  the 
problem  would  not  likely  be  fruitful.  Thus  it  was  decided  not  to 
consider  analytic  methods  further.  This  left  only  two  options  for 
attacking  the  problem:  (1)  heuristic  procedures;  or  (2)  branch  and 
bound  theory . 

The  branch  and  bound  approach  is  a  systematic  procedure  whereby 
each  point  in  the  entire  space  of  feasible  solutions  is  enumerated 
either  explicitly  or  implicitly.  It  is  desirable  to  minimize  the 
number  of  points  that  have  to  be  enumerated  explicitly.  Three  tech¬ 
niques  are  used  to  accomplish  this:  (1)  dominance;  (2)  feasibility; 

(3)  redundancy.  However,  Baker  (1974,  p.  276-277)  has  noted 
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Unfortunately,  all  implicit  enumeration  approaches  to  the 
determination  of  an  optimal  schedule  appear  to  be  susceptible 
to  the  combinatorial  nature  of  these  problems  when  they  are 
tested  on  the  large  versions  typically  found  in  practice.... 
there  is  no  evidence  that  such  techniques  can  reliably  handle 
multi-resource  versions  of  a  problem  that  contains  more  than 
50  activities. 

Thus  it  was  decided  not  to  consider  branch  and  bound  methods  as  candida¬ 
tes  for  the  solution  procedure.  This  left  heuristic  procedures  as  the 
only  viable  option  for  attacking  the  tactical  shipyard  planning  problem. 

A  review  of  heuristic  models  published  in  the  open  literature 
revealed  that  none  of  the  heuristic  approaches  to  the  mrmp  scheduling 
problems  published  in  the  open  literature  with  the  exceptions  of  a 
model  developed  by  Dar-El,  Behmoaram  and  Tur  (1978)  utilized  a  cost 
based  objective  function.  The  heuristic  use  in  Levy,  Thompson  and  Wiest 
(1963)  was  designed  to  reduce  peak  resource  requirements  and  smooth  out 
period-to-period  assignments  subject  to  a  constrain:  on  project  duration. 
SPAR-1  (Wiest,  1967)  utilizes  a  heuristic  which  places  primary  emphasis 
on  completing  the  job  as  early  as  possible.  RAMPS  which  was  developed 
for  proprietary  use  utilizes  a  heuristic  which  involves  minimizing  a 
weighted  function  of  variables  such  as  total  slack,  idle  resources, 
project  delay  cost,  and  number  of  successors  to  a  given  job  (Lambourne, 
1963;  Moshman,  Johnson  and  Laresen,  1963;  Wiest,  1963;  and  Wiest,  1969). 


Jennett  (1970)  lists  a  number  of  other  network  analysis  programs  which 
are  commercially  available  but  for  which  the  scheduling  heuristics 
are  kept  secret. 

The  model  considered  by  Dar-El,  Behmoaram  and  Tur  (1978)  was  to 
minimize  the  sum  of  the  following  costs. 

(1)  cost  of  productivity  lost  as  a  result  of  splitting  activities; 

(2)  cost  of  overtime  premiums; 

(3)  cost  of  project  lateness;  and 

(4)  total  cost  of  idle  resources; 
subject  to  the  following  constraints: 

(1)  resource  availabilities  cannot  be  exceeded;  and 

(2)  no  activity  can  be  started  until  all  of  its  precedent  require¬ 
ments  have  been  met. 

Note  that  this  problem  differs  from  the  tactical  shipyard  planning 
problem  in  that  it  does  not  consider  the  following  costs: 

(1)  cost  of  intertrade  interference; 

(2)  cost  of  intratrade  interference; 

(3)  cost  of  carrying  additional  material  inventories  to  protect 
against  late  deliveries;  and 

whereas  the  tactical  shipyard  planning  problem  does.  Another  difference 
is  that  the  Dar-El,  Behmoaram  and  Tur  (DBT)  model  includes  the  cost  of 
idle  resources  in  their  objective  function  whereas  the  tactical  shipyard 
planning  problem  does  not.  The  reason  for  not  including  this  cost  in  the 
tactical  shipyard  planning  model  is  that  such  costs  are  primarily  determined 
by  the  size  of  each  trade’s  workforce.  Another  difference  is  that  the  DBT 
model  assumes  that  the  resource  level  is  constant  for  all  time  periods. 

In  the  tactical  shipyard  planning  problem  it  is  possible  to  vary  resource 
levels  from  period  to  period.  These  variables  are  specified  by  strategic 


level  management.  Therefore,  it  is  not  reasonable  to  have  these  costs 
in  the  tactical  shipyard  planning  model.  However,  the  tactical  ship¬ 
yard  planning  model  does  report  the  percentage  of  idle  time  of  each 
resource  to  the  strategic  level  planner  who  evaluates  this  cost  with 
that  of  changing  appropriate  resource  levels.  The  DBT  model  was 
modified  to  handle  this  by  removing  the  cost  of  idle  resources  from  the 
objective  function  and  modifying  the  output  procedure  so  that  this 
variable  could  be  printed  out  separately. 

A  special  procedure  had  to  be  developed  for  handling  intertrade 
interferences.  This  procedure  consisted  of  the  following  steps: 

(1)  put  all  activities  in  the  following  categories  into  a 

scheduling  pool: 

(a)  activities  which  have  been  started  but  have  not  yet  been 
completed; 

(b)  activities  not  yet  started  that  have  no  precedents;  and 

(c)  activities  not  yet  started  for  which  all  precedents  have 
been  completed. 

This  is  illustrated  by  Figure  13  which  graphically  depicts  the  com¬ 
position  of  the  scheduling  pool. 

(2)  partition  the  scheduling  pool  into: 

(a)  an  interference  free  subset  which  contains  those  activities 
which  can  be  performed  with  any  other  member  of  the  subset 
without  creating  interference; 

(b)  a  number  of  mutually  exclusive  interference  subsets  such 
that  the  activities  within  a  given  interference  subset 
interfer  with  one  another,  while  activities  in  different 


started 


Scheduling 

pool 


FIGURE  13  .  Composition  of  scheduling  pool. 


subsets  do  not. 


This  is  illustrated  by  Figure  14  which  depicts  the  partitioning  of  the 
scheduling  pool  into  an  interference  free  subset  and  n  interference 
subsets.  The  activities  in  each  interference  subset  cannot  be  performed 
simultaneously  with  another  member  of  the  same  set  without  creating 
interference.  However,  any  combination  of  activities  consisting  of  all 
activities  in  the  interference  free  subset  and  at  most  one  activity 
from  the  n  interference  can  be  performed  simultaneously  without  creating 
interference. 

(3)  find  all  possible  combinations  of  activities  which  consists  of 
one  member  from  each  of  the  mutually  exclusive  interference 
subsets  and  all  of  the  activities  in  the  interference  free 
subset; 

This  is  illustrated  by  Figure  15  which  shows  all  possible  combinations 
of  activities  which  can  be  formed  when  the  interference  subset  consists 
of  activities  A,  B,  and  C;  the  first  interference  subset  consists  of 
D  and  E;  and  the  second  interference  subset  consists  of  F  and  G. 

(4)  calculate  resources  required  for  each  combination  of  activities, 
compare  resource  requirements  with  resource  availability  levels, 
and  eliminate  infeasible  combinations; 

(5)  calculate  value  of  objective  function  for  each  feasible 
resource  combination  and  keep  combinations  which  minimize 
value  of  objective  function  for  that  period. 

The  logic  for  dealing  with  intertrade  interference  is  to  first 
attempt  to  schedule  the  activities  without  intertrade  interference.  If 
all  projects  can  be  completed  by  their  respective  due  dates  then  the 
interference  free  solution  will  dominate  any  solution  which  has 
interference.  However,  if  lateness  results,  then  the  scheduler  should 
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Possible  Acitivy  Combinations 
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FIGURE  15.  Formation  of  all  possible  interference  free  combinations 
of  activities  where  D  and  E  interfere  and  F  and  G 
interfere. 


selectively  permit  intertrade  interference  to  arise.  This  is  accomplished 
by  removing  the  constraint  that  specifies  that  two  activities  cannot  be 
performed  together  and  increasing  the  times  required  to  perform  the 
activities . 

Intratrade  interference  can  be  handled  by  rerunning  the  model  for 
each  crew  size  which  can  be  used  to  perform  an  activity. 

In  doing  this  it  will  be  necessary  to  adjust  the  activity  duration  times 
to  account  for  the  effects  of  diminishing  returns  which  result  from 
intratrade  interference  and  the  cost  of  performing  the  activity  to 
reflect  the  size  of  the  crew. 

The  cost  of  late  deliveries  of  materials  can  be  handled  by  regarding 
the  arrival  of  a  material  as  an  activity  and  calculating  the  lateness 
penalty  for  that  activity.  In  this  case  the  penalty  function  will  be 
the  expected  cost  of  lateness  caused  by  material  shortages  which  corres¬ 
ponds  to  the  various  values  of  delivery  lead  time. 

An  example  will  be  presented  to  illustrate  the  logic  for  specifying 
the  lateness  penalty  function  and  the  associated  incremental  inventory 
carrying  cost  function. 

EXAMPLE.  The  main  engine  is  to  be  built  by  a  subcontractor.  If 
the  engine  is  not  delivered  in  time  to  be  installed  prior  to  launch, 
then  an  additional  cost  of  installation  equal  to  will  be  incurred. 

Let  t^  represent  the  latest  date  for  delivery  of  the  engine  which  will 
permit  it  to  be  installed  prior  to  launch.  If  the  ship^prd  specifies  a 
delivery  date  of  t^  then  there  will  be  a  50  percent  chance  of  receiving 
delivery  on  time,  33.3  percent  chance  of  being  one  period  late,  and  a 
16.7  percent  chance  of  being  two  periods  late.  The  shipyard  has  to 
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pay  M  dollars  for  the  engine  on  delivery.  If  the  delivery  date  denoted 
by  t^  is  specified  to  equal  t^,  then  there  will  be  a  50  percent  chance 
that  delivery  will  be  at  least  one  period  late.  The  expected  cost  of 
late  delivery  will  be  .5  Cp.  If  t^  is  specified  to  equal  then 

the  expected  cost  of  late  delivery  will  be  .167  Cp.  However,  there 
will  be  a  50  percent  chance  that  the  engine  will  be  delivered  one  period 
early.  If  this  occurs,  then  the  shipyard  will  incur  an  opportunity 
cost  of  iM  dollars  where  i  represents  the  shipyards  cost  of  capital. 

The  expected  opportunity  loss  will  be  .5  iM.  If  t^  is  specified  to 
equal  to  t^  2>  then  there  will  be  a  zero  percent  chance  of  late 
delivery.  However,  there  will  be  a  50  percent  chance  of  receiving 
delivery  two  periods  early  and  a  33.7  percent  chance  of  receiving 
delivery  one  period  early.  The  expected  opportunity  loss  will  be 
equal  to  .5  (2i  +  i^)  M  +  .337  iM. 

The  lateness  penalty  function  and  the  inventory  carrying  cost 
function  for  the  above  situation  is  shown  in  Figure  16. 

Phase  4:  Development  of  Computer  Programs 

A  list  of  an  experimental  computer  program  for  solving  the  tactical 
shipyard  planning  model  is  contained  in  Appendix  A. 


Charges  for  Period 

Incremental 

Delivery  Date  Specified  Lateness  Penalty  Inventory  Carrying  Case 
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FIGURE  16.  Lateness  penalty  function  and  inventory  carrying  cost 
function. 
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APPENDIX  A 


Listing  of  Computer  Program  for  Solving 
the  Tactical  Shipyard  Planning  Problem. 
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